Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,009,126	A→T	11.8%	Y428F (TAT→TTT)	rlmL →	23S rRNA m(2)G2445 and m(7)G2069 methyltransferases, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,009,126	0	A	T	11.8%	156.8 / 2.5	57	Y428F (TAT→TTT)	rlmL	23S rRNA m(2)G2445 and m(7)G2069 methyltransferases, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (35/12);  new base T (3/5);  total (40/17)
Fisher's exact test for biased strand distribution p-value = 9.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTAAGCTGCTTGCAGCTGCGTGCAGACAAACAGTACAAGGCGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACC  >  NC_000913/1008897‑1009355                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                         cTAAGCTGCTTGCAGCTGCGTGCAGACAAACAGTACAAGGCGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGAGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCCGAATATAACGTTGc                                                                                                                                                                                                                                 >  1:572032/1‑239 (MQ=255)         ctTGCAGCTGCGTGCAGACAAACAGTACAAGGCGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGCGAAGTTGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGGATGACGCCGATCTGCCCAGATAAAACCTTTCCCTTGaac                                                                                                                                                                                                                         >  2:27433/1‑237 (MQ=255)         ctTGCAGCTGCGTGCAGACAAACAGTACAAGGCGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGCTCTGCCAGAATATAACGTTGCCGTTGa                                                                                                                                                                                                                           >  1:163111/1‑237 (MQ=255)         ctTGCAGCTGCGTGCAGACAAACAGTACAAGGCGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGa                                                                                                                                                                                                                           >  1:815720/1‑237 (MQ=255)         ctTGCAGCTGCGTGCAGACAAACAGTACAAGGCGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGAcc                                                                                                                                                                                                                         >  1:23394/1‑239 (MQ=255)           tGCAGCTGCGTGCAGACAAACAGTACAAGGCGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGAAGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGt                                                                                                                                                                                                                       >  1:127434/1‑239 (MQ=255)                   cGTGCAGACAAACAGTACAAGGCGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATTTAACGTTGCCGTTGACCGTTATTCCg                                                                                                                                                                                                                >  1:356884/1‑238 (MQ=255)                                ggccaagggCAAAAACCGCCCCGGTGGCCTGCTAAAAAAAAAATCCCCAC‑TGGCCGAAAACCCCTCCAGCAAACAAGCCGGGGAGGGGGGGGAAAAACCCCACCACCCGCGTGCGTAAGAACCTAAAAAAATT‑CGAGAAGTGGGCTCGCCGGGAAGGGATTGAATGTTCCCGCCTGTATGACGCCGATCTGCCAGAATTTAACGTTGCCGTTGACCGTTATGCCGACTGGgtggtggtg                                                                                                                                                                                                  <  2:356884/232‑1 (MQ=255)                                      aGGCGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGggggggggggagggag                                                                                                                                                                                            >  2:402213/1‑231 (MQ=255)                                       ggCGAAAAACGGCCCGCTGGACTGCGTACAGCAAAATTACCACGTTGCCG‑AAAGCCCTCCAGACAGCAAACCGGCGATAGCAGCGGCAGACTACACCAACCCCCTGCGCAAGCACCTCACAAACGT‑CGAGTGGTGGGCCCGCCAGGCAGGGCTGGATTGTTACCGCCTGTATTCAGCCCACCGCCCAGAAAGTAACCTTGCACTTGACCTTTTGACCAACTGGGTGGGGGGGcgagaga                                                                                                                                                                                           >  2:897427/1‑233 (MQ=255)                                       ggCGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGGATGACGCCGATCTGCCAGAATATAACGTTGCCGTTTACCGTTATGCCGCCTGGGTGGTGGTGCAGGAGt                                                                                                                                                                                           >  1:239168/1‑239 (MQ=255)                                        gCGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACACCGACCTGCCGGAATTAAAAGTTGCCGTTGACCGTTATGCCGCCTGGGGGGGGGGGCAGGAgga                                                                                                                                                                                          >  2:188677/1‑237 (MQ=255)                                         cGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCAGCTATGGTAGCGGAAGCCTACACCAACCGCCTGCGTAAGAACATCCAAAAATT‑CGAGAGGTGGGCTCGCACGGACGGGTTTGAATGTTACCGCCTTCATGACGCCGACCTGCCAGAAACTCAAGTTGCCGTTGAACGTTCTGCCGCCTGGGGGGGGGGGCCGGAGtgg                                                                                                                                                                                         >  2:347089/1‑237 (MQ=255)                                            aaaaCGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGGGGTGCAGGGGAATgg                                                                                                                                                                                       >  1:383583/1‑237 (MQ=255)                                                                    aGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGAGGGCTCGCCAGGAAGGGCTTGACTGTTACCGCCTCCATGACTCCGATCTGCCAGAATTAAACGACGCCGTTCACCGTAATGCCGCCAGGTGGGGGGTGCCGGGGAATACG‑GCGCCCAAAAcatttcgagtccaa                                                                                                                                                             >  2:617211/1‑226 (MQ=255)                                                                                 cGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGCCTGGGTGGGGGGGCGGGAGTATTGG‑CCGCCAAACACTATTGATGCTCCCAAAGGGGGT‑CAGc                                                                                                                                               >  2:666132/1‑239 (MQ=255)                                                                                        g‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGACCACAAAGCGCGT‑CAGCGCCTGt                                                                                                                                         >  1:321571/1‑238 (MQ=255)                                                                                                 tCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATACG‑CCGCCAAAAACTATTGATGCTCACAAAGCCCGT‑CAGGGTCGTTTCGata                                                                                                                                   >  1:344947/1‑236 (MQ=255)                                                                                                        aGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCAAGACTATAACGTTGCCGTTGACCGTTATGCCGACTGGGGGGGGGT‑CCGGAGGTTGGGCCCGCCCAAAACTATTGATGTCAACAAAG‑GCGTCCAGCGtgctttcgaatataccccaga                                                                                                                         >  2:12798/1‑219 (MQ=255)                                                                                                          cAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTCTGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGGTCGGTATTATCGCTTCa                                                                                                                        >  1:278011/1‑238 (MQ=255)                                                                                                          cAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTa                                                                                                                                                                                          >  1:719503/1‑174 (MQ=255)                                                                                                           aaaCCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCACAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGGATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGCCGATGCGCCGGTGACGGAGTATGCC‑CCGCCACAAACAATTTATGCTAAAAAAGCGCGG‑CGGCCTCTGTTCGAtattatgcatcccac                                                                                                                      >  2:601511/1‑230 (MQ=255)                                                                                                             aCCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGGTGCCGTTGACCGTTATGCCGCCTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGG‑CAGCGTCTGTTCGCACTTATCGCTTCAACCa                                                                                                                    >  1:918798/1‑239 (MQ=255)                                                                                                                ggaggtgtggcgggAAGTCTACACACCCCCTCTGCGAAAGAACCCCAAAAATTT‑CGAAAAGGGGGCTCGCCAGGAAGGGTTTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCAttt                                                                                                                 <  1:450668/227‑1 (MQ=255)                                                                                                                gacGATGGTGGAGGAAGCCAACACAACCCGCCGGCATAAGAACCTCAAAAAATG‑CGAGAAGTGGGCTCGCCAGGAAGGTATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCAttt                                                                                                                 <  1:188677/237‑1 (MQ=255)                                                                                                                         aGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑gg                                                                                                       >  1:536920/1‑239 (MQ=255)                                                                                                                           cGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑gggg                                                                                                     >  1:63399/1‑239 (MQ=255)                                                                                                                               aGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGGGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CGCGGCCGGTTCGATATTATCGCTGCAACCATTTCGGTACG‑GGGGGtt                                                                                                  >  2:231670/1‑238 (MQ=255)                                                                                                                                cactccacaccctgCCTGTGTAAGACCCTAAAAAAATC‑CGAAAAGTGGCCCCGCCAGGAAGTGATTGATTGTAACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTgc                                                                                                <  1:27433/226‑1 (MQ=255)                                                                                                                                   tACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCggaa                                                                                                            >  1:709817/1‑223 (MQ=255)                                                                                                                                     cacCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCcaa                                                                                           >  1:534302/1‑239 (MQ=255)                                                                                                                                        cAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGCCTGGGTGGTGGTGCCGGAGTATGCG‑CCGCCAAAAACTATTTATGCCCCCAAAGCGCGT‑CGGCGCCTGTTCGATATTATCGCCGCAAACCTTTCggacacgggagttgcccacaaaa                                                                                         >  2:700662/1‑218 (MQ=255)                                                                                                                                          ccggcggggggaaGACCCTAAAAAAATG‑CGGGAGGGGGGCTGGCTAGGAAGGAATGGAGTGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCcaaacaaa                                                                                      <  1:142130/228‑1 (MQ=255)                                                                                                                                           ccGCCTGCGTAAAACCCTAAAAAATTG‑CGAGAAGTGGGCTCGCAAGGAAGGAATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCaaacaaac                                                                                     <  1:183097/239‑1 (MQ=255)                                                                                                                                                tGCGTAAACACC‑CAAAAAATTGAAAGGAGGGGGCTGCCCGGGAAGGGATTGGATGGTCCCGTCTTGATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTg                                                                                <  2:815720/239‑1 (MQ=255)                                                                                                                                                            tCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGCATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGGGCCAAACAAACTGGTGCTGAAAACCCGt                                                                    >  1:520629/1‑239 (MQ=255)                                                                                                                                                             cAAAAAATG‑CAGGAGGGGGGCTCGCAGGGAGGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTg                                                                   <  1:528717/239‑1 (MQ=255)                                                                                                                                                                      tgagggggggggCCCGCGAGGAGGGGATGTAATGTGCCCCCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAg                                                           <  2:127434/230‑1 (MQ=255)                                                                                                                                                                                    ggcagggatgggtttgttttttACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACGGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCa                                             <  2:23394/219‑1 (MQ=255)                                                                                                                                                                                    cgcagggAGGGGGTTAATTGTTACGGGCTGTATGCCGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCa                                             <  1:897427/234‑1 (MQ=255)                                                                                                                                                                                                          aCCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACTGGGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGaa                        >  1:282082/1‑239 (MQ=255)                                                                                                                                                                                                          aCCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAg                       >  1:475551/1‑239 (MQ=255)                                                                                                                                                                                                          aCCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGCCTGTTCGATATTATCGCTGCAACCATTTCGGCACG‑GGGGATTGCGCCAAACAAACTGGTGCGTAAAACCCGTGAACCCCAGAAGGGACAAAACCAAAACACAAAACAGGGGGGGAAg                       >  2:666981/1‑239 (MQ=255)                                                                                                                                                                                                              ccTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTgg                                                                                  >  2:306198/1‑176 (MQ=255)                                                                                                                                                                                                              ccTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTgg                                                                                  <  1:306198/176‑1 (MQ=255)                                                                                                                                                                                                               gggttttgccccGATCTGCCAGGAATTAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAATCTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGGCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGa                  <  2:738865/229‑1 (MQ=255)                                                                                                                                                                                                                          cgttCTTCCAAAATTTAACTTTGCCGTTGACCGTTATGCCGACTGGGTGGGGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAg       <  2:278011/236‑1 (MQ=255)                                                                                                                                                                                                                          cgcttttcccgAATTTAACTTTTCCCTTGACCGTTATGCCGACTGGGTGGTGGTGGAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAg       <  2:327462/229‑1 (MQ=255)                                                                                                                                                                                                                            aTCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCCAACAAACTGGTGCTGAAAACCCGTGAACCCCAGAAGGGCCAAAACCAATACCAGAACCTGGGCGAGAAGGGGCAATTTCTTGAAGta     >  2:311542/1‑238 (MQ=255)                                                                                                                                                                                                                            aTCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGACCGCCAGAAGGGCAAAAACAAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGtt     >  1:363499/1‑239 (MQ=255)                                                                                                                                                                                                                            aTCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCGGACGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAATTTCTTGAAgc      >  1:738865/1‑237 (MQ=255)                                                                                                                                                                                                                            aTCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCGGAAGGGCAAAAACCAAAACCAGAAACTGGGCGAGAAGGAGAGGTTTCTTAAAGta     >  2:155151/1‑238 (MQ=255)                                                                                                                                                                                                                            aTCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGtt     >  1:808104/1‑239 (MQ=255)                                                                                                                                                                                                                            aTCTGCCAGAATATAACGTTGCCGTCGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGCGAAGGGCGAGTTTCTTGACGtt     >  1:179307/1‑239 (MQ=255)                                                                                                                                                                                                                              cTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGGGAAGGGCGAGTTTCTTGAAGTTa    >  1:156126/1‑238 (MQ=255)                                                                                                                                                                                                                               tggtgaaaTTTAGCGGTGCCTGTGCCCGTTAGGCCGGCTGGGTGGGGGGTCCGGAATATGCG‑CCGCCAAAAACTATTGAGGCTCACAAAGCTCTG‑CAGCGTCTGGTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTAcc  <  2:890123/233‑1 (MQ=255)                                                                                                                                                                                                                                acaggAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTAcc  <  1:720146/234‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                         CTAAGCTGCTTGCAGCTGCGTGCAGACAAACAGTACAAGGCGAAAAACGGCCCGCTGGACTGCGTACAGAAAAATTACCACGTTGCCG‑AAAGCACTCCAGACAGCAAACCGGCGATGGTAGCGGAAGACTACACCAACCGCCTGCGTAAGAACCTCAAAAAATT‑CGAGAAGTGGGCTCGCCAGGAAGGGATTGAATGTTACCGCCTGTATGACGCCGATCTGCCAGAATATAACGTTGCCGTTGACCGTTATGCCGACTGGGTGGTGGTGCAGGAGTATGCG‑CCGCCAAAAACTATTGATGCTCACAAAGCGCGT‑CAGCGTCTGTTCGATATTATCGCTGCAACCATTTCGGTACT‑GGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACC  >  NC_000913/1008897‑1009355

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