Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,400,207 G→A T294I (ACC→ATC)  mreB ← dynamic cytoskeletal protein MreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,400,2070GA97.3% 210.2 / NA 75T294I (ACC→ATC) mreBdynamic cytoskeletal protein MreB
Reads supporting (aligned to +/- strand):  ref base G (1/1);  major base A (41/30);  minor base C (2/0);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCT  >  NC_000913/3399965‑3400452
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
aaaaTTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGcccc                                                                                                                                                                                                                                               >  2:420208/1‑248 (MQ=255)
 aaaTTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCAcc                                                                                                                                                                                                                                               >  1:429177/1‑250 (MQ=255)
       gCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCcgc                                                                                                                                                                                                                                        >  2:461983/1‑251 (MQ=255)
          tCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAAAGGTCAAGGTTACGCAGCAGGGCGCCACCACCGATTAGCACCATGCcgcgtt                                                                                                                                                                                                                                     >  2:359238/1‑249 (MQ=255)
                          ttCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGc                                                                                                                                                                                                                     >  1:154090/1‑251 (MQ=255)
                                 atcagaccatgttggctAAACTGACACTTCCCCTGCCTGCATCCTATCTCTCTTCGCTGAACATGTCTCCGCCGGGCATGTCGATCATTTCCAGCGCTATGCCGCCACTGCGCGCCACTCAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATTAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTc                                                                                                                                                                                                               <  1:187344/234‑1 (MQ=255)
                                 aTCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTc                                                                                                                                                                                                               <  1:161123/250‑1 (MQ=255)
                                      aCCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCcggcgg                                                                                                                                                                                                         >  1:156821/1‑251 (MQ=255)
                                         aggcaggGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTTTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGcc                                                                                                                                                                                                      >  1:378973/1‑250 (MQ=255)
                                         aggcaggGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGGGCACCATGCCGCGCTCGGAGATGTCGGCAGCCAGTTCCGGCGGGc                                                                                                                                                                                                       >  2:48864/1‑250 (MQ=255)
                                          ggcaggGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCTATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCAc                                                                                                                                                                                                     <  1:397643/251‑1 (MQ=255)
                                            caggGTAAACAGACACTTCCCCTGCCTGAATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGCAGATGTCGGTAGCCAGTTCCGGCTGGCACtt                                                                                                                                                                                                   >  1:11070/1‑250 (MQ=255)
                                               ggTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTc                                                                                                                                                                                                <  1:420208/251‑1 (MQ=255)
                                                             tCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCAtt                                                                                                                                                                                   <  1:147565/250‑1 (MQ=255)
                                                              ccccTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTAc                                                                                                                                                                                 >  1:235837/1‑251 (MQ=255)
                                                                cctgcctgCATCCGATTACGCGTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACcg                                                                                                                                                                               <  1:350050/251‑1 (MQ=255)
                                                                   gcctgcATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCgcgc                                                                                                                                                                            >  1:89525/1‑251 (MQ=255)
                                                                        cATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTcaca                                                                                                                                                                       <  1:461983/251‑1 (MQ=255)
                                                                             gATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACAATTACAGCGCTcaaactacc                                                                                                                                                                  >  1:172563/1‑244 (MQ=255)
                                                                                 aCTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCCCTGTTTCCGTGCAACCATTACCGCCCTCACCATACAGGTc                                                                                                                                                              >  2:8836/1‑251 (MQ=255)
                                                                                     ttCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATGCCGGTCAGCg                                                                                                                                                          >  1:99188/1‑251 (MQ=255)
                                                                                          tGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGCCAGCGGTTcc                                                                                                                                                     >  1:342111/1‑251 (MQ=255)
                                                                                            aaCAGGGCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTCCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCTCTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTg                                                                                                                                                   <  2:342111/251‑1 (MQ=255)
                                                                                             aCAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGATCACAATACCGGGCAGCGGTTCCTGc                                                                                                                                                  >  1:379190/1‑251 (MQ=255)
                                                                                               agggcgccgcCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCCCCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAg                                                                                                                                                <  1:250431/247‑1 (MQ=255)
                                                                                                 gTCGCCGCCGTGCATGTCGCTCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCAATGTTCCAGTGCCACCATTCCCGCGCCCACAAAACCGGTCATCGGTTCCTGCAGTg                                                                                                                                              >  2:402892/1‑251 (MQ=255)
                                                                                                 gTCGCCGCCGTGCATGTCGAGCATTTCCAGCGCTATGCCGCCACCGCGCGCCACACAGGTCAGCTGGTCTGCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTg                                                                                                                                              <  2:52704/251‑1 (MQ=255)
                                                                                                    gctgccgTGCTTGTGGATCGTTTCCAGCGCTTTGCCGCCCCCGCGCGCCCCACAGGTCAGCGGGTCTCCAGCAACAACGACTGGAATGCCGGTTCCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATTCCGGTCAGCGGTTCCTGCAGTGCtt                                                                                                                                           <  2:28838/248‑1 (MQ=255)
                                                                                                           ggCATGTCGATCATTTCCAGCGCTTTGCCGCCGCCGCGCGCCACACAGGTCAGCGGGTCTTCAGCCACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTGCCGGCGTGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGAt                                                                                                                                    <  1:457813/250‑1 (MQ=255)
                                                                                                            gCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGgaga                                                                                                                                                                                                                               >  2:125622/1‑159 (MQ=255)
                                                                                                            gCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGgaga                                                                                                                                                                                                                               <  1:125622/159‑1 (MQ=255)
                                                                                                             cATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACCCAGGTCAGCGGGACTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGCTAAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGACAACGGTTCCTGCAGTGCTTCGAGGAtgt                                                                                                                                  >  1:178873/1‑249 (MQ=255)
                                                                                                             cATGGCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGAttt                                                                                                                                  <  1:399352/251‑1 (MQ=255)
                                                                                                                      cATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAg                                                                                                                          <  1:46646/250‑1 (MQ=255)
                                                                                                                            cAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCTggg                                                                                                                   >  1:136741/1‑251 (MQ=255)
                                                                                                                            cAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAggg                                                                                                                   <  1:201157/251‑1 (MQ=255)
                                                                                                                                     gccgccACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACcgc                                                                                                          >  1:166500/1‑251 (MQ=255)
                                                                                                                                        gccACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGCTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCCCTGTTACAGTGCAACCATTACCGCGCTCACAATACAGGTACGCGGTTCCTGCAAGGCTTCGAGGATTTACTTGGTGTTCAGGGTAAAAcacagtg                                                                                                       >  2:315227/1‑245 (MQ=255)
                                                                                                                                            ccGCGCGCCACCCAGGTCAGCGGGTCTTCAGCAACGACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGGTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGTGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAACCCGCGTGGAac                                                                                                   <  2:158739/251‑1 (MQ=255)
                                                                                                                                             cgcgcgCCACAAAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACAGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGAGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACTGGTCAGCGGGTCCTGCAGTGCTTTGAGGATTTCATTGGAGTTTCTGGTAAAACCGCGAGGaat                                                                                                   >  1:7873/1‑249 (MQ=255)
                                                                                                                                              gcgcgcCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAacac                                                                                                 <  1:363384/251‑1 (MQ=255)
                                                                                                                                                  gcCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGt                                                                                                                                                                                                                            <  2:424820/124‑1 (MQ=255)
                                                                                                                                                  gcCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGt                                                                                                                                                                                                                            >  1:424820/1‑124 (MQ=255)
                                                                                                                                                    cacacaGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCCTTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGGGGATTTCACTTGAGTTCAGGGTAAAACCCCGTGGAACACCTTCTg                                                                                           >  2:94340/1‑251 (MQ=255)
                                                                                                                                                         aGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTcc                                                                                                                                                     >  1:357549/1‑188 (MQ=255)
                                                                                                                                                         aGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTcc                                                                                                                                                     <  2:357549/188‑1 (MQ=255)
                                                                                                                                                              aGCGTGTCTTCAGCAAAAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGGCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACg                                                                                 <  1:226161/251‑1 (MQ=255)
                                                                                                                                                                cGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGc                                                                               >  1:187427/1‑251 (MQ=255)
                                                                                                                                                                  ggTCTTCAGCAACTACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATCCCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCa                                                                             <  1:109240/251‑1 (MQ=255)
                                                                                                                                                                  ggTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGCAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTAAGGGTAAAACCGCGTGGAACCCCTTCCGCCCAGTTACGGccc                                                                             >  2:209543/1‑250 (MQ=255)
                                                                                                                                                                            aacaacGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGa                                                                   <  1:4240/251‑1 (MQ=255)
                                                                                                                                                                             acaacGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCCCCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCCGTGCAACCATTACCGCGCTCACAGTACCCGTCAGCGGTTCCt                                                                                                                                                    >  2:65771/1‑169 (MQ=255)
                                                                                                                                                                             acaacGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGAGGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAGTACCGGTCAGCGGTTCCt                                                                                                                                                    <  1:65771/169‑1 (MQ=255)
                                                                                                                                                                               aaagacaggaaaggcGGTTTCTTCCAAAAAAAAACGGTCAAGGATACGCAACAGAGGGCCGCCACAGATGAGCACCATGCCGCACACGGAGATTACGGAAGCCAGTGCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCGTCGAGGATGTCATTGGAGTTCAGGGTAAAACCGTGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTAGAttt                                                                <  2:277747/237‑1 (MQ=255)
                                                                                                                                                                                  gACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCAc                                                             >  2:137110/1‑251 (MQ=255)
                                                                                                                                                                                            ccGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATc                                                   <  1:373689/251‑1 (MQ=255)
                                                                                                                                                                                                     tcCATTAACAAACGGTCAAGGTTACGCAGCAGGGCGCCCCCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAAACATTACCGCGCTCACAATGCCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGAt                                           <  2:99188/250‑1 (MQ=255)
                                                                                                                                                                                                     cccATTAACAAACGGTCAAGGTTCCGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATa                                          <  2:154090/250‑1 (MQ=255)
                                                                                                                                                                                                            aCAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCCACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCCTTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTATGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCCTCGCCCGGATAAACCGa                                    >  2:316017/1‑250 (MQ=255)
                                                                                                                                                                                                            aCAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGGGCCACCATTACCGCGCTCACAATAACGGTCAGCGGTTCCTGCAGTGCCTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGccca                                    >  1:422777/1‑248 (MQ=255)
                                                                                                                                                                                                                        gggtACGCAGCAGTGCGCCCCCCCCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCg                          <  2:15457/245‑1 (MQ=255)
                                                                                                                                                                                                                        ggTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCg                          >  1:15457/1‑248 (MQ=255)
                                                                                                                                                                                                                        ggTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGc                       >  1:327624/1‑251 (MQ=255)
                                                                                                                                                                                                                         gTTACGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGTGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCAGCGCCCGGATAAGCCGAACCGATTTCGTGCt                      <  2:224601/251‑1 (MQ=255)
                                                                                                                                                                                                                          ttACGCAGCAGAGCGCCCCCACCGATGAGCACCATGCCGCGCTCGGAGATGTCAGAAGCCAGTTCCGGCGGACGCTGTTCCAGCGCAACCATTACCGCGCTCCCAATACCGGTCAGCGGATCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCCGGGAAAAACCGCGTGGAACACCTTCTGCCCGGGTACGGCCAAGAACTTCTATTTCCCGGACTTCATCGCCCCTATAAGCCGAACCGATTCCGTGCtt                     >  2:211104/1‑251 (MQ=255)
                                                                                                                                                                                                                             cGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGAt                  >  1:237053/1‑251 (MQ=255)
                                                                                                                                                                                                                             cGCAGCAGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGAGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGCGGATTTCATTGGAGTTAAGGGTAAAACCGCGTGGGACACCTTCTGCCAGGTTACGGCCACGAACTCTGTTTTCCCGGACTTCCTCGCCCGGGTAAGCCGAACCGGTTTAGTGCTTGAt                  >  2:405670/1‑251 (MQ=255)
                                                                                                                                                                                                                              gcagcaGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGGGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTTGAGTTCAGGGTAAAACCCCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGACCCGATTTCATGCTTGAt                  >  2:347926/1‑250 (MQ=255)
                                                                                                                                                                                                                              gcagcaGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCt                                                                                            <  2:265490/176‑1 (MQ=255)
                                                                                                                                                                                                                              gcagcaGTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCt                                                                                            >  1:265490/1‑176 (MQ=255)
                                                                                                                                                                                                                               cagcagTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGGAAAACCGAGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCACCCGGATAAGCCGAACCGATTTCGTGCTTGATAc                >  1:278036/1‑251 (MQ=255)
                                                                                                                                                                                                                                   agTGCGCCACCACCGATGAGCACCATGCCGCGCTCGGCGATGTCGGCAGCCAGTTCCGGCGGGCACTGTTCCCGTGCAACCATTACCGCGCTCACAATACCGGTCAGAGGTTCCTGCAGTGCTTCGAGGATTTCATTGGGGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTAACGGACTTCATCGCCCGGATAAACCGAACCGATTTCGGGGTTAATACGtt             >  2:338105/1‑250 (MQ=255)
                                                                                                                                                                                                                                        gcCACCACCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGCTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGCTTTCACGGGCTTCATCGCCCGGATAAACCGAACCGATTTCGTGCTTGATACGTTCTGCgg       >  1:275642/1‑251 (MQ=255)
                                                                                                                                                                                                                                          caccacCGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTg     >  1:92075/1‑251 (MQ=255)
                                                                                                                                                                                                                                             caccGATGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAAAACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGAATTCACGGACTTCATCGCCCGGCTAAGCCGAGACGAATTCCTGCTTGATACGTTCTGAGGTGGcg  >  2:248877/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
AAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCT  >  NC_000913/3399965‑3400452

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