Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,568,732	0	C	G	21.6%	85.0 / 10.4	37	P199P (CCG→CCC)	glgC	glucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (23/6);  new base G (0/8);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 7.78e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTTCCAGCGTACCCACATCGCGCCAGTACGGCTCGGCATCCGGGTCGGATTGTACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAACGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGG‑TTTTGGGTGACCGCATCT  >  NC_000913/3568484‑3568975                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       gCTTCCAGCGTACCCACATCGCGCCAGTACGGCTCGGCATCCGGGTCGGATTGTACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGAcggcgg                                                                                                                                                                                                                                                      <  1:417807/251‑1 (MQ=255)          gTACCCACATCGCGCCAGTACGGCTCGGCATCCGGGTCGGATTGTACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAgg                                                                                                                                                                                                                                             >  1:249119/1‑251 (MQ=255)                acaTCGCGCCAGTACGGCTCGGCATCCGTGTCGGATTGTACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCTGCTTCGGTTATCTTGGTAATCAAATCGTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACTTAGATACCCATAC‑TCGCCAGAGATTTGCTCGTATCGTTCGGCATTGACGGCGGGTTATAAGG‑TTTTTc                                                                                                                                                                                                                                      <  1:144425/251‑1 (MQ=255)                      cgcCAGTACGGCTCGGCATCCGGGTCGGATTGTACGCAAGAGAGCGGGCACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGGAGATCCCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTTACGGCGGGGTACCAGGGTTTTTCAACg                                                                                                                                                                                                                                  >  2:359405/1‑250 (MQ=255)                       gcCAGTGCGGCTGGGCATCCGGGTCGGATTGTACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACGCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTTGTATG‑GTTTTCAACGAA‑t                                                                                                                                                                                                                              <  2:172439/251‑1 (MQ=255)                           gTACGGCTCGGCATCCGGGTCGGATTGTACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAACGGA‑TTCgg                                                                                                                                                                                                                          >  2:216647/1‑250 (MQ=255)                           gTACGGCTCGGCATCCGGGTCGGATTGTACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAACGAA‑TTCGa                                                                                                                                                                                                                          <  2:138978/251‑1 (MQ=255)                                             gTCGGATTGTACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGGTAGCAGG‑GTTTTCAACGAA‑TTCGATAATTTTATCGTTTtcat                                                                                                                                                                                                        >  1:330524/1‑251 (MQ=255)                                               cGGATTGTACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTTAACCAA‑TTCGATAATTTTATCGTTCtcatc                                                                                                                                                                                                       >  2:70925/1‑250 (MQ=255)                                               cGGAGTGTACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAACGAA‑TTCGAGAATTTTATCGTTCtcatc                                                                                                                                                                                                       <  2:125417/250‑1 (MQ=255)                                                     gTACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAACGAA‑TTCGATAATTTTATCGTTCTCATCATCCGcc                                                                                                                                                                                                <  2:287396/251‑1 (MQ=255)                                                            agagagCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAACGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAAccc                                                                                                                                                                                         >  2:203506/1‑249 (MQ=255)                                                             gagagCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑GTTTTCAACGAA‑TTCGAATATTTTTTCGTTCTCATCAACCCCCCTAACGcc                                                                                                                                                                                        >  2:370404/1‑251 (MQ=255)                                                                  cGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAACGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATg                                                                                                                                                                                   >  1:99238/1‑251 (MQ=255)                                                                         gggTGAGCATAGCCCAGACCGACTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGCGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAACGACGTAGATACCCATAC‑TCGCCAGACATTTGCTCGGATCGTTCGGCATTGACGGCGGGGTCACCGG‑GTTTTCAACGAA‑ATCGCTAAATTTTTCGTTCTCATCAACCACCATAACGCCaaaagccgaagc                                                                                                                                                                            >  2:156831/1‑242 (MQ=255)                                                                          ggTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGGTAGCAGG‑GTTTTCAACGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAAAGCGGAGGct                                                                                                                                                                           >  1:80808/1‑251 (MQ=255)                                                                           gTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑GTTTTCAACGAA‑TTCGATAATTTTATCGTTCTCCTCAACCGCCATAACGCCCAATGCGGGGGctt                                                                                                                                                                          >  2:270892/1‑251 (MQ=255)                                                                                      ccAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAACGGA‑TTCGATAATTTTTTCGTTCTCATCAACCGCCATAACCCCAAATGCGGAGGCTTCTTCCATCggg                                                                                                                                                               >  2:205267/1‑250 (MQ=255)                                                                                             ggCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGGAGG‑TTTTTCAACGAA‑TTCGATAAATTTATTGTTCCCATCAACCGCCCTAACGCCAAATGCGGAGGCTTCTTCAATCGGTatt                                                                                                                                                            >  2:467005/1‑245 (MQ=255)                                                                                              gCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑GTTTTCAACGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCAt                                                                                                                                                       >  1:79249/1‑251 (MQ=255)                                                                                                    gTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAATCGGCGT‑CACAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGGTAGCAGG‑GTTTTCAACGCA‑TTTGATAATTTTATCGTTCTCATCAACCGCCCTAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGc                                                                                                                                                 >  1:385759/1‑251 (MQ=255)                                                                                                      gATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAAAGAA‑TTCGATAATTTTATCGTTCTTATCAACCGCCATAAAGCCAAATGCGGAGGCTTCTTAACTCGGTACTgggatcaaagca                                                                                                                                                >  1:244593/1‑240 (MQ=255)                                                                                                           tGGGAATAAAATCTTTGCAAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTTTTTTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑AAAAGACGTAGATACCCATAC‑CCGCAAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGTTTAGCAGG‑TTTTCCAACAAA‑TTCGATATTTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGt                                                                                                                                           <  1:270892/250‑1 (MQ=255)                                                                                                                         ttGCCAAAAGTGGTGGGTGGATTTCTCCCCGCGATCGTCTTCTTCC‑AGCATTTCAAACAAATAGCCGGCGT‑AAAAGACGTAGACCCCTATAC‑TCGCCAGAGTTTTGCTGGGTTCGTTCGTTTTGGGCGGGGGGTTAGCAGT‑TTTTACAACAAA‑TTCGATATTTTTATCGTTCTCACCAGCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCt                                                                                                                              <  2:271694/250‑1 (MQ=255)                                                                                                                              aaa‑GTCGTGGCTGGAGAGCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGCTAGTCGGCGT‑CAAAGACGGAGATACACATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTCAGCAGG‑CTTTTTAACCAA‑TTCGGTAATTTTATCGCTCTCAGCAACAGCCCTAACGACCAACGGGGAGGCGTCTTTCATTGGTACTGGCATACAAGCAACCGGGCAACATACGCCcttttttg                                                                                                                       >  2:125052/1‑248 (MQ=255)                                                                                                                                        ggtggtggtctccataggattctcttctTCCCAGCAGTTCTTCCAGAT‑GTCCGCGTGCAAATCCGTAGAT‑CCCCTACCTCGCCATAGTTTTGTTCGGATGGTTCGGCTTGGGCGGGTGGTAAGCAG‑‑TTTTTTAACCAAGTGCGTTATTTTTATCTTTCTCTCCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATct                                                                                                               <  2:328944/229‑2 (MQ=255)                                                                                                                                                       gcgATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAACGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCCATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGagg                                                                                               >  2:25631/1‑250 (MQ=255)                                                                                                                                                            cGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCCGAGATTTGCTCGGATCGTTCGCCATTGACGGCGGGTTAGCAGG‑GTTTTCCACGAA‑TTCGATAATTTTATCGTTCTCATCAACCCCCATAACCCCAAATGCGGAGGCTTCTTCAATCGGCACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGCCGCGCTCGGTAAGCCTACGCGAGTAGTct                                                                                          >  2:405177/1‑251 (MQ=255)                                                                                                                                                            cGTCTTCTCCC‑AGCAGTTAATACAGATAGTCGGCGT‑AAAAGACGTAAATCCCCATAC‑CCGCCAGAGTTTTGCTCGGATCGTTCGGAATTGACGGGGGGTTAGCAGT‑TTTTTAAACAAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTct                                                                                          <  1:216647/251‑1 (MQ=255)                                                                                                                                                                            gTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAACGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGGTATgg                                                                           >  1:33753/1‑251 (MQ=255)                                                                                                                                                                               cATACAGATAGTCGGCGT‑AAAAGACGTAGATACCCATAC‑CCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGGGGGTTAGCAGG‑TTTTTAAACGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTc                                                                         <  1:203506/250‑1 (MQ=255)                                                                                                                                                                                                                ccccattcccGCCAGAGTTTTGCTGGGTTCTTTCGGCATGGCCGGGGGTTTAGCAGT‑TTTTTAAACAAT‑TTCGATATTTTTTTCTTCCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTtat                                         <  1:290563/242‑1 (MQ=255)                                                                                                                                                                                                                          gCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGATTGGCAGG‑TTTTTCAACGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCAAAACGCCAAATGAGGAGGCTTCTTCAATCGGTACGGGCATACAAGCAACGGGGCAACGGGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGGAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGa                                >  1:145815/1‑250 (MQ=255)                                                                                                                                                                                                                             agagATTTGCTCGGATGTTTCGGATTTGAGGGGGGGTTAGCAGT‑TTTTTCAACAAA‑TTCGATATTTTTTTCGTTCCCACCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAAt                            <  1:408265/251‑1 (MQ=255)                                                                                                                                                                                                                                     gCTCGGATCGTTCGGCATTGCCGGGGGGTTAGCAGT‑TTTTTCAACGAA‑TTCAATATTTTTATCTTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGG‑t                   <  1:80825/251‑1 (MQ=255)                                                                                                                                                                                                                                        gggatgggtcggtatggACGGGGGGTTAGAATG‑TTTTACAACGAT‑TTCAATATTTTTTTCTTTCTAATCACCCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGG‑tttt                <  1:147994/235‑1 (MQ=255)                                                                                                                                                                                                                                              gTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAACGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCCAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTATTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAAcggcggggtattttgggggtttttgtgg           >  2:72500/1‑229 (MQ=255)                                                                                                                                                                                                                                                      ttGACGGCGGGTTAGCAGG‑GTTTTCAACGAA‑TTTGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTTTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGGAGTCTTGCTTGTAGATATGGTTGCCCGCCAGGATCACCCCGTATTCCGCTTTATAACGGCGGGTAAAGTCGAGGGTTTTGGG‑GACCGcctcc  >  2:331599/1‑247 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       GCTTCCAGCGTACCCACATCGCGCCAGTACGGCTCGGCATCCGGGTCGGATTGTACGCAAGAGAGCGGGAACGGGTGCGCATAGGCCAGACCGGCTTCGGTGATCTTGGGAATCAAATCTTTGCCAAA‑GTCGTGGCTGGAGTTCTCATCGCGATCGTCTTCTTCC‑AGCAGTTCATACAGATAGTCGGCGT‑CAAAGACGTAGATACCCATAC‑TCGCCAGAGATTTGCTCGGATCGTTCGGCATTGACGGCGGGTTAGCAGG‑TTTTTCAACGAA‑TTCGATAATTTTATCGTTCTCATCAACCGCCATAACGCCAAATGCGGAGGCTTCTTCAATCGGTACTGGCATACAAGCAACGGTGCAACGTGCGCCTTTTTCGACGTGATCGATAAGCATACGCGAGTAGTCTTGCTTGTAGATATGGTCGCCCGCCAGGATCACCACGTATTCCGCTTTATAACGGCGGATAATGTCGAGG‑TTTTGGGTGACCGCATCT  >  NC_000913/3568484‑3568975

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