Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,664,363	(G)6→7	coding (256/729 nt)	gadW ←	transcriptional activator of gadA and gadBC, repressor of gadX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,664,357	1	.	G	95.9%	254.1 / 2.6	73	coding (262/729 nt)	gadW	transcriptional activator of gadA and gadBC, repressor of gadX
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (39/31);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

ATTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAT  >  NC_000913/3664223‑3664479                                                                                                                                        |                                                                                                                            aTTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAt                                                                                                                              >  2:1222035/1‑135 (MQ=255) aTTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAt                                                                                                                              <  1:1222035/135‑1 (MQ=255) aTTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑gggg                                                                                                                        >  1:1132949/1‑139 (MQ=255) aTTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑gggg                                                                                                                        >  1:223689/1‑139 (MQ=255) aTTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑gggg                                                                                                                        >  2:48139/1‑139 (MQ=255) aTTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑ggg                                                                                                                         >  1:798480/1‑138 (MQ=255) aTTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑ggg                                                                                                                         <  2:798480/138‑1 (MQ=255)  ttGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑ggggg                                                                                                                       <  2:223689/139‑1 (MQ=255)       aaGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑gggggggtga                                                                                                                  <  1:698735/139‑5 (MQ=255)        aGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑gggggggtgac                                                                                                                 <  1:528216/139‑6 (MQ=255)         gTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACt                                                                                                                 <  2:685123/139‑1 (MQ=255)         gTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACt                                                                                                                 <  1:1102681/139‑1 (MQ=255)            aGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGt                                                                                                              >  1:117965/1‑139 (MQ=255)             ggggATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGtt                                                                                                             <  1:149587/139‑1 (MQ=255)                 aTCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGttaatt                                                                                                         >  1:738404/1‑139 (MQ=255)                 aTCACTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGttaatt                                                                                                         <  1:866961/139‑1 (MQ=255)                   cAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGttaattaa                                                                                                       >  1:1290066/1‑139 (MQ=255)                    aGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAAc                                                                                                      >  2:587315/1‑139 (MQ=255)                    aGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAAc                                                                                                      >  1:432424/1‑139 (MQ=255)                     gTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACa                                                                                                     >  2:121374/1‑139 (MQ=255)                      ttctttctTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACaa                                                                                                    <  2:690008/139‑1 (MQ=255)                       tctttcttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACaaa                                                                                                   <  1:587315/139‑1 (MQ=255)                           tcttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACt                                                                                                                 <  2:1239998/121‑1 (MQ=255)                           tcttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACt                                                                                                                 >  1:1239998/1‑121 (MQ=255)                           tcttATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTCAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGt                                                                                               >  2:869425/1‑139 (MQ=255)                                tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          <  2:760231/139‑1 (MQ=255)                                tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTg                                                                                          >  2:197442/1‑139 (MQ=255)                                tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑ggg                                                                                                                         <  1:1311125/107‑1 (MQ=255)                                tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑ggg                                                                                                                         >  2:1311125/1‑107 (MQ=255)                                tGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGc                                                                                         >  1:574593/1‑139 (MQ=255)                                     aaaCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     >  2:139767/1‑139 (MQ=255)                                     aaaCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     <  2:927115/139‑1 (MQ=255)                                       aCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTa                                                                                   >  2:1290106/1‑139 (MQ=255)                                       aCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTa                                                                                   <  2:1061355/139‑1 (MQ=255)                                                aaCAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAAt                                                                          <  1:197442/139‑1 (MQ=255)                                                   aGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGc                                                                       >  1:468829/1‑139 (MQ=255)                                                     gAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     <  1:1210880/123‑1 (MQ=255)                                                     gAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     >  2:1210880/1‑123 (MQ=255)                                                     gAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑gggggggtga                                                                                                                  >  1:985878/1‑89 (MQ=255)                                                     gAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑gggggggtga                                                                                                                  <  2:985878/93‑5 (MQ=255)                                                          gTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAg                                                                >  2:1240756/1‑139 (MQ=255)                                                          gTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAg                                                                <  1:121374/139‑1 (MQ=255)                                                            agagCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCt                                                              <  2:1132949/139‑1 (MQ=255)                                                                   atTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAat                                                       <  1:1290106/139‑1 (MQ=255)                                                                    ttGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAata                                                      >  1:230878/1‑139 (MQ=255)                                                                    ttGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAata                                                      >  1:1171333/1‑139 (MQ=255)                                                                    ttGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCGAAATTGCGTGGTAGCTGACGAata                                                      <  1:139767/139‑1 (MQ=255)                                                                     tGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACaa                                                                                                    >  1:1031219/1‑92 (MQ=255)                                                                     tGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACaa                                                                                                    <  2:1031219/92‑1 (MQ=255)                                                                     tGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAatat                                                     >  1:803828/1‑139 (MQ=255)                                                                        tCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATAtttt                                                  >  2:930033/1‑139 (MQ=255)                                                                          aGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTg                                                 >  2:347760/1‑139 (MQ=255)                                                                          aGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACaa                                                                                                    >  2:301172/1‑87 (MQ=255)                                                                          aGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACaa                                                                                                    <  1:301172/87‑1 (MQ=255)                                                                          aGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTgt                                                >  2:155672/1‑139 (MQ=255)                                                                          aGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTgt                                                >  1:300815/1‑139 (MQ=255)                                                                            tATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACTAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     >  1:1010743/1‑100 (MQ=255)                                                                            tATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCtgtg                                                                                     <  2:1010743/100‑1 (MQ=255)                                                                            tATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTgtaa                                              >  2:38437/1‑139 (MQ=255)                                                                                 gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAAc                                                                                                  <  1:514587/82‑1 (MQ=255)                                                                                 gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAAc                                                                                                  >  2:514587/1‑82 (MQ=255)                                                                                 gTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAAt                                         >  2:1096936/1‑139 (MQ=255)                                                                                      tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCaa                                    >  1:820568/1‑139 (MQ=255)                                                                                      tttACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCaa                                    >  2:1280297/1‑139 (MQ=255)                                                                                           cccGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGt                               >  2:932656/1‑139 (MQ=255)                                                                                                aTTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCt                          >  2:394703/1‑139 (MQ=255)                                                                                                aTTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTAGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCt                          <  1:869425/139‑1 (MQ=255)                                                                                                aTTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑GGGGGGTGACt                                                                                                                 <  1:1290160/51‑1 (MQ=255)                                                                                                aTTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑GGGGGGTGACt                                                                                                                 >  2:1290160/1‑51 (MQ=255)                                                                                                          gCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACaa                                                                                                    <  2:1039760/55‑1 (MQ=255)                                                                                                          gCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACaa                                                                                                    >  1:1039760/1‑55 (MQ=255)                                                                                                          gCTATCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATc                <  1:1240756/139‑1 (MQ=255)                                                                                                              tCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTACTTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGt                                                                     >  1:879325/1‑82 (MQ=255)                                                                                                              tCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAgt            >  1:598823/1‑139 (MQ=255)                                                                                                              tCAAAAATCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGAAACCCTAAATTGCGt                                                                     <  2:879325/82‑1 (MQ=255)                                                                                                                 aaaaTCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAgtcgt         <  2:738404/139‑1 (MQ=255)                                                                                                                 aaaaTCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAgtcgt         <  2:803828/139‑1 (MQ=255)                                                                                                                 aaaaTCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAgtcgt         <  2:300815/139‑1 (MQ=255)                                                                                                                    aTCGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTcgc      >  2:72092/1‑139 (MQ=255)                                                                                                                      cGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGt    >  2:605069/1‑139 (MQ=255)                                                                                                                      cGCCGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGt    <  2:432424/139‑1 (MQ=255)                                                                                                                        ccGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAt  <  2:1290066/139‑1 (MQ=255)                                                                                                                        ccGTCACCAGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAt  >  2:1069016/1‑139 (MQ=255)                                                                                                                           tCACCAGCGGCATGGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGa                                                         >  2:621954/1‑82 (MQ=255)                                                                                                                           tCACCAGCGGCATGGGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGa                                                         <  1:621954/82‑1 (MQ=255)                                                                                                                                aGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTa                                                                             >  1:226784/1‑56 (MQ=255)                                                                                                                                aGCGGCATG‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTa                                                                             <  2:226784/56‑1 (MQ=255)                                                                                                                                        |                                                                                                                            ATTGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCAT‑‑GGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTAT  >  NC_000913/3664223‑3664479

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑