Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,465,787	0	G	A	44.8%	74.8 / 39.0	67	R132C (CGC→TGC)	treB	trehalose‑specific PTS enzyme: IIB and IIC component
Reads supporting (aligned to +/- strand):  ref base G (31/6);  new base A (0/30);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.97e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCT  >  NC_000913/4465643‑4465927                                                                                                                                                     |                                                                                                                                               ccGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAAc                                                                                                                                           <  1:7624/150‑1 (MQ=255)    cTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑g                                                                                                                                        >  2:235846/1‑149 (MQ=255)       aCCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGAt                                                                                                                                   <  1:1137287/151‑1 (MQ=255)        ccAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATc                                                                                                                                  <  2:1278228/151‑1 (MQ=255)         cAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATc                                                                                                                                  >  1:956650/1‑150 (MQ=255)            cAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑c                                                                                                                             <  2:349944/151‑1 (MQ=255)             aaaTACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑cc                                                                                                                            >  1:76974/1‑151 (MQ=255)                    gaccgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCt                                                                                                                      >  1:1289816/1‑150 (MQ=255)                    gaccgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCt                                                                                                                      >  1:385313/1‑150 (MQ=255)                      cgggCAGGAAGAAGAATATCGCTTCCCCGACCACCCACAAAAAATCTTAGATCGTTTGCAGGGAAGGT‑T‑ACA‑TTTGCGCCGGCGTCTGCCCGTTCCTCGGGGGAAAA‑CCGCCGACCACTTTGCAAAACCC‑GAGAAAAAAA‑CCGCCGCTAAt                                                                                                                   <  2:26323/149‑1 (MQ=255)                        ggCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATc                                                                                                                  >  1:776615/1‑150 (MQ=255)                              tAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCgg                                                                                                           >  1:201626/1‑151 (MQ=255)                                 aagaagaTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAAcgcg                                                                                                            >  1:406228/1‑147 (MQ=255)                                 aagaagaTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGgc                                                                                                         >  2:915880/1‑150 (MQ=255)                                 aaaaagaTCGCTTCCCCGACAAACCAAAAAAAATCGTAGATGGTTTGCAGGAAAGGT‑T‑ACT‑TTTGCGCCAGCGTCTCCCCGTTCCTCAGGGAAAAA‑CCCCCGATCACATTGAAAAACCC‑GAGGATAAAA‑CCGCCGCTAATCAAcgcg                                                                                                            <  2:406221/144‑1 (MQ=255)                                 aaaaagaCCGCTTCCCCGATCACCCAAAAAAAATCGTAGATCGTTGGCGGGGAGGGG‑T‑ACT‑TTTGCGCCAGCGTCTGCCCGTTGCTCAGGGGAAAA‑CCGCCGATCACATGGCAAAAACC‑GGGGATAAAA‑CCGCCGCTAATCAACGCggg                                                                                                          <  1:383045/146‑1 (MQ=255)                                    aagaTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGACCACATTGCAAAACCC‑GAGGATAAAA‑CCGCCGCTAATCAACGCGGgcagc                                                                                                      <  1:1437209/150‑1 (MQ=255)                                        tCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACgg                                                                                                 >  2:137471/1‑151 (MQ=255)                                         cGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACggg                                                                                                >  1:1356280/1‑151 (MQ=255)                                             cccccGAACAACAACAAAAATTCGTAGATGGTTGGC‑GGGGAGGG‑T‑TCATTTTGCCCCGGCGTTTGACCGTTGCACAGGGGAAAA‑GCGCCGATCTAATGGCGAAACCA‑GAGGAACAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGgaaga                                                                                            <  2:1104032/148‑1 (MQ=255)                                              ccccGACCAACCAAAGAAAATCGTAGATCGTTTG‑AGGGAAGGGTT‑ACT‑TTTGCCCCAGCGTCTGCCCGTTGCTCATGGGAAAA‑TCGCCGACCACATGGCAAAACCG‑GAGGATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGgaagaa                                                                                           <  1:546099/149‑1 (MQ=255)                                                cgaaTCAACCAAAAAAAGTCGTAGAT‑GTTTGGCGGGAAGGG‑TAACT‑TGCGCCCCAGCTTCTTACGTTTCCCGGGGGGCAAA‑TCGCCAACCAAATGGAAAAACCG‑GAGTATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGgaaga                                                                                            <  1:1315930/145‑1 (MQ=255)                                                 cGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGgaaga                                                                                            >  1:670329/1‑147 (MQ=255)                                                      acccACAAAAAGTGGTAGTTCGTTGGCAGGGAGGGG‑T‑ACA‑TTCGCGCCAGCTTCTACCCGTTGCTCAGGGACAAA‑TCGCCGATCACATGGCAAAACCC‑GAGGATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCAc                                                                                     <  1:381468/147‑1 (MQ=255)                                                      aaCCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACcg                                                                                   >  1:663580/1‑151 (MQ=255)                                                         cacaGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACcgcg                                                                                 >  2:817864/1‑150 (MQ=255)                                                          acaGAAAGTCGTAGCTCTTTTGAGGGGAAGTG‑T‑ACT‑TTTCCCCCGGCGTCTGCCGGTGCCGAAGGGGAAAA‑TCGCGAACCACAGTGAAAAAACC‑GAGAAAAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGAt                                                                                        <  1:382857/142‑1 (MQ=255)                                                          acaGAAAATCGTAGATCGTTTGGGGGGAGGGT‑T‑ACT‑TTTGCGCCAGCCTCTGTCCGTTGCTGAGGGAAAAA‑TCGCTGACCAAATGGCAAAACCC‑AAGGATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGaa                                                                               <  1:84541/151‑1 (MQ=255)                                                            aaaaaaTCGTAGATC‑TTTGGAGGGAAGGGTT‑ACA‑TTTGCGCCAGCGTTTGACCTTTGCTCAGGGGAAAA‑CCGCCAACCACATTGCAAAACCC‑GAGAACAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGa                                                                                <  1:906701/146‑1 (MQ=255)                                                              aaaaTCGTATATTTTTTGGGGGGAGGGG‑T‑CTT‑TTTGCCCCAGCTTCTGCCCGTTGCTCATGGCAAAA‑CCGCAGACCCCATTGAAAAAACA‑GAGAATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAAt                                                                             <  2:1327265/149‑1 (MQ=255)                                                              aaaaTCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATg                                                                            >  2:555945/1‑150 (MQ=255)                                                              aaaaTCGTAGATCGTTGGGGGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGCCCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCAAAACCC‑GAGGATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATg                                                                            <  1:1070356/150‑1 (MQ=255)                                                                   cGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgag                                                                          >  2:606213/1‑147 (MQ=255)                                                                   cGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgag                                                                          >  1:133919/1‑147 (MQ=255)                                                                     tAGATCGTTTGCAGGGAGGGG‑A‑ACT‑TTTGCGCCGGCGTCTGCCCGTTGCTCAGGGACAAA‑GCCCCGACCAAATTGAAAAAACG‑GAGAAAAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGAt                                                                      <  2:1393779/149‑1 (MQ=255)                                                                      aaagtgTTGGCGGGGAGGGG‑T‑AAT‑TTGGCCCCAGCTTCTGCCGGTCGCTCA‑GGGCAAACCCGCCAAACCCATGGGAAAACCC‑GAGAATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagag                                                                        <  1:261253/141‑1 (MQ=255)                                                                          cGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGc                                                               >  1:178544/1‑151 (MQ=255)                                                                                cAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGc                                                         >  2:125550/1‑151 (MQ=255)                                                                                 aGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGcc                                                        >  2:1233607/1‑151 (MQ=255)                                                                                  ggaagggg‑T‑ACT‑TTTGCGCCCGCGTCTGCCCGTTGCTCAGGGCAAAA‑CCGCCGATCACATTGCAAAAACG‑GAGGATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGcc                                                        <  2:554023/145‑1 (MQ=255)                                                                                      agggtt‑ACT‑TTGGCGCCAGCGTCTGCCGTTTGCTAAGGGAAAAA‑CCGCCAATCAC‑TTGGAAAAACC‑CGGGAAAAAACCCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCAtt                                                     <  1:468412/145‑1 (MQ=255)                                                                                       ggg‑T‑ACT‑TTGGCGCCGGCGTCTGCCCTTTGCTCAGGGGAAAA‑TCCCCGACCACATTGCAAAACCC‑GAGAATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCAttt                                                    <  1:1061069/149‑1 (MQ=255)                                                                                        gg‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCat                                                 >  2:1236309/1‑151 (MQ=255)                                                                                               cttttcccccAGCGTCTCCCCTTTCCTCTTGGCCAAA‑CCCCCTATCACTTGGCAAACCCC‑GAGGACAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCAtatt                                              <  2:201623/142‑1 (MQ=255)                                                                                                 tttGCCCCAGCTTCTGCCGGTTCCTGGGGGAAAAA‑CCGCCGACCACATGGAAAAAACC‑GAGGATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCAtatt                                              <  2:1309772/148‑1 (MQ=255)                                                                                                     cgcCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCgg                                       >  2:49909/1‑151 (MQ=255)                                                                                                     cgacAGCGTCTGACCGTTGCTCCTGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCCAC‑CCGCCGCTAATCAACGCGGGCAGCAACGGGCAGAAGATCACCGCGAAATGAGAGATCAACTGCCCATGCCATTTCATATTATGCCgg                                       >  2:257851/4‑151 (MQ=255)                                                                                                     cccccgCGTCTGCCGGTTGCTCAGGGGAAAA‑TCGCCGATCACATGGAAAAAACC‑GAGGATAAAG‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCg                                        <  2:951301/145‑1 (MQ=255)                                                                                                       ccggCGTCGGCCCGTTGCTCGGGGGAAAA‑TCCCCGATCACTTGGAAAAAACG‑GAGAATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCg                                        <  2:477927/145‑1 (MQ=255)                                                                                                       ccAGCTTCGGACCGTTGCAAGGGGGAAAA‑CCGCCAATCACATGGAAAAAACG‑GAGAATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCgg                                       <  2:472349/149‑1 (MQ=255)                                                                                                           cGTCGGACCTTTGCTCAGGGGCAAC‑CCGCCGACCACATTGCAAAACCC‑GAGAATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGggcgg                                   <  2:685357/149‑1 (MQ=255)                                                                                                            gTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCtt                                >  1:233472/1‑151 (MQ=255)                                                                                                            gTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCtt                                >  2:663313/1‑151 (MQ=255)                                                                                                             tCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCttt                               >  1:1335918/1‑151 (MQ=255)                                                                                                                 cccGTTGCTAAGGGAAAAA‑TCGCCGACCAC‑TTGGAAAAACCGGAGAAAAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCttt                               <  2:765153/146‑1 (MQ=255)                                                                                                                  ccGTTGCTCATGGGCAAA‑CCGCCTATCCCATAGCAAAACCC‑AAGAACCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTa                          <  1:137474/151‑1 (MQ=255)                                                                                                                       gCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGc                     >  2:889688/1‑151 (MQ=255)                                                                                                                         tCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGctc                   >  1:633968/1‑151 (MQ=255)                                                                                                                             gggAAAA‑TCGCCGATCACATTGCAAAACCC‑GAGAAAAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGctct                  <  1:1353103/148‑1 (MQ=255)                                                                                                                             gggAAAA‑CCGCCGATCACATGGCAAAACCG‑GAGAATAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCttt                <  1:1302353/150‑1 (MQ=255)                                                                                                                             ggcaaaa‑CCGCCGCTCACATTGCAAAAACG‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCtt                 <  2:347428/145‑1 (MQ=255)                                                                                                                             ggacAAA‑CCCCCGACCACATTGCGAAACCC‑GGGGACCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCt                       <  2:713246/140‑1 (MQ=255)                                                                                                                             ggaaaaa‑CCGCCGACCACATGGCAAAAACC‑GAGAAAAAAC‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTg                      <  2:977464/140‑1 (MQ=255)                                                                                                                                  aaccccccGATCCCATTGCAAAAACC‑GAGAAAAAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCaa           <  2:956627/145‑1 (MQ=255)                                                                                                                                      cGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTg       >  2:265246/1‑151 (MQ=255)                                                                                                                                        ccGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTgcg     >  2:1196607/1‑151 (MQ=255)                                                                                                                                           aTCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCt  >  1:356485/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                               CCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGG‑T‑ACA‑TTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACC‑GAGGATCAAA‑CCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCT  >  NC_000913/4465643‑4465927

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