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

TTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGT  >  NC_000913/4465638‑4465934                                                                                                                                                        |                                                                                                                                                     tttCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCg                                                                                                                                                      >  2:974751/1‑150 (MQ=255)  ttCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGa                                                                                                                                                     >  2:634585/1‑150 (MQ=255)      ccGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAAcc                                                                                                                                                <  2:682286/151‑1 (MQ=255)       cGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCg                                                                                                                                               <  2:89675/151‑1 (MQ=255)           gACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑Agg                                                                                                                                           <  2:696308/150‑1 (MQ=255)           gACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑Agg                                                                                                                                           >  1:748346/1‑150 (MQ=255)                    aTACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAAccgc                                                                                                                                 >  2:41269/1‑151 (MQ=255)                         gaccgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCt                                                                                                                             >  1:106567/1‑150 (MQ=255)                         gaccgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCt                                                                                                                             >  1:20470/1‑150 (MQ=255)                             ggCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATc                                                                                                                         >  1:922438/1‑150 (MQ=255)                                aGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAAcg                                                                                                                     >  2:1159552/1‑151 (MQ=255)                                aGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAAcg                                                                                                                     >  1:175496/1‑151 (MQ=255)                                aGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAAcg                                                                                                                     >  1:102008/1‑151 (MQ=255)                                    agaagaagcTCCTTCCCCCGATCACCCACAAAAAATGGAATATCGTTGGAAGG‑GAGGGTTACTTTGCCG‑CCGGCGTCGGCCCGTTGCTCATGGGAAAACCCCCGATCTCATTGCAAAAACCG‑AGGAACAAACCGCCGCTAATCAAcgcg                                                                                                                   <  2:774342/149‑1 (MQ=255)                                      aagaagaTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGgc                                                                                                                >  2:336467/1‑150 (MQ=255)                                        gaagaTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGgcagc                                                                                                             >  1:781682/1‑151 (MQ=255)                                         aagaTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGgcagc                                                                                                             >  1:696389/1‑150 (MQ=255)                                           gaTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAAc                                                                                                          >  2:459263/1‑151 (MQ=255)                                           gaTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTATGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAAc                                                                                                          >  2:490376/1‑151 (MQ=255)                                               gCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACggg                                                                                                       >  2:869747/1‑150 (MQ=255)                                                 tcccccGATCAACCAAAAAAAATCGTAGATGGTTGGCGGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGCCCGTTGCTCTTGGGAAAATCGCCGACCACATTGCAAAACCCG‑AGGAAAAAACCGCCGCTAATCAACGCGGGCAGCAAc                                                                                                          <  1:497026/141‑1 (MQ=255)                                                   ccccGATCAACCACAAAAAATCGTAGATCGTTTGCAGG‑GAGGGGTACATTTGCG‑CCAGCGTCTGCCCGTTGCTCAGGGGAAAATCGCCGATCACATTGCAAAACCCG‑AGAAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGgaa                                                                                                     <  1:322975/146‑1 (MQ=255)                                                   cACCGAACAACCAAAAAAAACCGTAGATGGTTTGCGGG‑GAGGGGTACATTGGCG‑CCAGCGTCTGCCCGTTGCTCAGGGAAAAACCGCCGATCAATTTGAAAAACCCG‑AGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGgaa                                                                                                     <  1:1125797/148‑1 (MQ=255)                                                     ccgctaacccACAGAAAATCGTAGAGCGTTGGCAGG‑GAAGGTTACATTTCCG‑CCGGCGTCTGCCCGTTCCTAAGGGACAAATCGCGGATCACATAGCAAAACCCG‑AGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGgaagaa                                                                                                  <  2:1149656/141‑1 (MQ=255)                                                     ccGATCACCCAAAAAAATTCGTAGAT‑TTTTGGAGGAAAAGGAAACATTTGCC‑CCAGCGTCGGACCGTTGCTCAGGGGAAAACCCCCGACCAGATGGCAAAAACCG‑AGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGgaag                                                                                                    <  2:214575/147‑1 (MQ=255)                                                     ccGACCAACCAAAGAAAATCGTAGATCGTTTGCAGG‑GAGGGGTACTTTCGCG‑CCAGCGTCGCCCCTTTGCTCTGGGGAAAATCGCCAATCACATTGCAAAAACAG‑AGAAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGgaa                                                                                                     <  2:687344/146‑1 (MQ=255)                                                      cGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGgaaga                                                                                                   >  1:57913/1‑147 (MQ=255)                                                         tCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGGCCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCAc                                                                                            >  1:160211/1‑151 (MQ=255)                                                          cAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCAcc                                                                                           >  1:697099/1‑151 (MQ=255)                                                           acccaaaaaaaaTGGTAGATGGTTTGCGGG‑AAGGGGTACTTTGCCG‑CCAGCGTCTGCCCGTTGCTCAGGGCAAAACCGCCGATCACATTGAGAAACCCG‑AGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATc                                                                                              <  2:574741/139‑1 (MQ=255)                                                             ccACAGAAAATCGTAGATCTTTGGCAGG‑AAAGGGTACATTCGCG‑CCGGCGTCTGACCTTTCCCCAGGGACAAATCGCCGACCACATTGCAAAACCCG‑AGAATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACcgc                                                                                         <  2:210322/150‑1 (MQ=255)                                                             ccACAAAAAATCGTAGATCTTTTGAAGG‑GAGGGGTACATTTCCC‑CCAGCGTCTGCCGGTTCCTCAGGGCAAAACCGCCGACCACATTGCAAAACCCG‑AGAAACAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACcgcg                                                                                        <  1:994831/151‑1 (MQ=255)                                                                 aaaaaaTCGTAGATCGTTGGGGGG‑GAGGGGTACGTTCGCC‑CCACCGTCTGCCCGTTGCTAATGGAAAAATCGCCAATCAAATGGCAAAAGCCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGaaa                                                                                     <  2:1178466/148‑1 (MQ=255)                                                                 aaaaaaTCGAAGTTCGTTGGCAGG‑GAGGGGTCCTTTGGCG‑CCAGCGCCTGCCCGTTCCTCGTGGAAAAACCGCCAATCAAATTCCAAAAACGG‑AGGACAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACcgc                                                                                         <  2:967984/144‑1 (MQ=255)                                                                   aaaaTCGTAGTTCGTTGGCAGG‑GAGGGGTACATTCCCG‑CCGGCGTCTGCCCGTTGCTAAGGGGAAAACCCCCAACCACATTGAAAAAACGG‑AGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGaa                                                                                      <  2:323201/147‑1 (MQ=255)                                                                   aaaaTCGTAGATCTTTTGCAGG‑GAAGGGTACTTTGGCG‑CCGGCGTCTGCCCTTTGCTCAGGGGCAAACCCCCGATCACATTGCAAAACCCG‑AGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATga                                                                                  <  1:46148/151‑1 (MQ=255)                                                                    aaaTCGTAGATTGTTTGAGGG‑GAGGGGTACATGTGCC‑CCAGCGTCTCCCCTTTGCCCAGGGACAAATCCCCAAACACATTGAAAAAACGG‑AGGACAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgag                                                                                 <  2:466315/151‑1 (MQ=255)                                                                    aaaCCGTAGATCGTTTGCGGG‑GAGGGTTAAATTTGCG‑CCGGCGTCTGCCCTTTCCTAAGGGGAAAACCGCCAACCACATTGAAAAACCCG‑AGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgag                                                                                 <  1:224331/151‑1 (MQ=255)                                                                     aaTCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATg                                                                                   >  2:452622/1‑148 (MQ=255)                                                                      aTCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgag                                                                                 >  1:611471/1‑149 (MQ=255)                                                                          tAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATc                                                                            >  1:723143/1‑150 (MQ=255)                                                                          tAGATCGTTTGAGGG‑GAAGGTTCCTTTTGCC‑CCGGCGTCTGACGTTTGCCCATGGCAAAACCGCCAATCCCAGTG‑GAAAACCGAAGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagaga                                                                              <  1:1045009/148‑1 (MQ=255)                                                                           aGATCGTTGGCGGG‑GAAGGGTACATTGGCG‑CCAGCGTCTGCCCGTTGCTCAGGGGAAAACCGCCGACCACATTGAAAAAACGG‑AGGAACAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgagaga                                                                              <  2:307503/147‑1 (MQ=255)                                                                            gATCGTTGGCAGG‑GAGGGGTACATTGGCG‑CCAGCTTCTGACCTTTCCTCAGGGCAAAACCGCCGATCAAATGGAAAAAACAG‑AGAACCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCa                                                                           <  2:188399/149‑1 (MQ=255)                                                                             atggTTTGGGGG‑GAAGGTTACTTTGGCG‑CCAGCGTCTGACGGTTGCTCAGGGGAAATCCGCCGATCACATGGCAAAACCCG‑AGGAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACt                                                                        <  1:1002840/148‑1 (MQ=255)                                                                                  ttggaag‑gAAGGGTAACTTTGCGCCCGGCGTCTCCCCGTGGCCCAGGGGAAAACCCCCGACCACATGGAAAAAACCG‑AGGAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAAc                                                                         <  2:741195/140‑1 (MQ=255)                                                                                      aGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGcc                                                               >  2:519067/1‑151 (MQ=255)                                                                                       gg‑gAAGGTTACTTTTGCG‑CCAGCGTCTGCCCGTTGCTCATGGGAAAACCGCCGATCACATGGCAAAAACCG‑AGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGc                                                                <  1:369807/149‑1 (MQ=255)                                                                                            aGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTc                                                          >  1:167274/1‑151 (MQ=255)                                                                                             ggttACTTTGGCG‑CCAGCGTCGGCCCGTTGCTCAGGGGAAAATCGCCGATCACATTGCAAAAACCG‑AGAATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTTATGCCATTTCa                                                         <  1:227655/148‑1 (MQ=255)                                                                                             gggTACATTTGCG‑CCAGCTTCTGCCCGTTGCTCATGGGCAAACCCCCGATCACATTGCAAAACCCG‑AGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCAttt                                                           <  1:152179/149‑1 (MQ=255)                                                                                                     ttGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGc                                                 >  2:653860/1‑151 (MQ=255)                                                                                                      ttcg‑cCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATg                                                  >  2:995597/3‑149 (MQ=255)                                                                                                        cg‑cCAGCGTCTGACCGTTGCTCAGGGAAAAACCGCCGATCACATTGCAAAACCGG‑AGGAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCgg                                              <  1:1188615/151‑1 (MQ=255)                                                                                                         g‑cCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCggg                                             >  1:974161/1‑151 (MQ=255)                                                                                                           ccggcgttggACCGTCCCTCATGGACAAACCCCCTAACACTTGGCAAAACCCG‑AGTAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGggc                                            <  1:701196/142‑1 (MQ=255)                                                                                                            cAGCTTCTGCCCGTTGCTAATGGGCAAATCGCCAATCACATTGCAAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCgg                                              <  1:336472/148‑1 (MQ=255)                                                                                                                  cTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCtttt                                     >  1:943890/1‑151 (MQ=255)                                                                                                                     aCCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCttttttt                                  >  1:1002508/1‑151 (MQ=255)                                                                                                                     aCCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCttttttt                                  >  1:53130/1‑151 (MQ=255)                                                                                                                       cGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTAc                                >  1:494218/1‑151 (MQ=255)                                                                                                                       cGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTAc                                >  2:952054/1‑151 (MQ=255)                                                                                                                         tgcccctggggCAAATCGCCGATCACTTCGGAAAACCCG‑GGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCtt                                       <  2:550192/134‑1 (MQ=255)                                                                                                                                 gggAAAACCGCCGATCACATTGCAAAACCCG‑AGGAACAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTa                      <  1:41270/151‑1 (MQ=255)                                                                                                                                     aaaTCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCa                   >  2:284737/1‑150 (MQ=255)                                                                                                                                     aaaTCGCCAATCACATTGCAAAACCCG‑AGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTAt                     <  1:890932/148‑1 (MQ=255)                                                                                                                                     aaaCCGCCGATCACATGGAAAAACCGG‑AGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCaa                  <  2:792108/151‑1 (MQ=255)                                                                                                                                     aaaCCGCCGACCAATTTGAAAAACCCG‑AGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTg                             <  2:803135/140‑1 (MQ=255)                                                                                                                                     aaaCCGCCAATCACATGGAAAAAACCG‑AGGAACAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTAt                     <  2:20469/148‑1 (MQ=255)                                                                                                                                       aTCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATAAACTGCTCATGCCATTTCATATTATGCCGGACGGCTTTTTTTACCTGCTCTTTATCAAc                 >  1:90839/1‑150 (MQ=255)                                                                                                                                            cGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTgcgc           >  2:1000980/1‑151 (MQ=255)                                                                                                                                            cGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTgcgc           >  1:368771/1‑151 (MQ=255)                                                                                                                                                 acaTTGCAAAACCCG‑AGGAACAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTgc             <  1:812636/144‑1 (MQ=255)                                                                                                                                                     tGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGt  >  2:612618/1‑151 (MQ=255)                                                                                                                                                     tGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGt  >  2:363978/1‑151 (MQ=255)                                                                                                                                                     ggagAAACCCG‑AGGAACAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCgg   <  2:781663/147‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                     TTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGG‑GAAGGGTACATTTGCG‑CCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGT  >  NC_000913/4465638‑4465934

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: ‑