Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,400,872 C→A 6.2% intergenic (+35/‑134) 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,8720CA 6.2% 290.1 / 6.0 90intergenic (+35/‑134)ybgF/queETolQRA transport system periplasmic protein/7‑carboxy‑7‑deazaguanine synthase
Reads supporting (aligned to +/- strand):  ref base C (44/40);  new base A (2/4);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACAT  >  NC_002947/1400734‑1401018
                                                                                                                                          |                                                                                                                                                     
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                                                                                                                                    >  2:158448/1‑148 (MQ=255)
          cGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGt                                                                                                                                    >  3:259805/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)
                              ttCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAAAACCCGCGACTGGCGCGGGtttctttc                                                                                                                                     >  7:5768/1‑124 (MQ=38)
                              ttCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAAAACCCGCGACTGGCGCGGGtttctttc                                                                                                                                     <  8:5768/126‑3 (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                                                                                                          <  7:95958/138‑1 (MQ=255)
                                           gACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑t                                                                                                          >  8:95958/1‑138 (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                                                             <  1:81354/149‑1 (MQ=255)
                                                                             gCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGc                                                             <  8:8137/149‑1 (MQ=255)
                                                                             gCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGc                                                             <  1:307256/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                                                      >  4:311638/1‑138 (MQ=255)
                                                                                               aaGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACgg                                                      <  3:311638/138‑1 (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                                      <  5:118237/146‑1 (MQ=255)
                                                                                                       ggCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt                                      >  2:358384/1‑146 (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                                                                      <  2:380856/112‑1 (MQ=255)
                                                                                                         cGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCcg                                                                      >  1:380856/1‑112 (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:228576/1‑149 (MQ=255)
                                                                                                                          tCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAg                >  4:111747/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:372490/149‑1 (MQ=255)
                                                                                                                               agaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGc           <  7:286232/149‑1 (MQ=255)
                                                                                                                               agaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGc           <  6:390264/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)
                                                                                                                                          |                                                                                                                                                     
TGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCG‑‑‑TAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACAT  >  NC_002947/1400734‑1401018

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: