Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,400,874 A→T 6.3% intergenic (+37/‑132) ybgF → / → queE TolQRA transport system periplasmic protein/7‑carboxy‑7‑deazaguanine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,400,8740AT 6.3% 292.2 / 4.7 93intergenic (+37/‑132)ybgF/queETolQRA transport system periplasmic protein/7‑carboxy‑7‑deazaguanine synthase
Reads supporting (aligned to +/- strand):  ref base A (46/41);  new base T (2/4);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTAC  >  NC_002947/1400734‑1401021
                                                                                                                                            |                                                                                                                                                      
tGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCgg                                                                                                                                                 <  6:72226/148‑1 (MQ=255)
  ggCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGtt                                                                                                                                              >  7:90297/1‑149 (MQ=255)
    ccATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGtttt                                                                                                                                            <  5:330522/149‑1 (MQ=255)
     cATACCGACAAGGTAAAAGGCATCCTTGAGCAGGTCGTGACCCAGTACGGCGTCACCTCTGCCGCTCAACTGTCCCAGCTTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTTGAAAGAAACCCGCGCCAGTCGCGGGttttt                                                                                                                                           <  8:300037/149‑1 (MQ=255)
      aTACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCTCTGTACCGCTCGAAAGAAACCCGCGACTGGCGCGGGtttctt                                                                                                                                          <  5:293316/149‑1 (MQ=255)
       tACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTc                                                                                                                                         <  1:369674/149‑1 (MQ=255)
         ccGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGt                                                                                                                                       >  8:306345/1‑149 (MQ=255)
         ccGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAAAAACCCGCGACTGGCGCGGGtttctttc                                                                                                                                        <  5:159038/148‑3 (MQ=255)
          cGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGt                                                                                                                                       >  3:259805/1‑148 (MQ=255)
          cGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGt                                                                                                                                       >  2:158448/1‑148 (MQ=255)
             cAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAg                                                                                                                                    >  8:112933/1‑148 (MQ=255)
              aaGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGaa                                                                                                                                  <  3:280140/149‑1 (MQ=255)
                 gTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGtt                                                                                                                                      >  3:165086/1‑142 (MQ=255)
                 gTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGtt                                                                                                                                      <  4:165086/142‑1 (MQ=255)
                  tAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCAc                                                                                                                              >  7:11751/1‑149 (MQ=255)
                    aaaGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTg                                                                                                                            <  7:306345/149‑1 (MQ=255)
                      aGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGACTGGCGCTGGtttcttt                                                                                                                                         >  5:22959/1‑133 (MQ=255)
                      aGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGACTGGCGCGGGtttcttt                                                                                                                                         <  6:22959/134‑2 (MQ=255)
                            ccTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCttttt                                                                                                                    >  2:307256/1‑149 (MQ=255)
                             cTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCttttt                                                                                                                    <  8:90297/148‑1 (MQ=255)
                              ttCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAA‑AAACCCGCGACTGGCGCGGGtttctttc                                                                                                                                        <  8:5768/126‑3 (MQ=38)
                              ttCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAA‑AAACCCGCGACTGGCGCGGGtttctttc                                                                                                                                        >  7:5768/1‑124 (MQ=38)
                                cagcagGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTc                                                                                                                  <  8:261422/147‑1 (MQ=255)
                                  gcagGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑Taa                                                                                                           <  4:227572/149‑1 (MQ=255)
                                  gcagGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑Ta                                                                                                            >  5:205269/1‑148 (MQ=255)
                                   cagGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑Taaa                                                                                                          >  2:322905/1‑149 (MQ=255)
                                    agGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑Taaa                                                                                                          >  5:390264/1‑148 (MQ=255)
                                         gTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGTGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACgc                                                                                                     <  7:20555/148‑1 (MQ=255)
                                          tGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACgctg                                                                                                   >  7:8137/1‑149 (MQ=255)
                                           gACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑t                                                                                                             >  8:95958/1‑138 (MQ=255)
                                           gACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑t                                                                                                             <  7:95958/138‑1 (MQ=255)
                                           gACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACgctgc                                                                                                  >  6:72777/1‑149 (MQ=255)
                                                aGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCg                                                                                              >  2:81354/1‑148 (MQ=255)
                                                     cccGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAAc                                                                                        >  3:249275/1‑149 (MQ=255)
                                                     cccGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGACCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAAc                                                                                        >  2:333212/1‑149 (MQ=255)
                                                       cGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGc                                                                                      >  6:419599/1‑149 (MQ=255)
                                                            ccTCTGCCGCTCAACTGGCCCAGCGTGACCTGTAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGt                                                                                 <  7:153608/149‑1 (MQ=255)
                                                             ctctGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTc                                                                                >  5:392929/1‑149 (MQ=255)
                                                             ctctGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTc                                                                                >  8:298758/1‑149 (MQ=255)
                                                                  ccGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGa                                                                              <  5:165981/146‑1 (MQ=255)
                                                                  ccGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTc                                                                           <  1:158448/149‑1 (MQ=255)
                                                                    gCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACGCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCcg                                                                         <  4:259805/149‑1 (MQ=255)
                                                                     cTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTGGAGTCcgc                                                                        <  3:388234/149‑1 (MQ=255)
                                                                     cTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCcgc                                                                        >  5:372490/1‑149 (MQ=255)
                                                                     cTCAACTGGCCCACCGAGACCTGCAGAACCTCTAGGCGATGTCGCTGTACCGCTCGAAACAAACCCGCGCCAGACGCGGGTATAATCCATAGAATCACTGCCGTTTATTCG‑‑‑TAAACACGCTCCTGCGGATAACGCTATGTAGTGTCcgt                                                                        >  8:373744/1‑148 (MQ=255)
                                                                        aaCTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGa                                                                      <  6:205269/148‑1 (MQ=255)
                                                                             gCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATa                                                                                          <  4:227021/123‑1 (MQ=255)
                                                                             gCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATa                                                                                          >  3:227021/1‑123 (MQ=255)
                                                                             gCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGc                                                                <  8:8137/149‑1 (MQ=255)
                                                                             gCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGc                                                                <  1:307256/149‑1 (MQ=255)
                                                                             gCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGc                                                                <  1:81354/149‑1 (MQ=255)
                                                                               ccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCt                                                              >  5:24326/1‑149 (MQ=255)
                                                                               ccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCt                                                              >  6:118237/1‑149 (MQ=255)
                                                                                  gCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGAc                                                           >  7:81334/1‑149 (MQ=255)
                                                                                       aCCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGtttttt                                                                                                                                          <  7:284299/68‑1 (MQ=255)
                                                                                       aCCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGtttttt                                                                                                                                          >  8:284299/1‑68 (MQ=255)
                                                                                          tGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCg                                                    >  4:381378/1‑148 (MQ=255)
                                                                                          tGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCg                                                    >  4:388101/1‑148 (MQ=255)
                                                                                           gCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCgg                                                                                             >  6:147680/1‑106 (MQ=255)
                                                                                           gCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCgg                                                                                             <  5:147680/106‑1 (MQ=255)
                                                                                            cAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCgg                                                   <  6:392929/147‑1 (MQ=255)
                                                                                               aaGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACgg                                                         <  3:311638/138‑1 (MQ=255)
                                                                                               aaGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACgg                                                         >  4:311638/1‑138 (MQ=255)
                                                                                                  ctctAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTg                                           <  6:24326/149‑1 (MQ=255)
                                                                                                   tctAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGATATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGt                                          >  6:393376/1‑149 (MQ=255)
                                                                                                     tAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt                                         >  8:286232/1‑148 (MQ=255)
                                                                                                       ggCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt                                         >  2:358384/1‑146 (MQ=255)
                                                                                                       ggCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt                                         <  5:118237/146‑1 (MQ=255)
                                                                                                       ggCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTa                                       <  5:72777/148‑1 (MQ=255)
                                                                                                       ggCGTTGTAGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTATCATCACTGCCCTTTTTTCGCCCTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTg                                           <  1:333212/147‑1 (MQ=255)
                                                                                                        gCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGc                                     >  1:98492/1‑149 (MQ=255)
                                                                                                         cGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCcg                                                                         >  1:380856/1‑112 (MQ=255)
                                                                                                         cGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCcg                                                                         <  2:380856/112‑1 (MQ=255)
                                                                                                          gTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTg                                   >  7:298871/1‑149 (MQ=255)
                                                                                                                    tACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGt                         >  2:357002/1‑149 (MQ=255)
                                                                                                                         cTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGc                      >  1:252670/1‑147 (MQ=255)
                                                                                                                          tCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAg                   >  4:111747/1‑149 (MQ=255)
                                                                                                                          tCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAg                   >  4:228576/1‑149 (MQ=255)
                                                                                                                            gaaagaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCa                 >  3:28583/1‑149 (MQ=255)
                                                                                                                             aaagaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATACCGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCAt                <  7:112933/149‑1 (MQ=255)
                                                                                                                             aaagaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCAt                <  3:388101/149‑1 (MQ=255)
                                                                                                                             aaagaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCAt                <  1:357002/149‑1 (MQ=255)
                                                                                                                               agaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGc              <  6:390264/149‑1 (MQ=255)
                                                                                                                               agaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGc              <  6:372490/149‑1 (MQ=255)
                                                                                                                               agaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGc              <  7:286232/149‑1 (MQ=255)
                                                                                                                                 aaaCCCTCGCCAGGCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCTGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCaa            <  6:107807/149‑1 (MQ=255)
                                                                                                                                 aaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCaa            >  6:124533/1‑149 (MQ=255)
                                                                                                                                  aaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCaa            <  4:249275/148‑1 (MQ=255)
                                                                                                                                     ccGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt                                         <  4:246506/116‑1 (MQ=255)
                                                                                                                                     ccGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt                                         >  3:246506/1‑116 (MQ=255)
                                                                                                                                       gcgcCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGAcaca      >  8:6773/1‑149 (MQ=255)
                                                                                                                                        cgcCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACAt     <  1:358384/149‑1 (MQ=255)
                                                                                                                                           cAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTAc  >  1:220034/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                      
TGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTAC  >  NC_002947/1400734‑1401021

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