Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,869,637 C→A L87L (CTC→CTA cbrA → colicin M resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,869,6370CA90.0% 183.7 / 2.7 70L87L (CTC→CTAcbrAcolicin M resistance protein
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base A (39/24);  total (43/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGAT  >  NC_000913/3869396‑3869880
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
cgagtattGGCCGCGGCCCGTCAGGATCGGCGTTGGCACGAAAGGTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGGGGTACTGATGGTGTCAGCACACCCTGCGGCGGTCTGCTGGCACCTGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGAGCGCTAACCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAAc                                                                                                                                                                                                                                             <  2:420606/247‑1 (MQ=255)
cGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGTACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAAc                                                                                                                                                                                                                                             <  1:483945/251‑1 (MQ=255)
 gATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCCCCGGACGCGCCGCGGTCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGGGATCGCTAATCCGCAGATTTTTAGCGTCAAAACCGTCGGCGTCGCCGCCTCGCTCACGCGTAACt                                                                                                                                                                                                                                            >  2:5554/1‑251 (MQ=255)
     attagCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAAACCTGCGGCGGTCTGCTGGCACCGGACGAGCAACGGTCTTTTATTCGCGGTGGGCTGACGCTTCCTGTCGAGGGGATCGCTAAACCGCAGAGTTTCAGCGTCAAAACCGTCGACGACGCAGCATCGCGAACGCGTAACTAcc                                                                                                                                                                                                                                         >  2:517929/5‑250 (MQ=255)
       tGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGc                                                                                                                                                                                                                                      >  1:603838/1‑251 (MQ=255)
           cTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCCGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTTCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGGCGCCGCATAGCTAACGCGTAACTACCAGCGa                                                                                                                                                                                                                                    >  1:612418/1‑249 (MQ=255)
            ccGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGc                                                                                                                                                                                                                                 <  1:498263/250‑1 (MQ=255)
                    gCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCaa                                                                                                                                                                                                                         <  1:5554/251‑1 (MQ=255)
                      aGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAata                                                                                                                                                                                                                       <  1:150311/251‑1 (MQ=255)
                          ccGGCGTTGGCACCAAAGTTAGCAGGCAAAATGCAGGTTATCGCGCTGGATAAAAAGCACCAGTGTGGTACTTAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGATGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTaa                                                                                                                                                                                                                  <  2:682983/250‑1 (MQ=255)
                            ggCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATc                                                                                                                                                                                                                <  1:546000/251‑1 (MQ=255)
                              cGTTGGCACGCAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGAGACGCGCTCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGCTCGCTTCTCAGCAGCTTTTCAGCGTCAAAACCGTCGCCGACGCCGCCTAGATAAAGCGTAACTACACGCGAAGCTATATCAATATTAATcac                                                                                                                                                                                                              >  2:209320/1‑249 (MQ=255)
                               gTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGCAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGcc                                                                                                                                                                                                             >  1:680719/1‑251 (MQ=255)
                                 gggCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAATGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCa                                                                                                                                                                                                            <  2:224114/249‑1 (MQ=255)
                                   gCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGCTAAAAAGCACCAGTGTCCTACTGAAGGTTTCCGCAAACCCTGCGGAGGTCTGCTGGCCCCGGACGCGCCGCGCTCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTTAGCGTCCAAACCGCCGACGTCGCCGCATCGGTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATgg                                                                                                                                                                                                         >  1:336347/1‑250 (MQ=255)
                                       gAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTc                                                                                                                                                                                                     <  1:115044/251‑1 (MQ=255)
                                       gAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTAGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGACATGCtttt                                                                                                                                                                                                     >  2:527954/1‑250 (MQ=255)
                                        aacgTTAGCAGGCAAAATGCAGGCGATCGCGCGGGATAACAAGCACCAGTGAGGTATTGAAGGTTTCAGCAAACCCTGCGGCGATCTGCTGGCACCGGACTCGCAGCGATCTTTTATTCGCGATGGGCTCACGCGTCCTGTCGATGTGATCGCTCATCCGCAGATTTTCAGCGTCAAAACCGTCGACACCGCCGAATCGCTCACGCGACACTACCAGCAACGCTCTCTCACTATTAATTGCAATGCTTTCg                                                                                                                                                                                                    >  2:487313/4‑251 (MQ=255)
                                           gTTAGCCGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTAAGCAAACCCTGCGGCGGTCTGCTGGCACCGGAAGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAACACCGTCGACGTCCCCGCATCGCTAACGCGGAACTAACAGCGACGCTATATCAATATTAATCGCCATGCTTTCGACt                                                                                                                                                                                                 >  2:302389/1‑251 (MQ=255)
                                                    cAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATg                                                                                                                                                                                         >  1:451335/1‑250 (MQ=255)
                                                    cAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTACTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTTCCAGCGCAGCTATATCAATATTACACGCACTGCTTTCGACTTGGGGATgt                                                                                                                                                                                        >  2:70273/1‑250 (MQ=255)
                                                     aacaTGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGAACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCCCCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGATAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCCTCGCTCACGCGGAAATACCAGCGCAGCTATATCAATATTACTCGCCCTGCTTTCAGCTTGTGGATgc                                                                                                                                                                                        >  2:405337/4‑249 (MQ=255)
                                                               ggATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATc                                                                                                                                                                             <  2:182225/250‑1 (MQ=255)
                                                                              aaaaGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTg                                                                                                                                                                >  1:170726/1‑249 (MQ=255)
                                                                                aaGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGt                                                                                                                                                            >  1:571591/1‑251 (MQ=255)
                                                                                  gCACCAGAGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGAGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGAGTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCCATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCACTGATCCCCGCCAGCGTTGAGGttt                                                                                                                                                          >  1:484993/1‑251 (MQ=255)
                                                                                  cacccaggtgggTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCTGCGATCTTTTATTCGCGAGGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCATATTTCCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAAGATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTAGAGGttg                                                                                                                                                          <  2:447226/241‑2 (MQ=255)
                                                                                    aCCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTAt                                                                                                                                                        >  1:378440/1‑251 (MQ=255)
                                                                                          gtgGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGAc                                                                                                                                                  >  1:170489/1‑251 (MQ=255)
                                                                                                   aaGGCTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGCTCCTTTATTCGCGATGGGCTGACGCTTCCTGCCCATGTGATCGCTCATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGAATCGCTCCCCCGTCAATACCAGCGAAGCTATAACAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGCTTATCACGACAGCCTGTGc                                                                                                                                         >  1:36179/1‑251 (MQ=255)
                                                                                                           agcaaccccggggggggTCTGCTGGCCCCGGACGCGATGCGTTCTTTTATTCCCGATGGGCTGCCGCTTCCGTCCGATGTGATCCCTAATCCGCAGATTTTCAGCGTCAAAACCGTAGACGTCGCCGCATCGCTACCGCGTAACTCCCAGCGAAGTTATATCAATATTAATCGCATTGCTTTCGACTTGTGGATGAAATCGCCGATCACCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAAt                                                                                                                                 <  2:150353/236‑1 (MQ=255)
                                                                                                              aaaCCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCg                                                                                                                                                                         <  1:697246/208‑1 (MQ=255)
                                                                                                              aaaCCCTGCGGCGGTCTGCTGGCAACGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGAGAAGCTATATCAATATTAATCGCCCTGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTg                                                                                                                              >  1:441032/1‑251 (MQ=255)
                                                                                                                       ggcggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGCTGAAATCGCAGATCCCAGCCAGCGTTGAGGTTTATCACGCCAGCCTGTGACGAAAAATATGGCGTggga                                                                                                                      >  2:563475/1‑247 (MQ=255)
                                                                                                                          ggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGgatgt                                                                                                                  >  2:679570/1‑250 (MQ=255)
                                                                                                                          ggTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGgatg                                                                                                                   >  1:174010/1‑250 (MQ=255)
                                                                                                                                 tgGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATgg                                                                                                           >  1:531817/1‑251 (MQ=255)
                                                                                                                                        ggACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCtt                                                                                                                                                                                                       >  2:533879/1‑152 (MQ=255)
                                                                                                                                        ggACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAACACTGGCGTGCGGATGCTAAATGGCATGTc                                                                                                     >  2:278096/1‑250 (MQ=255)
                                                                                                                                        ggACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATAGTAATCGCCATGCtt                                                                                                                                                                                                       <  1:533879/152‑1 (MQ=255)
                                                                                                                                        ggACGCGCAGCGATCATTTATCTGCGATGGGCTGATGCTTCCGGGCGAGTGGATCCCTAATCCGCATATTTTCAGCGACAAAACCGTCGACGCCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAACCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCa                                                                                                    <  2:129870/251‑1 (MQ=255)
                                                                                                                                         gACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCAt                                                                                                   >  1:173259/1‑251 (MQ=255)
                                                                                                                                                aGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGTTAAATGGCATGTAATTTTTCgg                                                                                            >  2:502779/1‑250 (MQ=255)
                                                                                                                                                  cGATTTTTTATTCGCGATGGGGTGACGCTTCCTGTCGATTGGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACGACCAGCGAAGCTATATCAATATTAATCGCCATGCTGTCGACTTGTGGAGGAAATCGCTGACCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGc                                                                                          <  1:405337/251‑1 (MQ=255)
                                                                                                                                                    aTCTTTTATTCGCGAAGGGCTGACGCTTCCTGTCGAAGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGACGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATTTCATTTTTCGTGccg                                                                                        >  2:390094/1‑249 (MQ=255)
                                                                                                                                                        tttATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGc                                                                                          <  1:521299/245‑1 (MQ=255)
                                                                                                                                                        tttATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACg                                                                                     <  1:319493/250‑1 (MQ=255)
                                                                                                                                                              ggggaggggCGGACGCTCCCTGTTGATGTGACCGCTAAACCGCAGATTTTCAGCGACAAACCCGCGGAGGGCGCCCCATCGCTAACGCGTATCTACCAGCGAAGCTATATCAATATTAATCGCCATGTTTTCGATTGGGGGATGAAAGCGCTGACCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCAGGTCATTTTTCGTGCAGACGGCTGGGa                                                                              <  2:628186/245‑1 (MQ=255)
                                                                                                                                                                             cTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgc                                                               >  1:612867/1‑251 (MQ=255)
                                                                                                                                                                               cccTGTCGTTTTGATCGCTAAGCCTCAGATTGTCAACTTCAAAACCGCAGCCGTCGCCCCACCGTTACCGCGCAACTACCAGAGAAGTTTAACAATAATTAATCGCCATGCTTCCGACTTTGGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGAGGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgccc                                                             <  2:496321/250‑1 (MQ=255)
                                                                                                                                                                                  tGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCACAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTAACAGCGACGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGACATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAACCTGTGCCGCAAAATCTGGCGTAAGGAGGATAAATGGCCTGTCATTTTTCGTGCAGACGGCTGGGAgaaccataatcccgccagct                                                          >  2:179713/1‑232 (MQ=255)
                                                                                                                                                                                  tGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTAGTGCAGACGGCTGGGAGCAACATATTAccgcccgtt                                                          >  2:385607/1‑249 (MQ=255)
                                                                                                                                                                                      ggATGTGACCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATc                                                       <  1:385607/250‑1 (MQ=255)
                                                                                                                                                                                       gATGTGATCGCTACTCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTAGTGCAGACGGCTGGGAGCAACATATTACCCACCGCTATCt                                                      >  2:509042/1‑251 (MQ=255)
                                                                                                                                                                                           tgATCGCTAATCCGCAGCTTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGCCAGCCTCTGCCGAACAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGCCGGCTGGGCGCCACATATTACCTCCCGCTATCTGGTc                                                  >  2:363070/1‑251 (MQ=255)
                                                                                                                                                                                            gATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCg                                                 >  1:75880/1‑251 (MQ=255)
                                                                                                                                                                                            gATCGCTAACCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTACCTACCAGCGAAGCTATACCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTc                                                  <  1:140280/250‑1 (MQ=255)
                                                                                                                                                                                                cGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTg                                              >  1:506264/1‑251 (MQ=255)
                                                                                                                                                                                                   tAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTgacg                                           >  2:485194/1‑248 (MQ=255)
                                                                                                                                                                                                    aaTCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGCGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGGGCCgc                                          >  2:129903/1‑250 (MQ=255)
                                                                                                                                                                                                      tCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATg                                        >  2:12515/1‑251 (MQ=255)
                                                                                                                                                                                                        cGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGGTATATCAATATTAATCGCCAGGCTTTCGACTTGGGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGgc                                      <  1:658894/251‑1 (MQ=255)
                                                                                                                                                                                                        cGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGgc                                      >  1:14406/1‑251 (MQ=255)
                                                                                                                                                                                                             aTTTTCAGCGTCAAAACCGTCGCCGTCGCCGCTTCGCTAACGCGTCACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGCCAGCCTGTGCCGACAAATCTGGCGTGCGGATGATAAATGGCATGTCATTTTTCGTGCATCCGGCTGGGAGCAACATATTACCACCCGCTATCTGGTCGGTGCCGATggagcca                                  >  2:312521/1‑245 (MQ=255)
                                                                                                                                                                                                                ttCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGcc                                                                                                                                                                       >  2:511828/1‑114 (MQ=255)
                                                                                                                                                                                                                ttCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGcc                                                                                                                                                                       <  1:511828/114‑1 (MQ=255)
                                                                                                                                                                                                                         aaaaCCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTgcg                     >  1:510470/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACAtc             <  1:679570/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cATCGCTACCGCGTAACTACCAGCAAAGCTATAATAATATGAACCTCAATGCTTTCGATTGGTGGATAAAATCGCTGACCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGat  <  2:173259/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAACCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGat  <  1:517929/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
CGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGCTTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGAT  >  NC_000913/3869396‑3869880

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