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,8590TC90.7% 148.9 / NA 55Y110C (TAT→TGT) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/26);  minor base A (5/0);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 5.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATC  >  NC_000913/1358613‑1359104
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         
cAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGCTGACGCTCTTCAGATTCGCGCCGGGTAAAATAGCGGCGACGGTTACCCCCCCCACCCACGGACCCACTTCCGGGAAGAGGGCGGAGGGATACATTTTTGCAAACAAGGCGTTCATCAACGCCAAAAAAAGAAAATc                                                                                                                                                                                                                                                       >  1:81505/1‑251 (MQ=255)
  gACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCCCACCCCCGGCGGCACTTCCGGGGAGCGGGCGGGGAGATAAATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGAAAAtacc                                                                                                                                                                                                                                                     >  2:139615/1‑248 (MQ=255)
    gggTGACAACTCCTTCTCCTGTATTCAGTCGCTGAACTACCAGGAAGATAAACAACACCATGAGGGAGAATTGCACCAGTTCAAACAGGGCATTGAAGTTAGCGACCAGGTTGACGCTCGTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTGCCACTCACGCCCGCGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGAGTGTCGCCACCAGGACGTTAATCCTCGGCAAAAAGAGACAATCcagc                                                                                                                                                                                                                                                  <  2:636537/250‑1 (MQ=255)
     ggTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCcagca                                                                                                                                                                                                                                                 <  1:91385/251‑1 (MQ=255)
                cTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCGACGAAGGATAACCCCCACACCCACAGAGGCACATCCCGGAAGAGGGCGGAGAGATAAAATTTCTCCAACAAGACGTTAATCATTGGCAAAAAAAGAAAATCCAACAACGATggac                                                                                                                                                                                                                                       >  1:588511/1‑246 (MQ=255)
                 ttCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGGAGAGGGCGGAGAGATAGATTTTCGCCCACAAAACGTTAATCATCGGCAAAAAGAGACAATCCCGCAGCGATGACCAGc                                                                                                                                                                                                                                    >  1:413167/1‑251 (MQ=255)
                            gcaattcatGAACCACCAGGAAGATAAACCCCACCATGATGGAGATTTGCCCCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTGAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGACAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCATCGATGACCTGCCGACCATaaa                                                                                                                                                                                                                         <  2:645876/242‑1 (MQ=255)
                                            cAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTaa                                                                                                                                                                                                        <  1:354583/251‑1 (MQ=255)
                                                aaGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCATGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGTGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACTTTTATCATCGGCAAAAAGAGACATTCCAGCAGCGATGGCCAGCCGGCTAGAAATCCGACGTGCGGGTTAATCGa                                                                                                                                                                                                    <  1:334611/251‑1 (MQ=255)
                                                 aGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGGAGAGGGCGGGGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGAAAATCCAGCAGCGATGACCAGCCGACCATAAATCCGAACTACGGGTTAATAGAc                                                                                                                                                                                                   >  1:76309/1‑251 (MQ=255)
                                                               atgatgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGAGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGt                                                                                                                                                                                     <  1:418549/251‑1 (MQ=255)
                                                                ggatgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGta                                                                                                                                                                                    <  1:212054/250‑1 (MQ=255)
                                                                   gggAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGGTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCTCCAACAAGACTTTAAGCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAAGCGACTTTTGCGCGTAGGTATAg                                                                                                                                                                                 <  2:252732/250‑1 (MQ=255)
                                                                    ggAGATTTGCACCCGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGGGACGAAGGTTACCACCCCCACCCACCGCGGCACTTCCGGGAAGAGGGGGGGGGGATAGATTTTCGCCAACAACACGTTAAACCTCGGCAAAAAAAGACAATCCAGCCGCGGTTACCACCCGACCCAAAAACCGACGAGCGGGTTAATTGCCttttttgcggagggatagg                                                                                                                                                                                >  2:41984/1‑236 (MQ=255)
                                                                    ggAGATTTGCACCAGTACAATCAGGGTATTGAATTTAGCGACCAGGTTGACGGTCTGCAGATTCTCGGCGGTTAAAATGGCGACGAAGGATACCCCCCACACCCACGGCGGTACTTCCGGGAAGAGGGCTGAGAGATATATTTTCGCCTACAAGACTTTAGTCATCGTCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAgg                                                                                                                                                                                <  1:806369/251‑1 (MQ=255)
                                                                    ggAGATTTGCACCAGTACAAACAGGGTATTGAAGTGAGCGACCAGGGTGACGCTATTCAGATTCGCGTCTGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGTCCCTTCCGGGAAGAGGGCGGAGAGAGAGATTTTCGCCAACAAGACGTTAATCACCGGCAAAAAGAGACAAGCCAGCAGCGATGACCAGCCGTCCATAAATCCGACGTGCGGGTTAATCGTCTTTTGCGCGTAGGTATAgg                                                                                                                                                                                <  2:589700/251‑1 (MQ=255)
                                                                      aggtttGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGGGGAGGAGGTAACCACCCCCACCCCCGGAGGCACATTCGGGGAGAGGGGGGGGGAGAAAATTTTCGCCAACAAAACGTTCATCAACGGCACAAAGAGACAAAACCGCAGCGATTGACCGAAGACCATAAAACCGGCGTGCTGGTAAATAGACATTTgcgagatggtattgggc                                                                                                                                                                              >  2:93163/4‑237 (MQ=255)
                                                                      agaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGcc                                                                                                                                                                              <  1:435334/251‑1 (MQ=255)
                                                                           tGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAAc                                                                                                                                                                          <  1:271809/250‑1 (MQ=255)
                                                                            gccccAGTACAAACAGGGTATTGAAGTTAGCGACCAGTTTGACGCTCTCCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTTCCACCCACTCCCACGGCGGCACTTCCGGGAAGAGGGCGTAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGTTGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGta                                                                                                                                                                                    <  2:752835/236‑1 (MQ=255)
                                                                            gCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCg                                                                                                                                                                        <  1:41984/251‑1 (MQ=255)
                                                                                               aTTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTTACAAGGGGACGAAGGGTACCACCCAAAACCACGGCGGCAATTACGGGGAGAGGGCGGAGAGAAAGAATTTCGCCAACAAGAAGTTAAACATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAAGCGACAATAAATACGACGTGCGGGATAAACGACTTTTGTGCGTAGGGATAGGACGAACCGGCCGCCgggaaatggaga                                                                                                                                                      >  1:199901/1‑240 (MQ=255)
                                                                                                          cGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGGAGAGGGCGGGGAGATAGATTTTCGCCAACAAAAAGTTAATCATCGGCAAAAAGAGACAATCCAGAAGCGGTGAACAACCGACCATAAATCCGCCGTGCGGGTTAATCGACTTTTGCGCGTAGGGGTAAGCCGAAACGGCCTCCGGCAACAGGCGAACCAGttttgcg                                                                                                                                          >  1:165861/1‑247 (MQ=255)
                                                                                                                 gTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGGAGAGGGCGGAGAGAGAGATTTTTGCCAACAAGACGTTAATCATCGGCAAAAAGAGAAAATCCAACAGCGATGACCCACCGACACAAAATCCGCAGCGGCGGTGAAACGACTTTTGCGCGCAGGTAAAAGCCGGAACGGCCTCCAGGAACTGTCGAAACAGGTTTCCGTATCTg                                                                                                                                   >  2:572982/1‑250 (MQ=255)
                                                                                                                      cGCTCTTCAGATGCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCg                                                                                                                             <  1:253901/251‑1 (MQ=255)
                                                                                                                           gtCAGACTCGGGGCGGTGAAAAGGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGGGGAGAGAGAGATTTTCGCCAACAAGGCGGTAAGCATCGGCAAAAAGAGACAATCCAGCATCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTGGTGCGCGTAGGGATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCATTTGGCCGTAGCTGATAGCGGTaaa                                                                                                                        <  2:860604/250‑1 (MQ=255)
                                                                                                                              aGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAg                                                                                                                     <  1:155966/251‑1 (MQ=255)
                                                                                                                                        cggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgcc                                                                                                           >  1:388471/1‑251 (MQ=255)
                                                                                                                                                    cgacgaAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGACAGCaa                                                                                                >  2:25503/1‑250 (MQ=255)
                                                                                                                                                                    ccacacccacGGCGGCACTTCCGGGAAGAGGGCGGCGAGATAGATTTTCGCCAAAAAGACGTTCATCATCGGCAAAAAGAGACAATCCCGCAGCGATGACCAGCCGACCCTACATCCGACGTGAGGGTTAATCGACTTTTGCGCGTAGGTATAGGACCAACCGGCCTCCGGAAACTGACACACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCCGCCCATAGGCCGCCGGAACGt                                                                               >  2:845962/1‑251 (MQ=255)
                                                                                                                                                                          ccacGGCGGCACTTCCGGGAGGAGGGCGGAGAGATAGATTTTCGCAAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGCCCAGCCGACCATAAATCCGACGTGCGGTTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGACCCAGTTTGCCGTAGCTGATAGCGGTAAACAGCCCACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGt                                                                         <  1:228928/251‑1 (MQ=255)
                                                                                                                                                                           cccgggggcccTTCCGGGAAGAGGGCGGAGAGATAGATTTCCGCCAACAAGACTTTAATCACCGGAAAAAAGAGACATTCCAGCAGCGATTACCAGCCGACCATAAATCCGACGTGCGGGTAAACCGCCTTTTGCGCGGAGGTAAAGGCCGAACCGGCCTCCGGAAACTGGCGAACCATTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTc                                                                        <  1:588350/241‑1 (MQ=255)
                                                                                                                                                                            ccggcggCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAACCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCg                                                                       <  1:139615/250‑1 (MQ=255)
                                                                                                                                                                              ggcggcACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCt                                                                     >  1:13332/1‑251 (MQ=255)
                                                                                                                                                                                   cccTTCCGGGAAGAGGGAGGAGAGATAGATTTTCGCCAACAAGACGTTAACCATCGGAAAAAAGAGACAACCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACGGGCGAACCAGTTGGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGc                                                                <  1:348143/249‑1 (MQ=255)
                                                                                                                                                                                            gAAGATGGCGGAGAGATAGATTTTAGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGTCCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGAGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGCTAGCGGTAAACCGCACACCCGCCAGCGCCAGAAAATAGGTCGCCGGAACGCGACCGTCCCTATTGCCGGGGACa                                                        >  2:621925/1‑250 (MQ=255)
                                                                                                                                                                                             tagagGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTGTAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGCTAGCGGTCAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAt                                                      >  2:74077/2‑251 (MQ=255)
                                                                                                                                                                                             aaGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCAAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAt                                                      >  2:239629/1‑251 (MQ=255)
                                                                                                                                                                                              agagGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTAGCTAATGCCGGAGActaa                                                      >  2:668055/1‑247 (MQ=255)
                                                                                                                                                                                               ggggggggAGATATAGATTTCGCACAACAAGAAGTTAAGCATCGGAAAAAAGAGACATTCCAGCAGAGAGGACCATCCGACAATAAATCCGACGTGCGGGTTAACCGACTTTTGCGCGTAGGTATAGGCCGATCCGGCCTCCGGAAACGGGCGAACCATTTTGCCGTAGCTGAAAGCGGAAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATAc                                                    <  1:559830/245‑1 (MQ=255)
                                                                                                                                                                                               gagGGCGGAAAAATAGATTTCCGCCAACAAGACTTTAAACATCGGCAAAAAGAGCCAATCCAGCAGCGATGACCAGCCGACCATAAACCCGACGTGCGGGTTAACCGACTTTTGCGCGTAGGTATAGGCCGAACCGCCCCCCGGAAACTGGCGAACCAGTTTGCCGTATCTGTTAGCGGTAACCAGCACTCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATAc                                                    <  1:411680/251‑1 (MQ=255)
                                                                                                                                                                                                    cGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCACATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGCGACAATACCaaa                                                >  1:196779/1‑250 (MQ=255)
                                                                                                                                                                                                     ggAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCaaa                                                >  1:559701/1‑249 (MQ=255)
                                                                                                                                                                                                           ggacagtgtcgcgcaaacaaGATGTTAACCAGGGGCAAAAAGGGAAATCCCAGCACCAATGCCCACCCGCCAATAAACCCGCCGGGCGGGTTAATCGACTCGTGGCCGTAGGTATAGGCCGACCCGGCCCCGGGAAACGGACAACCCAGTTGGCGGAGCTAGATAGGGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATc                                          <  1:552906/236‑1 (MQ=255)
                                                                                                                                                                                                                 ttttCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGCCGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTTGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGACAGCAAATAGGACGCCGGACCGTGACCGTCGCTACTGCCGGAGACACTACCAAAAGTATCAAATACg                                   >  2:254010/1‑251 (MQ=255)
                                                                                                                                                                                                                   ttCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGGGACCGTCGCTAATGCCGGCAACCATACCAAAAGGATCAAATACggg                                 >  2:117207/1‑250 (MQ=255)
                                                                                                                                                                                                                    tCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATCGCGGTAAACAGAACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGGGACCGTCGCTAATGCCGGGAACAATACAACAAGGATAAAATACGGTc                                >  2:166938/1‑251 (MQ=255)
                                                                                                                                                                                                                          aaaaaGACGTTAATCATCGGAAAAAAGATACAACCCAGCAGCGATGACCAGCCAACCAAAAATCCGAAGTGGGGGTTAACCGCCTTTGGCGCGAAGGAATAGGCCGACCCGGCCCCCGGAAACGGGCGAACCATTTTGCCTTAGTTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCgg                           <  1:657141/248‑1 (MQ=255)
                                                                                                                                                                                                                                   ttAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCGAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTgagct                   >  2:240345/1‑248 (MQ=255)
                                                                                                                                                                                                                                    aaaTCATGGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTGGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAg                 <  1:572982/250‑1 (MQ=255)
                                                                                                                                                                                                                                     aaTCATCGGCAAAAAGAGACAACCCAGCAGCGATGCCCAGCCAACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAgg                <  1:476603/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ggCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGAcc         >  2:199276/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGAccc        >  1:795471/1‑251 (MQ=255)
                                                                                                                                                                                                                                               aaaaaGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCcat      >  1:213183/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    gagaCAATCCAGCAGCGATGACCAGCCGACCATACATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGAATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCatcatc  >  1:230208/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         
CAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATC  >  NC_000913/1358613‑1359104

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: