Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,169,260 Δ1 bp 8.5% intergenic (‑29/+24) ygiV ← / ← ygiW transcriptional repressor for mcbR biofilm gene/hydrogen peroxide and cadmium resistance periplasmic protein; stress‑induced OB‑fold protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,169,2600A. 8.5% 276.8 / 19.7 94intergenic (‑29/+24)ygiV/ygiWtranscriptional repressor for mcbR biofilm gene/hydrogen peroxide and cadmium resistance periplasmic protein; stress‑induced OB‑fold protein
Reads supporting (aligned to +/- strand):  ref base A (49/36);  new base . (3/5);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-01

CATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTGG  >  NC_000913/3169115‑3169407
                                                                                                                                                 |                                                                                                                                                     
caTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAAcg                                                                                                                                                    >  1:269031/1‑149 (MQ=255)
 aTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAAcgc                                                                                                                                                   <  2:65736/149‑1 (MQ=255)
     tAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTcc                                                                                                                                             >  1:151510/1‑149 (MQ=255)
      aaaTTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTccc                                                                                                                                            >  1:906945/1‑149 (MQ=255)
       aaTTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCg                                                                                                                                           <  1:852221/149‑1 (MQ=255)
       aaTTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATGTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCg                                                                                                                                           >  1:909579/1‑149 (MQ=255)
         tttCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGACGCGTTCCCggg                                                                                                                                         <  2:316029/149‑1 (MQ=255)
           tCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGTCTCCCTGCGAAATTATTCCCGGGACGCGTTCCCGGGaa                                                                                                                                       <  2:106457/149‑2 (MQ=255)
           tCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAg                                                                                                                                       <  1:520637/149‑1 (MQ=255)
            cGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGACGCGTTCCCGGGaat                                                                                                                                      <  2:397740/149‑3 (MQ=255)
            cGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGACGCGTTCCCGGGaat                                                                                                                                      <  1:531468/149‑3 (MQ=255)
              ccACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGACGCGTTCCCGGGaataa                                                                                                                                    <  1:364252/149‑5 (MQ=255)
               cccGCTGTAGTTCAGCAATTCAGGGCTACAGCGGGGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGt                                                                                                                                   <  2:1046371/147‑1 (MQ=255)
               cACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGt                                                                                                                                   <  2:417375/149‑1 (MQ=255)
               cACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGt                                                                                                                                   <  1:849745/149‑1 (MQ=255)
                 cGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTaa                                                                                                                                 <  1:1002628/149‑1 (MQ=255)
                  gCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAAc                                                                                                                                <  1:1229153/149‑1 (MQ=255)
                  gCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAAc                                                                                                                                <  1:172772/149‑1 (MQ=255)
                     gTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAAt                                                                                                                             >  1:284249/1‑149 (MQ=255)
                        gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTAc                                                                                                                          >  1:395161/1‑149 (MQ=255)
                        gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGGCCCGGGAGCGGTAACAATTAc                                                                                                                          <  2:642837/149‑1 (MQ=255)
                           gAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGa                                                                                                                       >  2:451330/1‑149 (MQ=255)
                            aGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGAt                                                                                                                      >  2:994867/1‑149 (MQ=255)
                                aTTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTAc                                                                                                                  >  2:874309/1‑149 (MQ=255)
                                 ttCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACt                                                                                                                 >  1:984173/1‑149 (MQ=255)
                                    aGGGCTACAGCGGTGCGGCAAGATCGCCACAGGTATTGATGGGAAATCGATAATGGTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACtttt                                                                                                              <  2:676320/149‑2 (MQ=255)
                                        cTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGTGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGa                                                                                                          >  1:786857/1‑149 (MQ=255)
                                         tACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTGATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGTGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGAt                                                                                                         >  2:847570/1‑149 (MQ=255)
                                           cAGCGGTGCGGCAACATCGCCACATGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGCACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCt                                                                                                       <  1:1064658/149‑1 (MQ=255)
                                           cAGCGGGGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCATGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCt                                                                                                       <  2:565951/149‑1 (MQ=255)
                                                gTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCATCCCAGAGAAATTATTCCCGGCAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTTCGGATCTGTTTg                                                                                                  >  1:649696/1‑149 (MQ=255)
                                                   cGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACg                                                                                               >  1:330403/1‑149 (MQ=255)
                                                    ggCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGGGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGt                                                                                              <  2:26300/149‑1 (MQ=255)
                                                     gCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTc                                                                                             >  2:471460/1‑149 (MQ=255)
                                                      cAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCa                                                                                            <  1:173303/149‑1 (MQ=255)
                                                       aaCATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCaa                                                                                           <  2:574391/149‑1 (MQ=255)
                                                        aCATCGCCACACGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCAATTTGCCTCGCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAt                                                                                          <  1:374419/149‑1 (MQ=255)
                                                          aTCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCAGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAttt                                                                                        >  1:1067553/1‑149 (MQ=255)
                                                               cacaGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGGCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGAGCTGTTTGACGTCAATTTCAACa                                                                                   <  1:559739/149‑1 (MQ=255)
                                                                     tATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTc                                                                             >  2:606681/1‑149 (MQ=255)
                                                                       ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCa                                                                           <  2:1006607/149‑1 (MQ=255)
                                                                       ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCa                                                                           <  2:151510/149‑1 (MQ=255)
                                                                         gATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGt                                                                         >  2:649998/1‑149 (MQ=255)
                                                                          aTGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTc                                                                        >  2:419507/1‑149 (MQ=255)
                                                                              gAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTa                                                                    >  2:188053/1‑149 (MQ=255)
                                                                                aaTCGATAATGTTTACATCGAGTGTCAGGGTTGTCATTTTGCCTCCCTGCGTAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAAATACGGATATACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATc                                                                  <  2:142023/149‑1 (MQ=255)
                                                                                    gATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGCCGTCAATTTCAACAGAATTCCAGTCTTTATCGACt                                                              >  2:796100/1‑149 (MQ=255)
                                                                                         tGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCAcc                                                         >  2:1320566/1‑149 (MQ=255)
                                                                                         tGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTGTGACGTCAATTTCAACAGAATTCCAGTCTGTATCGACTTCAcc                                                         >  2:430477/1‑149 (MQ=255)
                                                                                         tGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGAGCTGTGTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCAcc                                                         >  2:680774/1‑149 (MQ=255)
                                                                                         tGTTTACATCCAGTGTCAGGTTTGTCATTATGCCTCCCTGCGAAAATATGACCGGGAACGCGGCCCGGTAGCGGTAACAATTACGGATTTACTTTGCGGATATGTTTGACGTCAATTTCAACCGAATTCCAGTCTTTATCGACTTCAcc                                                         <  2:1230435/149‑1 (MQ=255)
                                                                                          gTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCAccc                                                        >  1:549064/1‑149 (MQ=255)
                                                                                            ttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCt                                                       <  1:1025852/149‑1 (MQ=255)
                                                                                              aCATCCTGTGTCTGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATATGTTTGACGTCAATTACAACAGAATTCCAGTCTTTATCGCCTTCACCCTGaa                                                    >  2:675844/1‑149 (MQ=255)
                                                                                               cATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAt                                                   <  2:1019488/149‑1 (MQ=255)
                                                                                               cATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAt                                                   <  1:388104/149‑1 (MQ=255)
                                                                                                aTCCAGTGTCAGGTTTGTCATTTTGCATCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGAGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtc                                                  >  1:831810/1‑149 (MQ=255)
                                                                                                 tCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtct                                                 >  2:1345475/1‑149 (MQ=255)
                                                                                                 tCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTAGACGACAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtct                                                 >  1:1034448/1‑149 (MQ=255)
                                                                                                     gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                             >  2:416568/1‑149 (MQ=255)
                                                                                                       gtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCg                                           >  1:1104255/1‑149 (MQ=255)
                                                                                                         cAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTa                                         <  2:984173/149‑1 (MQ=255)
                                                                                                          aGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGGCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTAt                                        <  1:874309/149‑1 (MQ=255)
                                                                                                            gTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCt                                      >  1:542632/1‑149 (MQ=255)
                                                                                                            gTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCt                                      >  2:357595/1‑149 (MQ=255)
                                                                                                                gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACTTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTCCg                                  <  2:786857/149‑1 (MQ=255)
                                                                                                                 tCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCgg                                 >  2:156313/1‑149 (MQ=255)
                                                                                                                    ttttGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgt                              >  2:340660/1‑149 (MQ=255)
                                                                                                                     tttGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtc                             <  2:154828/149‑1 (MQ=255)
                                                                                                                     tttGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtc                             <  2:477749/149‑1 (MQ=255)
                                                                                                                      ttGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtca                            <  2:986828/149‑1 (MQ=255)
                                                                                                                      ttGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtca                            >  2:964545/1‑149 (MQ=255)
                                                                                                                       tGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcac                           <  1:451330/149‑1 (MQ=255)
                                                                                                                        gCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTTTCTTTCGGcgtcacc                          >  2:395140/1‑149 (MQ=255)
                                                                                                                           tCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgtc                       >  1:576236/1‑149 (MQ=255)
                                                                                                                           tCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgtc                       >  1:640613/1‑149 (MQ=255)
                                                                                                                               tGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCCCACTGAATCTCAACCGTTTCTTTCGGCGTCACCGTCACGc                   >  1:693057/1‑149 (MQ=255)
                                                                                                                               tGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTAACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGc                   <  2:833807/149‑1 (MQ=255)
                                                                                                                               tGCGAAATTATTACCGGGAACGCGTCCCGGGAGCGGTAACAATTACTGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCATTATCGACTTCAACCTGAATCTCAACCGAATCTTTCGGCATCACCGTCACGc                   >  2:611906/1‑149 (MQ=255)
                                                                                                                                 cGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCg                 >  1:97825/1‑149 (MQ=255)
                                                                                                                                  gAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGt                <  2:524752/149‑1 (MQ=255)
                                                                                                                                  gAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCAGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGt                <  2:627004/149‑1 (MQ=255)
                                                                                                                                   aaaTTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGAGTCACCGTCACGCCGtt               >  1:997240/1‑149 (MQ=255)
                                                                                                                                     attattCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTcc             >  1:523771/1‑149 (MQ=255)
                                                                                                                                     attattCCCGGGAACGCGTCCCGGGAGCGGTAACAAGTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCACCCGTATCTTTCGGCGTCACCGTCACGCCGGTcc             >  1:528747/1‑149 (MQ=255)
                                                                                                                                       cgctCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCCAACGTATCTTTACGCGTCACCGTCACGCCGTTACAg           >  2:781058/4‑149 (MQ=255)
                                                                                                                                       cgctCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAg           >  2:621717/4‑149 (MQ=255)
                                                                                                                                        gctCCCGGGACCCCGTCCCGGGAGCGGTAAAAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAAAAAAATTCCAGTATTTTTAGACTTAACCATGAAAAAAAACCGTATCTTAAGGCGACCCCGACACGCCGTTCCAGc          >  2:703634/3‑149 (MQ=255)
                                                                                                                                        gctCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGc          >  2:25684/3‑149 (MQ=255)
                                                                                                                                        attCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATATCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGc          <  1:1261856/149‑1 (MQ=255)
                                                                                                                                            ccGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGACTCTCAACCGTATCTTTCGGCGTCCCCGTCACGCCGTTCCAGCGttt      >  2:1106535/1‑149 (MQ=255)
                                                                                                                                            ccGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt      >  2:8371/1‑149 (MQ=255)
                                                                                                                                             cGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTtg     >  1:322251/1‑149 (MQ=255)
                                                                                                                                                gAACGTGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTgg  <  1:847570/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTGG  >  NC_000913/3169115‑3169407

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