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

AGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATC  >  NC_000913/1358614‑1359101
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
aGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATcc                                                                                                                                                                                                                                                <  1:365150/251‑1 (MQ=255)
  aCGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCcag                                                                                                                                                                                                                                              <  1:503804/251‑1 (MQ=255)
   cGGTGCAAACGCATTCTCCTTTATGCAGTCCAGGAACTATCAGGAAGATAAACACCGCCATGATGGAGATTTGCACCAGTACAATCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATGCGCGGCGGTTAAAATGGCGACGAATGTTACCACCCACACCCACGGCCGCACTTCCGGGAAGAGGGCGGAGAGAGAGATGTTCGCCAACAAGACGTGAATCATCGGCAAAAAGAGACAATCcagc                                                                                                                                                                                                                                             <  2:32277/251‑1 (MQ=255)
     gTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCcagcag                                                                                                                                                                                                                                           <  1:274000/251‑1 (MQ=255)
                 tctcCTTGATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTGGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTGGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGc                                                                                                                                                                                                                                <  2:759207/250‑1 (MQ=255)
                             aGTCCCTGAACTACCAGGAAGATAAACTCCACCATGATGTAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACTCTCTGCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATcc                                                                                                                                                                                                                   <  2:262157/251‑1 (MQ=255)
                                cccTGAACTACCAGGAAGATAAACCCCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGAc                                                                                                                                                                                                                <  2:326012/251‑1 (MQ=255)
                                cccTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGAc                                                                                                                                                                                                                <  1:244609/251‑1 (MQ=255)
                                    gAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTGCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTg                                                                                                                                                                                                             <  1:116257/250‑1 (MQ=255)
                                     aaCTACCAGGAAGATAAACACCACCATGCTGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTTAAATGGCGACGAAGGGTACCACCCCAACCCACGACGGCACTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACAAGAAGATAATCCTCGGCAAAAAGAGAAAATAACGCAGCGGTGACCAGCCCACCATaaaatcgacgggcg                                                                                                                                                                                                           >  1:661227/1‑240 (MQ=255)
                                         cccAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGtt                                                                                                                                                                                                       <  2:285500/250‑1 (MQ=255)
                                            cggAAGTTAAACACCACCACTAGGGAAATTTGCACCTGTACAAACAGGGTATTGAATTGAGCCACCTGGTTTACGCTCTTCAGTATCGAGGCGGTTAAAAGGGCGTCGAAGTTTACCACCCCCACTCACGGCGTCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTGAAGCACCGGCAAAAAGAGACAATCCAGCTGCGGTGACCAGCCGACCCTACAACCGACGTGCTGGTTAAt                                                                                                                                                                                                    <  2:789076/250‑1 (MQ=255)
                                                aGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGAGCGGGGTAATCGAc                                                                                                                                                                                                >  1:605987/1‑251 (MQ=255)
                                                    aaaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACtttt                                                                                                                                                                                            >  1:289961/1‑251 (MQ=255)
                                                    aaaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACtttt                                                                                                                                                                                            <  1:647463/251‑1 (MQ=255)
                                                                gatgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGGCGAAGGTTACCACCCACACCCACGGCGGCACCTCCGGGGAGAGGGCGGAGAGGTAGATTTTCGCCCACAAGACGTTTATCATCGGGCAAAAGAGACAATCCAGCAGCGCTGACCAGCCGACCAAAAATCCGACGTGCGGGGTAATCGGCTTTTACGCGCAggca                                                                                                                                                                                 >  2:186015/1‑248 (MQ=255)
                                                                gatgGAGATTTGCACCAGTACAAACAGGGGATTGAAGTTAGCGCCCAGGTTGACGCTCGTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGta                                                                                                                                                                                 <  2:449721/250‑1 (MQ=255)
                                                                   ggAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAgg                                                                                                                                                                             <  1:655370/251‑1 (MQ=255)
                                                                           tcACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGGCGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTCCCACCCACACCCACGGCGGCGCTTCCGGGAAGAGGACGGAGATATAGAGTTTCGCCAACAAGACGTGAATCATCGGCAAAAAGAGACCATCCAGCAGCGATGACCAGCCGATCATAAATTCGACGTGCAGGTTAATCGTCTTTTGCTCGTAGGTATAGGCCGAAcc                                                                                                                                                                      <  2:330587/249‑1 (MQ=255)
                                                                                  tACAACCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGATAAAATGGCGGCGAAGGGTACAACCCACACCCACGGCGGGACTTCCGGGAAGAGGGGGGAGAGATAGATTTTTGACAACAAAACGTTAATCATCGGCAAAAAGAGACAATACAGCAGCGGTGAACACCCGACCATAAAACCGACGAGCGGGGTAATAGACTATTGCGCGTAGGTATAggacggagcagccaccg                                                                                                                                                              >  1:29446/1‑236 (MQ=255)
                                                                                       aCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGGGCGGGGTAATCGACTTTTGCGCGGAGGGATAGGCCGAACCGGCCTCCGGAAAc                                                                                                                                                        >  1:184555/1‑251 (MQ=255)
                                                                                           ggTATTGAAGGTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGc                                                                                                                                                    <  1:330556/251‑1 (MQ=255)
                                                                                                 gAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGGGACGAAGGGTACCACCCCCCCCCACGGCGGCCCTTCCGGGGAGAGGGGGGGGAGATAAATTTTTGCCAACAAGACGTTAATCATCCGCAAAAAGAGACAAACCAGCAGCGGCGACCCGCCGCCCATAAAACCAAACTGCGGGGTAATCGACTTTTGCGCGTAGGTAAAGGCCGGACCCGCCTCCGGACACAGGCGCAccc                                                                                                                                              >  2:380732/1‑250 (MQ=255)
                                                                                                  aaGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTTAAATGGGGGCGAAGGTTACCACCCCCACCCACGGGGGGACTTTCGGGAAGAGGGAGGGGAGATAGATTTTTGCCAACAAGAAGGTAATCATCGGCAAAAAGAGAAAAACCAGCAGCGATGACCAGCCGACCAAAAATCCGAAGTGCGGGTTAATAGACTTTTGTGCGTAGGTATAGGCCGAACCGGCCTCCGGGAACCGGCGAAcc                                                                                                                                               >  2:731425/1‑249 (MQ=255)
                                                                                                       aGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAACATGGCGGCGAAGGTTAACACCCACACCCCCGGCGGCACTTCCGGGGAGAGGGCGGAGAGATAGATTTTTGCCCACAAGACGGTAATCATCGGCCAAAAGAGACAAACCAGCAGCGGTGACAAGCCGACCATAAATCCGACGTGCGGGTTAAACGACTTTTGCGCGAAGGGATAGGCCGAACCGGCCTCCCGAAACTGGCGAACCAGGTTGc                                                                                                                                        >  1:697199/1‑251 (MQ=255)
                                                                                                           aCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCGACGAAGGGTACCACCCACCCCCACCGCGGCCCTTCCGGGCAGGGGGCGGAGAGATAGATTTTCGCCAACCAAACGCTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTTCGCGTAGGGATAGGCCGAACCCGCCTCCGGAAACTGGCGAACCAGTTTGCCGtt                                                                                                                                    >  1:415922/1‑250 (MQ=255)
                                                                                                             cAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCAAACCCACGGCGGCACTTCCGGGAAGAGGGGGGGGAGATAGATTTTCGCCAAAAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCAAAAAACCGACGTGCGGGGTAATCGACTTTTTCGCGTAGGGATAGGCCGAACCGGCCTCCGGGAACTGGCGAAACAGTTTGCCGTAgg                                                                                                                                  >  1:1740/1‑250 (MQ=255)
                                                                                                               ggTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCCCACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCCACGTGCGGGTTAATCGACTTTTTCGCGTAGGTATAGGCCGAACCGGCCTCCCGAAACTGGCGAACCAGTTTGCCGTAGCTg                                                                                                                                >  1:332649/1‑251 (MQ=255)
                                                                                                                    aCGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCATCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGc                                                                                                                           <  1:118017/251‑1 (MQ=255)
                                                                                                                              gATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACACGACGTTAAACATCGGCAAAAAAAGACAAACCAGCAGCGATGAACAGCCGACCATAAATCCGACGGGCGGGGTAATCGACTTTTTCGCGTAGGTATAGGCCGAACCGGCCTCCGGGAAACGGCGAACCAGTTTGCCGTAGCTGAAAGCGGTAAAAAGc                                                                                                                 >  1:698217/1‑251 (MQ=255)
                                                                                                                              gATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGc                                                                                                                 <  1:304662/251‑1 (MQ=255)
                                                                                                                               aTTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGca                                                                                                                <  1:441298/251‑1 (MQ=255)
                                                                                                                                    cggcggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACAccc                                                                                                           <  1:731425/251‑1 (MQ=255)
                                                                                                                                       cggTTAAAATAGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATCGGCCGAACCGGCATCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgcc                                                                                                        >  2:803673/1‑251 (MQ=255)
                                                                                                                                        ggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGATAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGGTTGCCGTAGCTTATAGCGGTAAACAGCACACCcgcc                                                                                                        >  2:752021/1‑250 (MQ=255)
                                                                                                                                         gatAAAATGGCGACGAAGGTTACCACACACACCCACGGCGGCACTTACGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCCTAAATCCGACGTGCGGGTTAATCGCCTTTTGCGCGTAGGGATAGGCCGAACCGGCCGCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAACGGTGAActgaaacccccccc                                                                                                       >  2:19772/3‑237 (MQ=255)
                                                                                                                                              aaTGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgccagcgccc                                                                                                 >  2:628021/1‑250 (MQ=255)
                                                                                                                                                        aaGGTTACCACCCACACCCACGGCGGCACTTACGGGAGGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAgg                                                                                        >  1:217998/1‑250 (MQ=255)
                                                                                                                                                          ggTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGcac                                                                                                               <  2:59151/225‑1 (MQ=255)
                                                                                                                                                          ggTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGcac                                                                                                               >  1:59151/1‑225 (MQ=255)
                                                                                                                                                           gTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGAAAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGc                                                                                    >  1:721777/1‑251 (MQ=255)
                                                                                                                                                            ttACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGcc                                                                                   >  1:793630/1‑251 (MQ=255)
                                                                                                                                                              cccacccacCCCCACGGCGGCACTCCCGGGAAGAGGGCGGAGAGATAGATTTCCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGTTAAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCgg                                                                                 <  1:112989/250‑1 (MQ=255)
                                                                                                                                                               ccccccacccccccGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTACGCAAAAAAGACGTTAATCATCGGAAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTTATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGa                                                                                <  1:717468/248‑1 (MQ=255)
                                                                                                                                                                        cccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGCCGCCGGAACGTGACCg                                                                       >  2:168486/1‑251 (MQ=255)
                                                                                                                                                                          cacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTc                                                                     >  2:136583/1‑251 (MQ=255)
                                                                                                                                                                               cggcACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACAGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTa                                                                 >  1:638651/1‑250 (MQ=255)
                                                                                                                                                                                 gcACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAAAAAGACGTTAATCACCGGAAAAAAGAGACAATCCAGCAGCGATGCCCAGCCGACAATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGTTATAGGCCGACCCGGCCCCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATg                                                              <  1:380732/251‑1 (MQ=255)
                                                                                                                                                                                    cgtCCGGGAAGAGGGCGGAGAGTAAAATTGTCGCCACAAAGCCGTAAACAACCGGAAAAAAGAGACAACCCAGAAGCGATGCCCAGCCGACCATAAATCCGACGTGGGGGTTAATCGCCTTTTGCGGGTAGGTAAAGGCCGAACCGGCCTCCGGAACCGGGCGAACCAGTTTGCCGTAGGTGAGAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCg                                                           <  2:332649/249‑1 (MQ=255)
                                                                                                                                                                                       ccGGGAAGAGGGCGGAAAAATAGATTTCCGCCAACAAGACGTTAATCATCGGAAAAAAGAGACAATCCAGCAGCGAGGACCAGCCGACAAAAAATCCGACGTGCGGTTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCCCCGGAAACGGGCGAACCAGTTTGCCGTAGCGGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGgag                                                        <  1:651657/251‑1 (MQ=255)
                                                                                                                                                                                         gggAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGAc                                                      >  2:242128/1‑251 (MQ=255)
                                                                                                                                                                                           gAAGAGGGGGGAGAAATAAATTTCCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGCCCAGCCGCCCATAAACCCGACGGGCGGGTTAATCGACTTTTGCGCGTGTGTATAGGCCGAACCGGCCCCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCTCTCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACaa                                                    <  1:541039/251‑1 (MQ=255)
                                                                                                                                                                                            aaGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAt                                                   >  1:570403/1‑251 (MQ=255)
                                                                                                                                                                                            aaGAGGGCGGAGAGAAGAATTTCCGCCAACAAGACGGTATACAGCGGCAAAAAAGATAAACCCAGCAGCTATGACCAGCCGACGATAAACCCGGCGTGCGGGTTAGTCTACTTTTGCGCGTAGGTATGGGCGGAACCGGCCTCGTGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGAAAACAACGCACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAt                                                   <  1:791131/251‑1 (MQ=255)
                                                                                                                                                                                                gggCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACACTCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATAcc                                                >  1:336023/1‑250 (MQ=255)
                                                                                                                                                                                                            aagtttttCGCCAAAAAGACTTTAATCAACGGAAAAAAAGGACACTCCACCAGCGATGACCGTCCGACCAAAAATCCGCCGTGCGTGTTAATCGCTTTTTGCGCGTAGTTATAGGCCGACCCGGCCTCCGGAAACTGGCGAACCAGTTTGCCTTAGCGTATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCaaa                                    <  1:590334/247‑1 (MQ=255)
                                                                                                                                                                                                                    cGCCAACAAGACGTTAATCATCGGCAACAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTCGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGAATCAAATACGGTc                             >  1:625111/1‑250 (MQ=255)
                                                                                                                                                                                                                          caaGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCg                      >  1:691845/1‑251 (MQ=255)
                                                                                                                                                                                                                              aCGTTAATCATCGGCAAAAAGAGACAATCCAGCAGAGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGAAAACAGCACTCCCGCCAGCTCCAGCAAATAGGACGCCGAAAAGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAAATCGGTCATCGGCGTg                    >  1:609685/1‑249 (MQ=255)
                                                                                                                                                                                                                                    atcatcGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCACATAGGACGCCGGACCGTGACCGTCGCTAACGCCGGAGACACTACCAAAAGTATCAAATACGGCCATCGGCGTGAGAAAGGc            >  2:97369/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aaaaGATACAATCAAGCAGCGATGACCACCCGACCATAAATCCGCCGTGCGGTTTAACCGATTTTTGCGCTTAGTTATAGGCCGAACCGGCCTCCGGAAACTGGCGACCCATTTTGCCGTAGCTAATAGCGGTAAACAGCACACCCGCCATCGCCAGCAAATAGGACGCCGGAACGTGACCGGCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCcatc  <  1:475197/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aaaaGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCacac                                                                                                             >  2:502594/1‑144 (MQ=255)
                                                                                                                                                                                                                                              aaaaGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCacac                                                                                                             <  1:502594/144‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
AGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATC  >  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: