Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,811,185 T→C L139P (CTC→CCC)  mprA → transcriptional repressor of microcin B17 synthesis and multidrug efflux

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,811,1850TC98.9% 346.2 / NA 94L139P (CTC→CCC) mprAtranscriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (44/49);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTAT  >  NC_000913/2811037‑2811331
                                                                                                                                                    |                                                                                                                                                  
cgTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAAccct                                                                                                                                                  <  1:578711/151‑1 (MQ=255)
    ttGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTc                                                                                                                                              >  2:863987/1‑151 (MQ=255)
       ccGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCcgc                                                                                                                                           <  2:24467/151‑1 (MQ=255)
       ccGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCcgc                                                                                                                                           >  1:798504/1‑151 (MQ=255)
        cGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCcgcg                                                                                                                                          <  2:641499/151‑1 (MQ=255)
          aTGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTg                                                                                                                                                 >  1:71859/1‑142 (MQ=255)
          aTGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTg                                                                                                                                                 <  2:71859/142‑1 (MQ=255)
           tGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTc                                                                                                                                       >  2:218263/1‑151 (MQ=255)
           tGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTc                                                                                                                                       >  1:347765/1‑151 (MQ=255)
            gAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTc                                                                                                                                       >  2:296354/1‑150 (MQ=255)
            gAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTc                                                                                                                                       <  2:555677/150‑1 (MQ=255)
             aaCTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATACCTGCCTGCATCAACCCTGGTCCGCGCTCAg                                                                                                                                     <  1:275889/151‑1 (MQ=255)
              aCTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGc                                                                                                                                    <  2:249489/151‑1 (MQ=255)
               cTGGAAACACGCGGTTGGATCGCACGTCGTGAAAGCGATCACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCATCATAACTGCCTGCATCAACCCTGGTCCGCGCTc                                                                                                                                       >  1:603513/1‑147 (MQ=255)
                tGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGcac                                                                                                                                  >  1:709054/1‑151 (MQ=255)
                 ggAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGcaca                                                                                                                                 <  2:266767/151‑1 (MQ=255)
                 ggAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGcaca                                                                                                                                 <  1:726815/151‑1 (MQ=255)
                  gAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGcaca                                                                                                                                 >  1:662918/1‑150 (MQ=255)
                   aaaaaCGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCacaac                                                                                                                               <  2:1131566/151‑1 (MQ=255)
                        cgcgGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGaa                                                                                                                           <  1:1183233/150‑1 (MQ=255)
                          cgGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGaaaa                                                                                                                         >  1:377125/1‑150 (MQ=255)
                           ggTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGaaaa                                                                                                                         >  1:170171/1‑149 (MQ=255)
                              tGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGaaaa                                                                                                                         >  2:171508/1‑146 (MQ=255)
                              tGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGAt                                                                                                                     <  1:616642/150‑1 (MQ=255)
                               ggATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCa                                                                                                                   <  2:432724/151‑1 (MQ=255)
                                     aaCGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCg                                                                                                              >  1:642208/1‑150 (MQ=255)
                                      aCGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAg                                                                                                            <  1:154578/151‑1 (MQ=255)
                                        gtcgtGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCa                                                                                                          <  1:953231/151‑1 (MQ=255)
                                          cgtGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCaaa                                                                                                        <  2:825497/151‑1 (MQ=255)
                                          cgtGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCaa                                                                                                         <  2:1092209/150‑1 (MQ=255)
                                           gtGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCCATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAAt                                                                                                       >  2:249751/1‑151 (MQ=255)
                                           gtGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAAt                                                                                                       <  2:63814/151‑1 (MQ=255)
                                            tGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATc                                                                                                      >  2:387936/1‑151 (MQ=255)
                                             gAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCa                                                                                                     <  2:25335/151‑1 (MQ=255)
                                                aGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCAccc                                                                                                  >  1:644325/1‑151 (MQ=255)
                                                aGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCAccc                                                                                                  >  1:561338/1‑151 (MQ=255)
                                                aGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCAccc                                                                                                  <  1:16919/151‑1 (MQ=255)
                                                    aTAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCa                                                                                               >  1:940469/1‑150 (MQ=255)
                                                       aCGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAAtt                                                                                           <  1:1175630/151‑1 (MQ=255)
                                                       aCGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAAtt                                                                                           <  1:516407/151‑1 (MQ=255)
                                                       aCGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAAtt                                                                                           <  1:246503/151‑1 (MQ=255)
                                                               cgcTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTccc                                                                                   >  2:425828/1‑151 (MQ=255)
                                                                gcTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCg                                                                                  >  2:457474/1‑151 (MQ=255)
                                                                  tgcctgcATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGtc                                                                                >  2:414813/1‑151 (MQ=255)
                                                                   gcctgcATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGtct                                                                               <  1:36957/151‑1 (MQ=255)
                                                                     ctgcATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCg                                                                             >  2:683892/1‑151 (MQ=255)
                                                                     ctgcATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCg                                                                             >  1:930923/1‑151 (MQ=255)
                                                                       gcATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGa                                                                            <  2:377408/150‑1 (MQ=255)
                                                                        cATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGAcc                                                                          <  1:686053/151‑1 (MQ=255)
                                                                          tCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGAcc                                                                          >  2:37626/1‑149 (MQ=255)
                                                                             gCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAg                                                                        <  2:1073904/148‑1 (MQ=255)
                                                                              cAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATgg                                                                    <  2:436811/151‑1 (MQ=255)
                                                                                   aaCGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACa                                                                >  1:801156/1‑150 (MQ=255)
                                                                                        aaaaaGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACgg                                                          >  1:18487/1‑151 (MQ=255)
                                                                                           aaGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGgtgt                                                       >  1:800708/1‑151 (MQ=255)
                                                                                           aaGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGgtgt                                                       >  2:491748/1‑151 (MQ=255)
                                                                                            aGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGtgtg                                                      >  1:376271/1‑151 (MQ=255)
                                                                                              gTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGt                                                    >  2:980607/1‑151 (MQ=255)
                                                                                                  cGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTATACCGTCTCGACCAGATGGAACAAGACGGTGTGGTtctc                                                <  1:192814/151‑1 (MQ=255)
                                                                                                   gAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTtctca                                               <  1:442752/151‑2 (MQ=255)
                                                                                                    aGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCg                                               >  1:846331/1‑150 (MQ=255)
                                                                                                     gTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGa                                              >  1:302080/1‑150 (MQ=255)
                                                                                                     gTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGa                                              <  2:742900/150‑1 (MQ=255)
                                                                                                      tttttGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGa                                              >  2:922665/1‑149 (MQ=255)
                                                                                                         ttGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCg                                          >  2:403904/1‑150 (MQ=255)
                                                                                                           gcgcGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGAt                                        <  2:900887/150‑1 (MQ=255)
                                                                                                              cgAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGc                                    >  2:757658/1‑151 (MQ=255)
                                                                                                                   ttttACCACCGCAGCATAGCTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAAcg                               >  2:1120740/1‑151 (MQ=255)
                                                                                                                   ttttACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAAcg                               >  1:1144599/1‑151 (MQ=255)
                                                                                                                    tttACCACCGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAAcg                               <  1:552915/150‑1 (MQ=255)
                                                                                                                      tcccaccGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAAc                                <  2:368116/145‑1 (MQ=255)
                                                                                                                        ccaccGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCa                          <  2:685478/151‑1 (MQ=255)
                                                                                                                         caccGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCAt                         <  2:297026/151‑1 (MQ=255)
                                                                                                                           ccGCAGCATAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCAtct                       <  1:893839/151‑1 (MQ=255)
                                                                                                                              cagcaTAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGc                    >  2:720894/1‑151 (MQ=255)
                                                                                                                               agcaTAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCt                   >  1:411909/1‑151 (MQ=255)
                                                                                                                                gcaTAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTc                  >  1:1016085/1‑151 (MQ=255)
                                                                                                                                gcaTAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTc                  <  2:1185953/151‑1 (MQ=255)
                                                                                                                                gcaTAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTc                  <  1:572510/151‑1 (MQ=255)
                                                                                                                                 caTAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCa                 <  1:678170/151‑1 (MQ=255)
                                                                                                                                  aTAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCa                 >  1:867008/1‑150 (MQ=255)
                                                                                                                                   tAACTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCaa                >  2:830888/1‑150 (MQ=255)
                                                                                                                                    aaCTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCaaaa              <  2:900512/151‑1 (MQ=255)
                                                                                                                                     aCTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCaaaaa             <  1:62379/151‑1 (MQ=255)
                                                                                                                                     aCTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCaaaaa             <  1:672909/151‑1 (MQ=255)
                                                                                                                                     aCTGCCTGCATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCaaaaa             <  1:190890/151‑1 (MQ=255)
                                                                                                                                       tgcctgcATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATc           >  2:380210/1‑151 (MQ=255)
                                                                                                                                        gcctgcATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCCGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATc           <  1:83022/150‑1 (MQ=255)
                                                                                                                                        gcctgcATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAAt            <  1:965326/149‑1 (MQ=255)
                                                                                                                                             cATCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGAtt     <  2:54351/151‑1 (MQ=255)
                                                                                                                                              aTCAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAGTTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGAtt     <  1:191495/150‑1 (MQ=255)
                                                                                                                                                cAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCTTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTta   <  1:258026/150‑1 (MQ=255)
                                                                                                                                                cAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTtat  <  2:333512/151‑1 (MQ=255)
                                                                                                                                                cAACCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTtat  <  2:1025878/151‑1 (MQ=255)
                                                                                                                                                 aaCCCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTtat  >  2:1129074/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACAAGACGGTGTGGTTCTCGAAGCGATGAGCTAACGCGTCATCTCGCTCAAAAATCCAGATTTAT  >  NC_000913/2811037‑2811331

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