Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,162,216 G→A Q152* (CAG→TAG)  ftsP ← cell division protein FtsP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,162,2160GA96.8% 199.3 / ‑1.7 64Q152* (CAG→TAG) ftsPcell division protein FtsP
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (29/32);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

AGCAGTCGCAAGCGCACCCAGCCACGCGAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATC‑GTCGCCCTTCCAGACGCGGATAGTCGGCCCCAGGTAACGACCATTG  >  NC_000913/3161969‑3162460
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
agcagTCGCAAGCGCAACCAGCCACGCGAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACctactg                                                                                                                                                                                                                                                      >  2:186127/1‑251 (MQ=255)
    gTCGCAAGCGCACCCAGCCACGCGAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGcc                                                                                                                                                                                                                                                  >  1:252894/1‑251 (MQ=255)
         aggcgcACCCAGCCACGCGAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCg                                                                                                                                                                                                                                             >  1:442429/3‑251 (MQ=255)
         aaGCGCACCCAGCCACGCGAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCg                                                                                                                                                                                                                                             >  1:442437/1‑251 (MQ=255)
              cACCCAGCCACGCGAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTgg                                                                                                                                                                                                                                        <  1:208540/251‑1 (MQ=255)
                 ccAGCCACGCGAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTggg                                                                                                                                                                                                                                       >  1:105562/1‑249 (MQ=255)
                       ccgcgAGACTTCACCTTACGGGCTTTGTACCCCGTTAACCATCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTgg                                                                                                                                                                                                                               <  2:310143/250‑1 (MQ=255)
                         gcgAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCg                                                                                                                                                                                                                             <  2:426460/251‑1 (MQ=255)
                          cgAGACTTCAACGCAAGGGCTTTGTACACCGTTCACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCACGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAACATCATCCACACCATAATGGTTGGGGATCGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCACGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGt                                                                                                                                                                                                                            >  2:41642/1‑251 (MQ=255)
                               cgtCAAAGTACGGGCTTTGTACCCCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAAAGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATAc                                                                                                                                                                                                                       <  2:6240/249‑1 (MQ=255)
                                 acAACGAACGGGCGTTGTCCACCGTTAACCAGCAGCGTAGCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATAc                                                                                                                                                                                                                       <  2:442980/248‑1 (MQ=255)
                                                 gTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGtt                                                                                                                                                                                                      <  1:140543/250‑1 (MQ=255)
                                                        gTTAACCAGCAGCTTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGaa                                                                                                                                                                                              <  1:471010/251‑1 (MQ=255)
                                                           aaCCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATgg                                                                                                                                                                                           >  1:132254/1‑251 (MQ=255)
                                                           aaCCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATgg                                                                                                                                                                                           >  2:398919/1‑251 (MQ=255)
                                                            aCCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATggg                                                                                                                                                                                          <  2:252894/251‑1 (MQ=255)
                                                             ccAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGc                                                                                                                                                                                         >  1:89/1‑251 (MQ=255)
                                                               agctgccacgcACCAACAAAGCCTCCGCTTCCCGGTTCGTTGTATTCGGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAg                                                                                                                                                                                       <  2:423837/241‑1 (MQ=255)
                                                                                   gccgccGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGt                                                                                                                                                                  <  2:107187/251‑1 (MQ=255)
                                                                                   gccgccGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTGATCCAGCCGTTTAGCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGt                                                                                                                                                                  <  1:567842/251‑1 (MQ=255)
                                                                                      gccgcTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTg                                                                                                                                                               >  1:209320/1‑251 (MQ=255)
                                                                                      gccgcTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑GACGGGTGCCCAGTCAGCGTTTg                                                                                                                                                               >  2:328874/1‑251 (MQ=255)
                                                                                           ttgccGGTTCGTAGTATTCTGGCGTACCAACGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGATTGGGGATAGGCAGCGAATTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGCATGGGCAG‑TACGGGTGCCCAGGCAGCGTTTGG‑CGAc                                                                                                                                                         >  2:550289/4‑251 (MQ=255)
                                                                                             cccGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTGCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGAcat                                                                                                                                                       <  2:415871/251‑1 (MQ=255)
                                                                                                        gTATTCTGGCGTACCAAAGTTATCCAGCCGATTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGAAGCGACTTGCTGACTTCCTCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGGGTACCTGCGTTATGGGGAA‑GGGCAGGTATGGGGGCCCAGT‑AGGGGTTGGTGGACATCAACCGTGCCg                                                                                                                                            >  1:25747/1‑251 (MQ=255)
                                                                                                               ggCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGAcc                                                                                                                                        <  1:559354/248‑1 (MQ=255)
                                                                                                                    aCCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAACTCATCCACACCATAGTGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGGTCTGGCGAATGGGCAG‑TACGGGGGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATcc                                                                                                                                >  2:103982/1‑250 (MQ=255)
                                                                                                                        aaGTTATCCAGCCGTTTATCCTGGAGAATGACCGGAAAATCATCCCCACCATAATGGTTGGGTATAGGCTGCGACGTGCTGACTTCATCTTCTACCAGTCGCATTCCGGCAAGGCCGTGAGAGTCCTACTGGGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TTCGGGTGCTCATTTAGCGTTTTG‑CGACATCATCCGTGCCGGACCGCCCATCAGCgg                                                                                                                            <  2:371330/251‑1 (MQ=255)
                                                                                                                         agttattaaGCCGTTTATCCTGGATAATGATCGGATAATCATCCACCCCACAATGGTTGGGGATAGGCAGACCCCTGCTGACTTCCCCCTCCACCAGCCACATTCCGGCAAGGCCGTTAAAGACCTACTGCGCCGTGCGGTTGGGAGCATTGGCGTGCTACCACAGGTTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGCCAGCGTTTAG‑CGCCATCATCCGTGACGGACCGCCCACCACCggg                                                                                                                           >  2:223578/9‑251 (MQ=255)
                                                                                                                          gTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGCCTTGCTGACTTCATCTTCCAACAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACCTCATCCGTGCCGGACCGCCCATCAGCGGGc                                                                                                                          >  2:219309/1‑251 (MQ=255)
                                                                                                                           ttATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCCGCGGGcc                                                                                                                         >  1:432758/1‑251 (MQ=255)
                                                                                                                             aTCCAGCCGTTTATCCTGTAAAATGACCGTAAAATCATCCAAACCAGAATGGTTGGGGATAGGCAGCGATGTGCTGCCGTCACCTTCCTCCAGCCACATTCCGGCAAGGCCGTTATAGACCTACAGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCCGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCTGGCCt                                                                                                                        <  2:194192/250‑1 (MQ=255)
                                                                                                                                  gCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAAGGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGGCAGCGTTTGG‑CGCCATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTAcc                                                                                                                  <  2:328620/251‑1 (MQ=255)
                                                                                                                                       ttATCCTGGATAATGACCGGAAAATCATCCACGCCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTa                                                                                                              <  2:209156/250‑1 (MQ=255)
                                                                                                                                       ttATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTa                                                                                                              <  2:494194/250‑1 (MQ=255)
                                                                                                                                              ggagaaTGACCGGAAAATCATCCACACCATAAGGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGATCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGTCTGGTACCTGTAGCCCGGc                                                                                                       <  2:316402/246‑1 (MQ=255)
                                                                                                                                               gATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGAAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGcc                                                                                                      >  1:276256/1‑250 (MQ=255)
                                                                                                                                                     gACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACAGTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTc                                                                                                <  1:144208/250‑1 (MQ=255)
                                                                                                                                                      aCCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCAt                                                                                              <  1:248140/251‑1 (MQ=255)
                                                                                                                                                      aCCGGAAAATAATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGAACGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCAt                                                                                              >  2:382284/1‑251 (MQ=255)
                                                                                                                                                         ggAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTg                                                                                            >  1:208494/1‑250 (MQ=255)
                                                                                                                                                              atcatcCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACAt                                                                                      <  1:304135/251‑1 (MQ=255)
                                                                                                                                                                   ccACACCATACTGCCTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCAAATTCCGGCAAGGCCGTTATAGACCTCCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATCCCCCAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCCTTGAGCCATTTTcc                                                                                 >  2:263663/1‑250 (MQ=255)
                                                                                                                                                                               ggTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGACGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGCATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTg                                                                     >  1:82570/1‑251 (MQ=255)
                                                                                                                                                                                       aTAGGCGGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAAt                                                             <  2:105562/251‑1 (MQ=255)
                                                                                                                                                                                       aTAGGCAGCGACTTGCTGACTTCATCTTCCCCCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTGCTGTCAGGCGGTTGCTGTAAAt                                                             <  2:445199/251‑1 (MQ=255)
                                                                                                                                                                                                gACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAAc                                                    <  1:393407/251‑1 (MQ=255)
                                                                                                                                                                                                gACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGCTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGAGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGAACCTGTAGCCCGGCCACCGTCATTGAGACCTTTTCTGTCAGGCGGTTGCTGTAAATACGCTTAAc                                                    >  1:348675/1‑251 (MQ=255)
                                                                                                                                                                                                   ttGCTGACTTCATCTTCCACCAGCCACATTCCGGCACGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGCGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACAtc                                                 >  2:443197/1‑251 (MQ=255)
                                                                                                                                                                                                     gCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACAtc‑g                                               <  1:328874/250‑1 (MQ=255)
                                                                                                                                                                                                        gACTTCATCTTCCCCCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACAtc‑gtc                                             <  2:432758/249‑1 (MQ=255)
                                                                                                                                                                                                        gACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGCGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGCCATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATC‑GTCGc                                           >  2:365553/1‑251 (MQ=255)
                                                                                                                                                                                                        gACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACAtc‑gtc                                             <  1:39544/249‑1 (MQ=255)
                                                                                                                                                                                                        gACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATC‑GTCGc                                           >  1:258300/1‑251 (MQ=255)
                                                                                                                                                                                                         aCTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATC‑GTCGcc                                          <  1:398919/251‑1 (MQ=255)
                                                                                                                                                                                                           ttCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTACCATC‑GTCGCCCt                                        >  2:242990/1‑251 (MQ=255)
                                                                                                                                                                                                                       caGCTACATTCCGGCAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATC‑GTCGCCCTTCCAGACGCGGa                            <  2:320232/251‑1 (MQ=255)
                                                                                                                                                                                                                             caTTCCGGCAAGGCCGTTATAGACCGACGGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCGGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTGTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATC‑GTCGCCCGTCCAGACGCGGATAGTCg                      <  2:89/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGCAAGGCCGTTATAGACCTAATGAGCAGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACCTCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATC‑GTCGCCCTTCCAGACGCGGAGAGTCGGcccc                 >  2:360783/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cAAGGCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATATCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATC‑GTCGCCCTTCCAGACGCGGCTAGTCGGCCCCAgg              >  2:207446/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ggCCGTTATAGACCTACTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATC‑GTCGCCCTTCCAGACGCGGATAGTCGGCCCCCGGTa            >  2:294608/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ccGTTATAGACCTACTTATCCGTGCGGTATGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCTTTTTG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTTGTACCTTGAGCCCGGCCACCGTCAGTGAGACGTTTTCTGTCAGGCGTTTGCCGTGAATAAGCTTAACATCTGT‑GCCCTTCCAGACTCGGATAGTCGTCGCCATGTAACg         <  2:559023/251‑1 (MQ=255)
                                                                                                                                                                                                                                               ataGACCTGCTGCGACGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTGCCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGtt                                                                      <  2:407695/186‑1 (MQ=255)
                                                                                                                                                                                                                                               ataGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGtt                                                                      >  1:407695/1‑186 (MQ=255)
                                                                                                                                                                                                                                                 aGACCTACTAAGCCGTGCGGTTGGGAGTATTGGCGTGTTACCACAGAGTAGCTGAGTTATGGCGAATGGGCAG‑TGCGGGTGCCCAGTCAGCGTGTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGGAGCACGTCCACCGTCATTGAGCCATTTTCTGTCAGGCGGTTGCTGGAAATAAGCTTAACATC‑GTCGCCCTTCCATACGCGGATAGTCGGCCCCAGGCAACGACCATTg  >  2:271980/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
AGCAGTCGCAAGCGCACCCAGCCACGCGAGACTTCAACGTACGGGCTTTGTACACCGTTAACCAGCAGCGTATCACCAACAAAGCCGCCGCTTCCCGGTTCGTTGTATTCTGGCGTACCAAAGTTATCCAGCCGTTTATCCTGGATAATGACCGGAAAATCATCCACACCATAATGGTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAG‑TACGGGTGCCCAGTCAGCGTTTGG‑CGACATCATCCGTGCCGGACCGCCCATCAGCGGGCCTGGTACCTGTAGCCCGGCCACCGTCATTGAGACATTTTCTGTCAGGCGGTTGCTGTAAATAAGCTTAACATC‑GTCGCCCTTCCAGACGCGGATAGTCGGCCCCAGGTAACGACCATTG  >  NC_000913/3161969‑3162460

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