Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,358,859 T→C Y110C (TAT→TGT)  puuP ← putrescine:H(+) symporter PuuP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,358,8590TC96.8% 191.5 / NA 63Y110C (TAT→TGT) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base C (30/30);  minor base A (2/0);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

AGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATC  >  NC_000913/1358614‑1359101
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
aGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑Tcc                                                                                                                                                                                                                                                <  1:365150/251‑1 (MQ=255)
  aCGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCcag                                                                                                                                                                                                                                              <  1:503804/251‑1 (MQ=255)
   cGGTGCAAACGCATTCTCCTTTATGCAGTCCAGGAACTATCAGGAAGATAAACACCGCCATGATGGAGATTTGCACCAGTACAATCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATGCGCGGCGGTTAAAATGGCGACGAATGTTACCACCCACACCCACGGCCGCACTTCCGGGAAGAGGGCGGAGAGAGAGA‑TGTTCGCCAACAAGACGTGAATCATCGGCAAAAAGAGACAA‑TCcagc                                                                                                                                                                                                                                             <  2:32277/251‑1 (MQ=255)
     gTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCcagcag                                                                                                                                                                                                                                           <  1:274000/251‑1 (MQ=255)
                 tctcCTTGATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTGGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTGGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGc                                                                                                                                                                                                                                <  2:759207/250‑1 (MQ=255)
                             aGTCCCTGAACTACCAGGAAGATAAACTCCACCATGATGTAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACTCTCTGCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATcc                                                                                                                                                                                                                   <  2:262157/251‑1 (MQ=255)
                                cccTGAACTACCAGGAAGATAAACCCCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGAc                                                                                                                                                                                                                <  2:326012/251‑1 (MQ=255)
                                cccTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGAc                                                                                                                                                                                                                <  1:244609/251‑1 (MQ=255)
                                    gAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTGCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTg                                                                                                                                                                                                             <  1:116257/250‑1 (MQ=255)
                                     aaCTACCAGGAAGATAAACACCACCATGCTGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTTAAATGGCGACGAAGGGTACCACCCCAACCCACGACGGCACTTCCGGGAAGAGGGCGGGGAGATAGA‑TTTTCGCCAACAAGAAGATAATCCTCGGCAAAAAGAGAAAA‑TAACGCAGCGGTGACCAGCCCACCATaaaatcgacgggcg                                                                                                                                                                                                           >  1:661227/1‑240 (MQ=255)
                                         cccAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGtt                                                                                                                                                                                                       <  2:285500/250‑1 (MQ=255)
                                            cggAAGTTAAACACCACCACTAGGGAAATTTGCACCTGTACAAACAGGGTATTGAATTGAGCCACCTGGTTTACGCTCTTCAGTATCGAGGCGGTTAAAAGGGCGTCGAAGTTTACCACCCCCACTCACGGCGTCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTGAAGCACCGGCAAAAAGAGACAA‑TCCAGCTGCGGTGACCAGCCGACCCTACAACCGACGTGCTGGTTAAt                                                                                                                                                                                                    <  2:789076/250‑1 (MQ=255)
                                                aGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGAGCGGGGTAATCGAc                                                                                                                                                                                                >  1:605987/1‑251 (MQ=255)
                                                    aaaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACtttt                                                                                                                                                                                            >  1:289961/1‑251 (MQ=255)
                                                    aaaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACtttt                                                                                                                                                                                            <  1:647463/251‑1 (MQ=255)
                                                                gatgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGGCGAAGGTTACCACCCACACCCACGGCGGCACCTCCGGGGAGAGGGCGGAGAGGTAGA‑TTTTCGCCCACAAGACGTTTATCATCGGGCAAAAGAGACAA‑TCCAGCAGCGCTGACCAGCCGACCAAAAATCCGACGTGCGGGGTAATCGGCTTTTACGCGCAggca                                                                                                                                                                                 >  2:186015/1‑248 (MQ=255)
                                                                gatgGAGATTTGCACCAGTACAAACAGGGGATTGAAGTTAGCGCCCAGGTTGACGCTCGTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGta                                                                                                                                                                                 <  2:449721/250‑1 (MQ=255)
                                                                   ggAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAgg                                                                                                                                                                             <  1:655370/251‑1 (MQ=255)
                                                                           tcACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGGCGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTCCCACCCACACCCACGGCGGCGCTTCCGGGAAGAGGACGGAGATATAGA‑GTTTCGCCAACAAGACGTGAATCATCGGCAAAAAGAGACCA‑TCCAGCAGCGATGACCAGCCGATCATAAATTCGACGTGCAGGTTAATCGTCTTTTGCTCGTAGGTATAGGCCGAAcc                                                                                                                                                                      <  2:330587/249‑1 (MQ=255)
                                                                                  tACAACCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGATAAAATGGCGGCGAAGGGTACAACCCACACCCACGGCGGGACTTCCGGGAAGAGGGGGGAGAGATAGA‑TTTTTGACAACAAAACGTTAATCATCGGCAAAAAGAGACAA‑TACAGCAGCGGTGAACACCCGACCATAAAACCGACGAGCGGGGTAATAGACTATTGCGCGTAGGTATAggacggagcagccaccg                                                                                                                                                              >  1:29446/1‑236 (MQ=255)
                                                                                       aCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGGGCGGGGTAATCGACTTTTGCGCGGAGGGATAGGCCGAACCGGCCTCCGG‑AAAc                                                                                                                                                        >  1:184555/1‑251 (MQ=255)
                                                                                           ggTATTGAAGGTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGc                                                                                                                                                    <  1:330556/251‑1 (MQ=255)
                                                                                                 gAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGGGACGAAGGGTACCACCCCCCCCCACGGCGGCCCTTCCGGGGAGAGGGGGGGGAGATAAA‑TTTTTGCCAACAAGACGTTAATCATCCGCAAAAAGAGACAA‑ACCAGCAGCGGCGACCCGCCGCCCATAAAACCAAACTGCGGGGTAATCGACTTTTGCGCGTAGGTAAAGGCCGGACCCGCCTCCGG‑ACACAGGCGCAccc                                                                                                                                              >  2:380732/1‑250 (MQ=255)
                                                                                                  aaGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTTAAATGGGGGCGAAGGTTACCACCCCCACCCACGGGGGGACTTTCGGGAAGAGGGAGGGGAGATAGA‑TTTTTGCCAACAAGAAGGTAATCATCGGCAAAAAGAGAAAA‑ACCAGCAGCGATGACCAGCCGACCAAAAATCCGAAGTGCGGGTTAATAGACTTTTGTGCGTAGGTATAGGCCGAACCGGCCTCCGG‑GAACCGGCGAAcc                                                                                                                                               >  2:731425/1‑249 (MQ=255)
                                                                                                       aGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAACATGGCGGCGAAGGTTAACACCCACACCCCCGGCGGCACTTCCGGGGAGAGGGCGGAGAGATAGA‑TTTTTGCCCACAAGACGGTAATCATCGGCCAAAAGAGACAA‑ACCAGCAGCGGTGACAAGCCGACCATAAATCCGACGTGCGGGTTAAACGACTTTTGCGCGAAGGGATAGGCCGAACCGGCCTCCCG‑AAACTGGCGAACCAGGTTGc                                                                                                                                        >  1:697199/1‑251 (MQ=255)
                                                                                                           aCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCGACGAAGGGTACCACCCACCCCCACCGCGGCCCTTCCGGGCAGGGGGCGGAGAGATAGA‑TTTTCGCCAACCAAACGCTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTTCGCGTAGGGATAGGCCGAACCCGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGtt                                                                                                                                    >  1:415922/1‑250 (MQ=255)
                                                                                                             cAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCAAACCCACGGCGGCACTTCCGGGAAGAGGGGGGGGAGATAGA‑TTTTCGCCAAAAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCAAAAAACCGACGTGCGGGGTAATCGACTTTTTCGCGTAGGGATAGGCCGAACCGGCCTCCGG‑GAACTGGCGAAACAGTTTGCCGTAgg                                                                                                                                  >  1:1740/1‑250 (MQ=255)
                                                                                                               ggTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCCCACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCCACGTGCGGGTTAATCGACTTTTTCGCGTAGGTATAGGCCGAACCGGCCTCCCG‑AAACTGGCGAACCAGTTTGCCGTAGCTg                                                                                                                                >  1:332649/1‑251 (MQ=255)
                                                                                                                    aCGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCATCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGc                                                                                                                           <  1:118017/251‑1 (MQ=255)
                                                                                                                              gATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGA‑TTTTCGCCAACACGACGTTAAACATCGGCAAAAAAAGACAA‑ACCAGCAGCGATGAACAGCCGACCATAAATCCGACGGGCGGGGTAATCGACTTTTTCGCGTAGGTATAGGCCGAACCGGCCTCCGGGAAAC‑GGCGAACCAGTTTGCCGTAGCTGAAAGCGGTAAAAAGc                                                                                                                 >  1:698217/1‑251 (MQ=255)
                                                                                                                              gATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGc                                                                                                                 <  1:304662/251‑1 (MQ=255)
                                                                                                                               aTTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGca                                                                                                                <  1:441298/251‑1 (MQ=255)
                                                                                                                                    cggcggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACAccc                                                                                                           <  1:731425/251‑1 (MQ=255)
                                                                                                                                       cggTTAAAATAGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATCGGCCGAACCGGCATCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgcc                                                                                                        >  2:803673/1‑251 (MQ=255)
                                                                                                                                        ggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGATAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGGTTGCCGTAGCTTATAGCGGTAAACAGCACACCcgcc                                                                                                        >  2:752021/1‑250 (MQ=255)
                                                                                                                                         gatAAAATGGCGACGAAGGTTACCACACACACCCACGGCGGCACTTACGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCCTAAATCCGACGTGCGGGTTAATCGCCTTTTGCGCGTAGGGATAGGCCGAACCGGCCGCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAACGGTGAActgaaacccccccc                                                                                                       >  2:19772/3‑237 (MQ=255)
                                                                                                                                              aaTGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgccagcgccc                                                                                                 >  2:628021/1‑250 (MQ=255)
                                                                                                                                                        aaGGTTACCACCCACACCCACGGCGGCACTTACGGGAGGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAgg                                                                                        >  1:217998/1‑250 (MQ=255)
                                                                                                                                                          ggTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGcac                                                                                                               <  2:59151/225‑1 (MQ=255)
                                                                                                                                                          ggTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGcac                                                                                                               >  1:59151/1‑225 (MQ=255)
                                                                                                                                                           gTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGAAAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGc                                                                                    >  1:721777/1‑251 (MQ=255)
                                                                                                                                                            ttACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGcc                                                                                   >  1:793630/1‑251 (MQ=255)
                                                                                                                                                              cccacccacCCCCACGGCGGCACTCCCGGGAAGAGGGCGGAGAGATAGA‑TTTCCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGTTAAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCgg                                                                                 <  1:112989/250‑1 (MQ=255)
                                                                                                                                                               ccccccacccccccGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTACGCAAAAAAGACGTTAATCATCGGAAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTTATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGa                                                                                <  1:717468/248‑1 (MQ=255)
                                                                                                                                                                        cccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGCCGCCGGAACGTGACCg                                                                       >  2:168486/1‑251 (MQ=255)
                                                                                                                                                                          cacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTc                                                                     >  2:136583/1‑251 (MQ=255)
                                                                                                                                                                               cggcACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACAGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTa                                                                 >  1:638651/1‑250 (MQ=255)
                                                                                                                                                                                 gcACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAAAAAGACGTTAATCACCGGAAAAAAGAGACAA‑TCCAGCAGCGATGCCCAGCCGACAATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGTTATAGGCCGACCCGGCCCCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATg                                                              <  1:380732/251‑1 (MQ=255)
                                                                                                                                                                                    cgtCCGGGAAGAGGGCGGAGAGTAAAA‑TTGTCGCCACAAAGCCGTAAACAACCGGAAAAAAGAGACAA‑CCCAGAAGCGATGCCCAGCCGACCATAAATCCGACGTGGGGGTTAATCGCCTTTTGCGGGTAGGTAAAGGCCGAACCGGCCTCCGG‑AACCGGGCGAACCAGTTTGCCGTAGGTGAGAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCg                                                           <  2:332649/249‑1 (MQ=255)
                                                                                                                                                                                       ccGGGAAGAGGGCGGAAAAATAGA‑TTTCCGCCAACAAGACGTTAATCATCGGAAAAAAGAGACAA‑TCCAGCAGCGAGGACCAGCCGACAAAAAATCCGACGTGCGGTTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCCCCGG‑AAACGGGCGAACCAGTTTGCCGTAGCGGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGgag                                                        <  1:651657/251‑1 (MQ=255)
                                                                                                                                                                                         gggAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGAc                                                      >  2:242128/1‑251 (MQ=255)
                                                                                                                                                                                           gAAGAGGGGGGAGAAATAAA‑TTTCCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGCCCAGCCGCCCATAAACCCGACGGGCGGGTTAATCGACTTTTGCGCGTGTGTATAGGCCGAACCGGCCCCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCTCTCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACaa                                                    <  1:541039/251‑1 (MQ=255)
                                                                                                                                                                                            aaGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAt                                                   >  1:570403/1‑251 (MQ=255)
                                                                                                                                                                                            aaGAGGGCGGAGAGA‑AGAATTTCCGCCAACAAGACGGTATACAGCGGCAAAAA‑AGATAAACCCAGCAGCTATGACCAGCCGACGATAAACCCGGCGTGCGGGTTAGTCTACTTTTGCGCGTAGGTATGGGCGGAACCGGCCTCGTG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGAAAACAACGCACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAt                                                   <  1:791131/251‑1 (MQ=255)
                                                                                                                                                                                                gggCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAC‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATAcc                                                >  1:336023/1‑250 (MQ=255)
                                                                                                                                                                                                            aagtttttCGCCAAAAAGACTTTAATCAACGGAAAAAAAGGACAC‑TCCACCAGCGATGACCGTCCGACCAAAAATCCGCCGTGCGTGTTAATCGCTTTTTGCGCGTAGTTATAGGCCGACCCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCTTAGCGTATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCaaa                                    <  1:590334/247‑1 (MQ=255)
                                                                                                                                                                                                                    cGCCAACAAGACGTTAATCATCGGCAACAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTCGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGAATCAAATACGGTc                             >  1:625111/1‑250 (MQ=255)
                                                                                                                                                                                                                          caaGACGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCg                      >  1:691845/1‑251 (MQ=255)
                                                                                                                                                                                                                              aCGTTAATCATCGGCAAAAAGAGACAA‑TCCAGCAGAGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGAAAACAGCACTCCCGCCAGCTCCAGCAAATAGGACGCCGAAAAGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAAATCGGTCATCGGCGTg                    >  1:609685/1‑249 (MQ=255)
                                                                                                                                                                                                                                    atcatcGGCAAAAAGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCACATAGGACGCCGGACCGTGACCGTCGCTAACGCCGGAGACACTACCAAAAGTATCAAATACGGCCATCGGCGTGAGAAAGGc            >  2:97369/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aaaaGATACAA‑TCAAGCAGCGATGACCACCCGACCATAAATCCGCCGTGCGGTTTAACCGATTTTTGCGCTTAGTTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGACCCATTTTGCCGTAGCTAATAGCGGTAAACAGCACACCCGCCATCGCCAGCAAATAGGACGCCGGAACGTGACCGGCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCcatc  <  1:475197/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aaaaGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCacac                                                                                                             >  2:502594/1‑144 (MQ=255)
                                                                                                                                                                                                                                              aaaaGAGACAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCacac                                                                                                             <  1:502594/144‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
AGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGA‑TTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAA‑TCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGG‑AAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATC  >  NC_000913/1358614‑1359101

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