Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,370,192	T→G	23.1%	H74P (CAC→CCC)	yjeH ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,370,192	0	T	G	23.1%	62.7 / 2.2	56	H74P (CAC→CCC)	yjeH	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (15/23);  new base G (6/10);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCACTCA  >  NC_000913/4369993‑4370411                                                                                                                                                                                                          |                                                                                                                                                                                                                             ccccggattACCGATATGCCCCCCCAGCGCCGGCGAACCGGATACCGCACACAACCGGTGCCAGCTGTGCCAGCCAAACAG‑AGCCTGGCCGAACCGGGCGGCAATTTGTAGTGCGGCAGGCAAACCCCCGGGAATGACCGATAAACACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTg‑g                                                                                                                                                                                                <  2:466010/221‑1 (MQ=255)   cgcgAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCCAACGCCATACCGACGAAGTGCGCGACGCCGCCTCCGCTGGGATAGCG‑Gc                                                                                                                                                                                               >  1:112169/1‑228 (MQ=255)                   ccACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCCAACCCACGGGAAAGACCGATAAAAAAA‑GCCAACCGGGGACCGGCCCAAGCCGCGGACCAAACGCCAAAACGACGAAAGGCGCGACCCCCCCCGCGCCGGGAAAGTG‑GGCACCCAGaaatccaaca                                                                                                                                                                              >  2:425931/1‑221 (MQ=255)                    caccCCAGGCCCAGCGTACCGAGTTCCGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGAAATTTGTAGTGCGGCAGGCAACCCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAac                                                                                                                                                                              <  2:9272/228‑1 (MQ=255)                      cacCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATGGCAAAcaca                                                                                                                                                                           >  1:246537/1‑229 (MQ=255)                         cagcgccagcgTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGCATGACCGATAAAAACA‑GCCAGCCGGTGCATCGCTCAAGCCGCGACCCAAACGCCAAACCCACGAAAGCCGCGCCGCCGCCCGCGCCGGGCACGGG‑GGGCCCCAGCATCGCAAACCcaaac                                                                                                                                                                        >  2:523966/1‑227 (MQ=255)                         cagcgccagcgTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGACTGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCCAACGCCATACCGACGCAGTGCGCGACGCCCCCTGCCCTGGGTTAATG‑GCGACCCCGAATCGCAAACAcaata                                                                                                                                                                        >  1:600650/1‑228 (MQ=255)                         cagcgccagcgTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCCGAATCGCAAACCcaaac                                                                                                                                                                        >  1:83362/1‑227 (MQ=255)                         cagcgccagcgTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAATG‑CCGACCCCGACTCCGAAACAcaatc                                                                                                                                                                        >  1:257537/1‑229 (MQ=255)                         cagcgccagcgTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAACCGCGAACCAAAAGCCATACCGCAGAAAGGCGCGACGCCCCCCGCGCGGGGAAAGGG‑GCACCCCAAAAACCCAAACaaaaac                                                                                                                                                                        >  2:325015/1‑227 (MQ=255)                            cgccagcgTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGCCCGATAAAAACA‑GCCAGCCGGTGACCCGCTCAAGCCGCGAACCAAACCCCAAACCGACGAAAGGCGCGACGCCCCCTGCCCCGGGATAGTG‑GCGCACCACAAACCCAAACAcaacccc                                                                                                                                                                      >  2:722554/1‑224 (MQ=255)                            cgccagcgTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACAAAACGCCATACCGACACAGAGCGCGACGCCGCCGCCCCTGGGATACTG‑GCAACCCAGAATCCCaacaccaaccgc                                                                                                                                                                      >  1:335696/1‑218 (MQ=255)                                  cgTACCGAGTTCTGCCAACAACAGTTGCAAGCTAAGCCAACCAAACAT‑CCCCTGGCCGAACCCGGCGGAAATTTGTAGTGAGGCAGGCAAACCCACGGGAATGCCCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCggg                                                                                                                                                               <  1:141216/229‑1 (MQ=255)                                  cgTACCGAGTTCGGCAAACAAAAGTTGCCAGCTATGCCCGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGCAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCggg                                                                                                                                                               <  1:346231/229‑1 (MQ=255)                                     aCCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGGGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAACAACAATCGCAATCGGGaaa                                                                                                                                                            >  1:171625/1‑228 (MQ=255)                                           ttCTGCACACACCGGTTGCCAGCTATGCCGGCCAAACAT‑CGCCGGGCCGAACCCGGGGGCAATTGGTAGGGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAg                                                                                                                                                      <  2:171625/229‑1 (MQ=255)                                                ccAACAACGTTTGCGAGTTAGGCCAG‑CACACATACGCCTGGCCGAGCCCGGCGGCAATGTGTGGTGCGGCAGGCAAACCCACGGGAAAGACCGAACAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAAt                                                                                                                                                 <  2:107272/229‑1 (MQ=255)                                                     aacaGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAACCCCACGGGAAGGACCGATAAAAACA‑GCCAGCCGGCGGCTCGCTCAAGCCGCGAACCAAACGCCAAACCGACGAAGTGCCGGACGCCCCCCGCCCTGGGGTTAGG‑GCGCCCCAAACTCGCCACCAAACACCCACTCGGGACCACCaaataacacaaaa                                                                                                                                            >  2:544330/1‑218 (MQ=255)                                                          ttGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGCCGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAAACGCAATCGGGAACCCTAAGATAATCAAAACGGGc                                                                                                                                       >  1:153065/1‑229 (MQ=255)                                                          ttGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCCAAACGGGc                                                                                                                                       >  1:11949/1‑229 (MQ=255)                                                            gCCAGTTATGCCAGCCAAACAC‑TCCCTGGCCAAACCCGGCGGCAATTTGTAGTGCGGCAGGCAACCCCACGGGAAGGACCGATAAAACCA‑GCCAGCCGGTGACCCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCa                                                                                                                                     <  1:482270/229‑1 (MQ=255)                                                               aGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGCG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGc                                                                                                                                  >  1:17110/1‑229 (MQ=255)                                                                  tATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCCCAAAAACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCca                                                                                                                               >  1:385488/1‑229 (MQ=255)                                                                  tATGCCAGCCAAACAC‑GGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGTAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCca                                                                                                                               <  1:425931/229‑1 (MQ=255)                                                                    ggccacccaAACAT‑CGCCTGGCCGAGCCCGGCGACAATTTGCAGCGCGGGAGGCAAACCCCCGGGAATGACCAAAAAAACCA‑GCCAGCAGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCTTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCcaca                                                                                                                             <  2:375128/228‑1 (MQ=255)                                                                     gccagccaAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGCTAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCACAATCGCAAACACAATCCCAATCGGGAACACCAAGATAAACAAAACCGGCCACACCcacc                                                                                                                             >  2:727364/1‑227 (MQ=255)                                                                      ccagccaAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGAGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAAACAAAACGGGCAACGCCCACAg                                                                                                                            >  1:637298/1‑228 (MQ=255)                                                                       aagccaAACCT‑CGCCGTGGCAAACCCGGCCGTAATTTTGAGTGGGGCAGACACACCCAGGGGAATGGCCGA‑AAAAACACACCGGCCGGGAGCTCGGCCAAGCCGCGAACCAAACGCCATACCGACGAAGGGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACTCAAAATCGCAAACACAAGCGGAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGc                                                                                                                          <  2:391140/228‑1 (MQ=255)                                                                                                  ggcggcAATTTGTAGTGGGGAAGGCAAACCCACGGGAACACCCGAAAAAAACC‑GCCCGCCGGGGACTCGCTCAAGCCGCGAACCAAACGCCATCCCGACGAAGTGCGCGACGCCGGCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCt                                                                                                <  2:199728/229‑1 (MQ=255)                                                                                                           tggcggggCGGCAGGCAAACCCCGGGAAAAGACCGATAAAACCG‑CCCGGCCGGTGACCCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGa                                                                                       <  2:11949/222‑1 (MQ=255)                                                                                                            ttggggccGGCAGGCAACCCCAGGGGAACGCCCGAAAAAACCA‑CCACGCCGTGGACTCGCACAAGCCGCGAACCAAACGCCATACCGGCGAAGTGCCCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGTCTGTTATTGCCCGCTACCAGAGCAGCTAACGCAGGaa                                                                                      <  2:385488/222‑1 (MQ=255)                                                                                                                   cGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCACCCACGCCAACCCCGGAACGGcca                                                                                >  1:188559/1‑226 (MQ=255)                                                                                                                    ggcaggcaAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGGGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCCCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCCGTTATTGCCCCCAACAACCGCAGCAACAGCAGGAACgcaaaaa                                                                              >  2:433842/1‑223 (MQ=255)                                                                                                                      caggcaAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGGACACTAAGATAACCAAAACGGGCCACGCCCCCAGGCGGTTATTGCCCGCAACCAGCCCAACTAACCCCGGAACGGCAAacac                                                                            >  2:217976/1‑229 (MQ=255)                                                                                                                         gcaAACCCACGGGAATGCCCGACAAAAACA‑GCCAGCCGGTGACTCGCGCAAGCCGCGAACCAACCGCCATGCCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAACCGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGc                                                                          <  1:433842/228‑1 (MQ=255)                                                                                                                          caAACCCACGGGAATGACCCAAAAAAACA‑GCCAGCCGGTGACCCGCTCAAGCCGCGAGCCAAACGCCAGACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCa                                                                        <  1:671066/229‑1 (MQ=255)                                                                                                                           aaaCCCACGGGAATGAACGATAAAAACG‑GCCAGCCGGTGACTGGCCCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCa                                                                        <  1:550344/228‑1 (MQ=255)                                                                                                                            aaCCCAGGGGAATGACCGATAAAAACA‑GCCAGCCGGTGCCTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCa                                                                        <  1:325015/227‑1 (MQ=255)                                                                                                                                     gacggACAGAAAAAACCA‑GCCAGCCGGGGCACCGCCCAAGCGGCGCACCAAACGCCAGCCCGACGAAGGCCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCCCGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCtaata                                                              <  2:415306/224‑1 (MQ=255)                                                                                                                                       atcaCCGAAAAAAACA‑GCCAGCCGGGGCCCCGCTCCAGCCGGTAACCAAACCCCAACCCGCAGAAGGGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAAAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATg                                                           <  2:246537/226‑1 (MQ=255)                                                                                                                                       agaaCCGAAAAAAACC‑GCCGGCCGGGGCCCCGCTCAACCCGCAAACCAACCGCACAACCGGCGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAACCCCAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATg                                                           <  2:185417/226‑1 (MQ=255)                                                                                                                                                  aaaaacccccgccGGGGTCCCGCCCAAGCCCAGAACCAAACGGCAGTCCGAGGAGGGGCGCGACGCCGCCTGGGCGGGGATAGGG‑GCGACCCAGAAGCGCAAACCAAATCGGAATCGGGAACACTAAAATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAg                                                 <  2:636048/219‑1 (MQ=255)                                                                                                                                                              cGGTGCCCCGCCAAAGCGGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAACCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGcc                                     <  1:544330/229‑1 (MQ=255)                                                                                                                                                               ggTGACCCGCCCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACCCAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGccc                                    <  1:722555/229‑1 (MQ=255)                                                                                                                                                                gTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAACGAAGCTAACCCAGGAACGGCAAACACGCCCGGGCCATAAAATGAACGCCATAACAGGCCaaatccc                                    >  2:614110/1‑223 (MQ=255)                                                                                                                                                                  gACTCCCGCAAGCCGGCAACCAAACGCCCATCCGACGAAGGGGCCGACGCCGCCTGCCCTGGGATAGTG‑GCGACCCAGAATCGGAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTgg                                 <  2:17110/229‑1 (MQ=255)                                                                                                                                                                     tccccaaggCCGGGAACCAAACGCCAAACCGGCGAAGGGCCGGACGCCGCCTGCGCGGGGATAGGG‑GCGACCCAGAATCGCAAACACAATCGAAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGcc                              <  2:188559/221‑1 (MQ=255)                                                                                                                                                                     ccgccccaGCCGCGAACCAAAAGCCATCACGACGAAGGGCACGACCCCGCCGGCGCTGGG‑TAGAGTGCGACCCAGAACCGCAACCCCAAACGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGcc                              <  2:600228/222‑1 (MQ=255)                                                                                                                                                                        cTCAAGCCGCGAACCAAACCCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCAACCCAGAACCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGcc                              <  1:697910/226‑1 (MQ=255)                                                                                                                                                                           aaGCCGCTAACAAAACGCCACACCGACGAAGTGCGCCACGCCGCCTGCGGTGGGATAGTG‑GCGACCCAGAAGCGCAAACACAATCGCAAGCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGATACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGccc                         <  2:112169/228‑1 (MQ=255)                                                                                                                                                                             gCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCCGGCCAAAACCCCGGGCCCACCCCAg                      >  2:126820/1‑229 (MQ=255)                                                                                                                                                                                    aCCAAACGCCATACCGACGAAGGGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGt               <  1:126820/229‑1 (MQ=255)                                                                                                                                                                                     ccAAACGCCACACCGACGAATGGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCCCGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGtt              <  1:341102/229‑1 (MQ=255)                                                                                                                                                                                           gCCATCCCGACGAAGTGCGCGACCCCGCCTGCGCGGGGATAGTG‑GCGCCCCAGAACCGCAAACACAACCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTc        <  1:280435/229‑1 (MQ=255)                                                                                                                                                                                              agcgcgacggAGGGCGCGACGCCGCCGGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCa      <  1:508447/224‑1 (MQ=255)                                                                                                                                                                                                 ccGACGAAGTGCCCGACGCCGCCTGCCCCGGGATAGGG‑GCAACCCAGAATCGCAAACACAATCGCAAGCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCACTCa  <  1:523966/229‑1 (MQ=255)                                                                                                                                                                                                          |                                                                                                                                                                                                                             ACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAT‑CGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGGCAGGCAAACCCACGGGAATGACCGATAAAAACA‑GCCAGCCGGTGACTCGCTCAAGCCGCGAACCAAACGCCATACCGACGAAGTGCGCGACGCCGCCTGCGCTGGGATAGTG‑GCGACCCAGAATCGCAAACACAATCGCAATCGGGAACACTAAGATAATCAAAACGGGCCACGCCCACAGGCTGTTATTGCCCGCTACCAGCGCAGCTAACGCAGGAACGGCAAACACGCCAGTGCCTAATAATGACGTCGATAGCAGGCCAATGCCCTGGGCCAGCCCCAGTTCTTGTTTGAGTCCACTCA  >  NC_000913/4369993‑4370411

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