Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,676,805	0	T	G	39.8%	108.2 / 21.5	83	T367P (ACC→CCC)	pntA	pyridine nucleotide transhydrogenase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (22/28);  new base G (33/0);  total (55/28)
Fisher's exact test for biased strand distribution p-value = 9.61e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGTTTGTGCCGTAAAGCTGTGAGGATTGCGTCGGCAGAC  >  NC_000913/1676656‑1676951                                                                                                                                                       |                                                                                                                                                     cATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑a                                                                                                                                                    <  1:1218148/151‑1 (MQ=255)       gCCAAAAAGAATGCTTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGAACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑CAtttt                                                                                                                                               >  1:1215525/1‑149 (MQ=255)        ccAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCCGGG‑AATTTCg                                                                                                                                              >  2:448722/1‑150 (MQ=255)        ccAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AAttt                                                                                                                                                >  2:291212/1‑148 (MQ=255)         cAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGcc                                                                                                                                            <  2:717328/151‑1 (MQ=255)          aaaaaGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTGCCTCAGTTTTCGCTTCCGGTGCGGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGccc                                                                                                                                           <  2:996487/151‑1 (MQ=255)             aaGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑GCTGGT‑AATTTTGccccc                                                                                                                                         >  1:484150/1‑148 (MQ=255)             aaGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCCG‑CCAGGT‑AATTTCTccc                                                                                                                                           >  1:538012/1‑148 (MQ=255)                 aTGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACg                                                                                                                                      >  1:15281/1‑149 (MQ=255)                 aTGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTGCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGa                                                                                                                                    <  1:726887/151‑1 (MQ=255)                  tGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTCGCCCGCACGGAt                                                                                                                                   >  1:1315172/1‑151 (MQ=255)                    aTTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTTGCCCGCACgg                                                                                                                                     >  2:864579/1‑147 (MQ=255)                       gCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGCGTCCG‑CCAGGG‑AATTTCGCCCCCCCGGGT‑CACgg                                                                                                                             >  2:850382/1‑151 (MQ=255)                        ccAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCCGGG‑AATTTTGCCCGCCCGGGT‑CACgg                                                                                                                             >  2:1207835/1‑150 (MQ=255)                        ccAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTCGCCCGCACGGAT‑CACgg                                                                                                                             >  1:1004650/1‑150 (MQ=255)                         cAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AAATTCGCCCGCCCGGAT‑CACgg                                                                                                                             >  1:284524/1‑149 (MQ=255)                             gcCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTCGCCCGCACGGAT‑CACGGTCAc                                                                                                                         >  1:791275/1‑149 (MQ=255)                               cATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCGGCTTTTTGCGCCGCCTGCCGCTGAGCTGATACCTGAATCCGCGGCGCCGG‑GCCGGG‑AATTTCCCCCCCCCGGAT‑CCCGGTCTcgcc                                                                                                                      >  1:413593/1‑150 (MQ=255)                                aTCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCCGGG‑AATTTCGCCCGCCCGGAT‑CACGGTCAcgccg                                                                                                                     >  2:1011530/1‑150 (MQ=255)                                  cAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCcgcg                                                                                                                   <  1:265814/150‑1 (MQ=255)                                     cgcgTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATc                                                                                                               <  2:1261227/151‑1 (MQ=255)                                     cgcgTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTCGCCCCCACGGAT‑CACGGTCACGCCGCGaa                                                                                                                 >  1:622634/1‑149 (MQ=255)                                          aTTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGTCGG‑GCAGGG‑AATTTTGCCCCCACGGAT‑CAAgggaacccc                                                                                                                      >  2:1179424/1‑132 (MQ=255)                                           tttACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AAATTTACCCCCGCCGAT‑TACAGGGAAACCTCGAATcacca                                                                                                           >  2:801034/1‑149 (MQ=255)                                                gCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCAcatc                                                                                                       >  1:958546/1‑148 (MQ=255)                                                 ccACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACAtcatca                                                                                                   <  1:225130/151‑1 (MQ=255)                                                   aCGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTTGCCCGCCCGGAT‑CACGGTCACGCCGCGAATCACCCCCTCATCaa                                                                                                  >  2:153746/1‑150 (MQ=255)                                                    cGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCaaa                                                                                                 <  1:1354351/150‑1 (MQ=255)                                                     ggTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGGCCCGGGT‑AATTT‑GCCCGCCCGGAT‑CACGGTGCAGCCCCGAAAAACCACATCATCaaaa                                                                                                >  2:1238403/1‑150 (MQ=255)                                                     ggTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAAt                                                                                               <  1:393056/151‑1 (MQ=255)                                                         aGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCaa                                                                                            >  1:712414/1‑150 (MQ=255)                                                          gCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTTGCCCGCACGGAT‑CACGGTTACGCCGCGAATCACCACATCATCAAAATCaa                                                                                            >  1:738742/1‑149 (MQ=255)                                                           cAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAg                                                                                         <  2:415323/151‑1 (MQ=255)                                                              gTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGTGATTTTCGCCC‑CACGGATCCACGG‑CACGCCGCGCATTACCACACCATCAAAATCcaca                                                                                          >  2:1146912/1‑147 (MQ=255)                                                                aCATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AAATTCCCCCGCCCGGGG‑CACTGGCAACCGGCGAAAGAACCCATCATCcaaaacta                                                                                            >  2:54110/1‑139 (MQ=255)                                                                 cATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGatat                                                                                   <  1:958809/151‑1 (MQ=255)                                                                  aTTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCCGGG‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATAtt                                                                                  >  1:347485/1‑151 (MQ=255)                                                                      ttCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCg                                                                              >  1:970913/1‑151 (MQ=255)                                                                      ttCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCCG‑CCCGGT‑AAATTTCCCCCCCCGGAA‑TACGGTCACGCCGCGAATCACCCCATCATCAAAATCAACAGTGATATTGcc                                                                               >  1:883827/1‑150 (MQ=255)                                                                       tCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCg                                                                              <  1:851029/150‑1 (MQ=255)                                                                        ccTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTTGCCCGCACGGGA‑CACGGGCACGCCGCGAATCACCACATCCTCAAAATCAACAGTGATATTGcc                                                                               >  2:946481/1‑148 (MQ=255)                                                                         cTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGGCGG‑CCCGGG‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATTACCACATCATCAAAATCAACAGGGATATTGCCGt                                                                             >  1:581047/1‑149 (MQ=255)                                                                          tCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCCCGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAAACAACAGTGATATTGCCGTCtt                                                                          >  1:1248855/1‑151 (MQ=255)                                                                          tCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGGCACGCCGCGCATCACCACATCATCAAAATCAACAGTGATATTGCCGTCtt                                                                          >  1:1243065/1‑151 (MQ=255)                                                                           cAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCCGGG‑AATTTTGCCCGCCCGGAT‑CACGGGCACGCCGCGAATCACCACATCATCAAAAACAACAGTGGTATTGCCGTCttt                                                                         >  2:144323/1‑151 (MQ=255)                                                                             gTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTtct                                                                       <  2:1100134/151‑1 (MQ=255)                                                                             gTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTACCGG‑CCAGGG‑AATTTCTCCCCCACGGAT‑CACGGGCACGCCGCGAATCAACACATCATCAAAATCAACAGAGATATTTTCGTCctt                                                                         >  2:800119/1‑149 (MQ=255)                                                                              ttttCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCCCATCATCAAAATCAACAGTGATATTGCCGTCTTtct                                                                       >  1:1369889/1‑150 (MQ=255)                                                                                ttCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTctct                                                                     <  1:1157191/150‑1 (MQ=255)                                                                                    cTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCCGGG‑AATTTCGCCCCCCCGGGT‑CACGGTCACGCCCCGAATCACCACACCATCAAAATCAACAGGGATATTTCCGTCTTTCTCttt                                                                   >  2:995662/1‑148 (MQ=255)                                                                                    cTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCCCACCGAT‑CACCGTTACGCCGCGAATCACCACATCCTCAAAAACAACAGTGATATTTCCGTCTTTCTtttt                                                                   >  1:649162/1‑148 (MQ=255)                                                                                     ttCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGcac                                                               >  1:250350/1‑151 (MQ=255)                                                                                      tCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTCGCCCACCCGGGG‑CACGGTCACGCCCCGAAACACCACATCATCAAAAACAAAAGGGATATTTTCCTCTTTTTCTTTg                                                                  >  2:716057/1‑147 (MQ=255)                                                                                        cGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCCGGT‑AATTTTGCCCGCACGGAA‑TACGGGCACGGCGAGAAACAACACATCATCAAAATCAACAGTGATATTTTCGTCTTTCTCTTTGCacaa                                                             >  2:576399/1‑150 (MQ=255)                                                                                         ggTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCacaaca                                                           <  2:132053/151‑1 (MQ=255)                                                                                          gTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAg                                                          <  2:1315099/151‑1 (MQ=255)                                                                                               gCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTca                                                     <  2:1020758/151‑1 (MQ=255)                                                                                                  ttttGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTcagc                                                   <  1:118778/150‑1 (MQ=255)                                                                                                    ttGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGAT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTcagca                                                  <  2:67601/149‑1 (MQ=255)                                                                                                     tgtgCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑GCACGG‑AATTTTTCCCCCACGGGT‑CACCGGCACGCCGCGAATCACCACAATATCAAAATCAACAATTATATTGCCGGCTTTCTCTTTGCATAAAAGCTTcagcag                                                 >  2:491977/1‑149 (MQ=255)                                                                                                     tgtgCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCCGGG‑AATTTCGCCCGCACGGAT‑CACGGTTACGCCGCGAAACACCACATCATCAAAAATAACAATGATATTTCCGTCTTTTTCTTTTCACAAAAGTTTTAGCAAAt                                               >  2:600631/1‑151 (MQ=255)                                                                                                       tgCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCCG‑CCAGGG‑AATTTCGCCCGCACGGAT‑CACGGGCACGCCGCGAATCACCACATCATCAAAAACAACAGGGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGAtt                                              >  1:844904/1‑150 (MQ=255)                                                                                                         ccgccTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCCGGG‑AATTTCCCCCGCACGGGT‑CACGGTCCCGCCGCGAAACACCACCACCTCCAAAACAACAGGGATATTGCCGTCTTTCTCTTTTCACAACAGGTTCAGCAGATTaa                                            >  2:834043/1‑150 (MQ=255)                                                                                                          cgccTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTTGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATTATCAAAATCAACAGTGATATTGCCGTCTTTTTCTTTGCACAACAGTTTCAGCCGAtt                                              >  1:1301701/1‑147 (MQ=255)                                                                                                           gccTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCCTCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTaa                                            >  1:1359326/1‑148 (MQ=255)                                                                                                            ccTGCGGCTGAGCTGATACCTGAATCGGCGGCGCCGG‑CCAGGT‑AATTTCGCCCGCCCGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAAACAACAGTGATATTTCCGGCTTTTTCTTTGCACAACAGTTTTAGCAGATTAACg                                          >  2:307886/1‑149 (MQ=255)                                                                                                             cTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTTCACAACAGTTTCAGCAGATTAACGAgg                                       >  1:834132/1‑151 (MQ=255)                                                                                                              tGCGGCTGAGCTGATACCTGAATCGGCGGCTCCCG‑GCAGGT‑AATTTTGCCCGCACGGGT‑CACGGGCACGCCGCGAATTACCACCTCATCCAAATACACAATGATAATGCCCTCTTTCCCTTTGCACAACAAGTTCAGCAGAATAAAGa                                         >  2:1188238/1‑148 (MQ=255)                                                                                                               gCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGt                                      <  1:996109/150‑1 (MQ=255)                                                                                                               gCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGG‑AATTTCGCCCGCACGGGT‑CACGGTCAACCCGCGCATCACCACCTCATCAAAATCAACAGTGAAATTGCCCCCTTTCTCTTTGCACAAAAATTTCAGCAGATTAAAGAGGtt                                     >  1:1140996/1‑151 (MQ=255)                                                                                                                 ggCTGAGCTGATACCTGAATCGGCGGTGGCGG‑GCCGGT‑AATTTTGCCCGGCCGGAT‑CACGGTCACGCCGCGCATCACCCCCTCATCAAAATCAACAGAGATATTTCCGTCTTTCTCTTTTCACAACAGGTTCAGCAGATTAAAGAGGttt                                    >  2:1308734/1‑150 (MQ=255)                                                                                                                  gctgagctgaTACCTGAAGCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGTTtgt                                  <  2:1215460/151‑1 (MQ=255)                                                                                                                    tgagctgaTACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTTTCTTTGCACAACAGTTTCAGCAGATTAACGAGGGTtgtg                                 >  1:477310/1‑150 (MQ=255)                                                                                                                          gaTACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGTTTGTGCCGTaaa                          <  1:1207900/151‑1 (MQ=255)                                                                                                                          gaTACCTGAACCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CTCGGTCACGCCGCGAATCGCCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTATGCACAGCAGTTTCAGCAGATTAACGAGGTTTGGGCCGTaaa                          <  2:753719/151‑1 (MQ=255)                                                                                                                                 gAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGTTTGTGCCGTAAAGCtgt                     <  1:667098/149‑1 (MQ=255)                                                                                                                                     cggcggTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGTTTGTGCCGTAAAGCTGTGAgg                 <  1:854184/149‑1 (MQ=255)                                                                                                                                        cggTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGTTTGTGCCGTAAAGCTGTGAGGATTGc            >  2:693254/1‑151 (MQ=255)                                                                                                                                        cggTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGTTTGTGCCGTAAAGCTGTGAGGATTGc            <  1:1207802/151‑1 (MQ=255)                                                                                                                                        cggTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGTTTGTGCCGTAAAGCTGTGAGGATTGc            <  2:1097342/151‑1 (MQ=255)                                                                                                                                           tGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGTTTGTGCCGTAAAGCTGTGAGGATTGCGTc         >  2:531488/1‑151 (MQ=255)                                                                                                                                                g‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGTTTGTGCCGTAAAGCTGTGAGGATTGCGTCGGCAg    <  2:161690/151‑1 (MQ=255)                                                                                                                                                   cAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGTTTGTGCCGTAAAGCTGTGAGGATTGCGTCGGCAGAc  >  1:523214/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                     CATCCAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAATCGGCGGTGCCGG‑CCAGGT‑AATTTCGCCCGCACGGAT‑CACGGTCACGCCGCGAATCACCACATCATCAAAATCAACAGTGATATTGCCGTCTTTCTCTTTGCACAACAGTTTCAGCAGATTAACGAGGTTTGTGCCGTAAAGCTGTGAGGATTGCGTCGGCAGAC  >  NC_000913/1676656‑1676951

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