Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,489,370 A→G 8.3% intergenic (+169/‑34) ybgI → / → cysD metal‑binding hydrolase/oxidase/sulfate adenylyltransferase subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,489,3700AG 8.3% 237.3 / 7.1 84intergenic (+169/‑34)ybgI/cysDmetal‑binding hydrolase/oxidase/sulfate adenylyltransferase subunit 2
Reads supporting (aligned to +/- strand):  ref base A (42/35);  new base G (2/5);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTCATAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGTACTCGATCGGC  >  NC_002947/1489226‑1489511
                                                                                                                                                 |                                                                                                                                               
gCCTCATAGGAGCGGGTTCCCCCGCGACGAGGCCGGCCCTGCCAGACCAAACCCCCAGTCATATCGTGAGCCTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGcc                                                                                                                                             >  1:372633/1‑149 (MQ=255)
  cTCATAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAgg                                                                                                                                               >  4:179105/1‑145 (MQ=255)
  cTCATAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAgg                                                                                                                                               <  3:179105/145‑1 (MQ=255)
  cTCATAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGt                                                                                                                                           <  6:249390/149‑1 (MQ=255)
     aTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTc                                                                                                                                          <  6:356533/147‑1 (MQ=255)
     aTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCa                                                                                                                                        >  1:162277/1‑149 (MQ=255)
     aTAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCa                                                                                                                                        >  1:86086/1‑149 (MQ=255)
      tAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCTGCGGGCCGtgttc                                                                                                                                       <  4:328074/149‑5 (MQ=255)
      tAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCTGCGGGCCGtgttc                                                                                                                                       <  8:153734/149‑5 (MQ=255)
        ggAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTATACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATaa                                                                                                                                     <  8:350255/149‑1 (MQ=255)
        ggAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGAATATTCGTTAGATTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCTGCGGGCCGtgttcg                                                                                                                                      <  1:334428/148‑6 (MQ=255)
         gAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGt                                                                                                                                           <  2:272815/142‑1 (MQ=255)
         gAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGt                                                                                                                                           >  1:272815/1‑142 (MQ=255)
             gggTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCg                                                                                                                                <  6:139349/149‑1 (MQ=255)
                   aCCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTc                                                                                                                          >  1:81470/1‑149 (MQ=255)
                   aCCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGAGCGTTTTTc                                                                                                                          >  3:6533/1‑149 (MQ=255)
                     ccGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCGGGCCGTCCATAAGATCGTTTTTCGt                                                                                                                        >  2:240063/1‑149 (MQ=255)
                      cgcgAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTg                                                                                                                       <  3:1599/149‑1 (MQ=255)
                       gcgAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGa                                                                                                                      <  1:387960/149‑1 (MQ=255)
                       gcgAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGa                                                                                                                      <  1:244544/149‑1 (MQ=255)
                          aaGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAg                                                                                                                     >  1:25009/1‑147 (MQ=255)
                            gagGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGc                                                                                                                 <  1:367928/149‑1 (MQ=255)
                             agGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGc                                                                                                                 >  6:238949/1‑148 (MQ=255)
                                 cggccCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTCTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGtt                                                                                                                              <  2:152465/131‑1 (MQ=255)
                                 cggccCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTATTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGtt                                                                                                                              >  1:152465/1‑131 (MQ=255)
                                 cggccCTGCCAGCCCAAACCCCCAGCCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATgg                                                                                                            >  6:235177/1‑149 (MQ=255)
                                          cAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATTAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAAc                                                                                                 >  8:190444/1‑149 (MQ=255)
                                           ggCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACt                                                                                                >  5:12912/2‑149 (MQ=255)
                                              ccAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGAc                                                                                             >  3:12092/1‑149 (MQ=255)
                                              ccAAACCCCCAGCCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGAc                                                                                             <  5:235177/149‑1 (MQ=255)
                                               cAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACg                                                                                            >  4:205008/1‑149 (MQ=255)
                                                   cccccAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACt                                                                                        <  7:327676/149‑1 (MQ=255)
                                                     cccAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTg                                                                                                                       <  6:197157/118‑1 (MQ=255)
                                                     cccAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTg                                                                                                                       >  5:197157/1‑118 (MQ=255)
                                                     cccAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACtt                                                                                       >  5:359103/1‑148 (MQ=255)
                                                        aGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCGGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGaaa                                                                                   <  1:240063/149‑1 (MQ=255)
                                                        aGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGaaa                                                                                   <  6:12912/149‑1 (MQ=255)
                                                            atatCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGc                                                                               <  5:229021/149‑1 (MQ=255)
                                                              atCGTTAGTCTTTTTCGATCTAGCCAGCCTCCTAAATAGTAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCTTGGTCGACAAACTGACGCACTTGAAACAGCt                                                                              <  6:264259/148‑1 (MQ=255)
                                                                tCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTg                                                                             <  6:203371/148‑1 (MQ=255)
                                                                  gTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGa                                                                           <  2:25009/148‑1 (MQ=255)
                                                                       aCTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCg                                                                       >  6:167910/1‑147 (MQ=255)
                                                                       aCTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGga                                                                     <  4:296476/149‑1 (MQ=255)
                                                                           tttCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGagag                                                                  >  5:111622/1‑148 (MQ=255)
                                                                              cGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTGGAAACAGCTGGAGGCGGAGAGCATc                                                              >  8:87952/1‑149 (MQ=255)
                                                                              cGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCCCTCGAACACGGCCCGCAGGCAGTCCATAAGATCGTTTTTCGTGAGTAGCAATGGTCGACAATCTGACGCACTTCAAACAGCTCGACGCGCAGAGCATc                                                              >  5:397240/1‑149 (MQ=255)
                                                                               gATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATc                                                              <  2:198213/148‑1 (MQ=255)
                                                                                    agccagccTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCAcat                                                         >  7:146338/1‑148 (MQ=255)
                                                                                     gccagccTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCcaca                                                          <  3:205008/146‑1 (MQ=255)
                                                                                           ctcctAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCcgcg                                                 <  8:146338/149‑1 (MQ=255)
                                                                                             cctAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAg                                               >  3:171018/1‑149 (MQ=255)
                                                                                              ctAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGa                                                <  2:86086/147‑1 (MQ=255)
                                                                                               tAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGt                                             >  1:215506/1‑149 (MQ=255)
                                                                                                 aaTAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTgg                                           >  5:287599/1‑149 (MQ=255)
                                                                                                 aaTAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCAGCGAGGTgg                                           >  8:4918/1‑149 (MQ=255)
                                                                                                 aaTAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGgc                                          >  7:305663/1‑149 (MQ=255)
                                                                                                   tAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATAGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGgcc                                         >  7:329967/1‑149 (MQ=255)
                                                                                                          agcTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGa                                    <  5:248954/147‑1 (MQ=255)
                                                                                                          agcTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGa                                    >  6:248954/1‑147 (MQ=255)
                                                                                                           gctgctGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCACGAGGTGGCCGCCGAGtt                                 >  7:50470/1‑149 (MQ=255)
                                                                                                           gctgctGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCACGAGGTGGCCGCCGAGtt                                 >  7:56438/1‑149 (MQ=255)
                                                                                                             tgctgTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCcac                                                           <  2:105349/121‑1 (MQ=255)
                                                                                                             tgctgTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCcac                                                           >  1:105349/1‑121 (MQ=255)
                                                                                                               ctgTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGa                              >  3:169189/1‑148 (MQ=255)
                                                                                                                  tgATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAAc                                                                                  <  8:73362/93‑1 (MQ=255)
                                                                                                                  tgATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAAc                                                                                  >  7:73362/1‑93 (MQ=255)
                                                                                                                    aTAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAAccc                        >  4:1737/1‑149 (MQ=255)
                                                                                                                     tAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGa                                    <  4:248684/136‑1 (MQ=255)
                                                                                                                     tAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGa                                    >  3:248684/1‑136 (MQ=255)
                                                                                                                       aaGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCATAGCAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCg                       >  1:125451/1‑149 (MQ=255)
                                                                                                                       aaGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGACACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGGCAACCCGGt                     >  8:165162/1‑149 (MQ=255)
                                                                                                                       aaGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGt                     <  7:87952/149‑1 (MQ=255)
                                                                                                                         gTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGa                   >  4:357433/1‑149 (MQ=255)
                                                                                                                            gctcgctcGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATg                 <  7:53960/148‑1 (MQ=255)
                                                                                                                             ctcgctcgAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATg                 >  7:345378/1‑147 (MQ=255)
                                                                                                                             ctcgctcgAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATg                 <  8:345378/147‑1 (MQ=255)
                                                                                                                             ctcgctcgAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCt               >  5:383292/1‑149 (MQ=255)
                                                                                                                                gctcgAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGt             <  5:238949/148‑1 (MQ=255)
                                                                                                                                    atggacGGCCTGCGGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGa                   <  1:353068/134‑1 (MQ=255)
                                                                                                                                    atggacGGCCTGCGGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGa                   >  2:353068/5‑138 (MQ=255)
                                                                                                                                       cacGGCCCGCAGGCCGTCCATAATATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGTACTCGATc     <  4:12092/149‑1 (MQ=255)
                                                                                                                                       cacGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGTACTCGATc     <  6:287599/149‑1 (MQ=255)
                                                                                                                                         cGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGTACTCGATCgg   <  6:111622/149‑1 (MQ=255)
                                                                                                                                           gCCCGCAGGCCGTCCATAAGATCGTTTTTCCTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGTACTCGATCGGc  >  1:34389/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCCTCATAGGAGCGGGTTCACCCGCGAAGAGGCCGGCCCTGCCAGCCCAAACCCCCAGTCATATCGTTAGACTTTTTCGATCTAGCCAGCCTCCTAAATAGAAGCAGCTGCTGTGATAAAGTGGCTCGCTCGAACACGGCCCGCAGGCCGTCCATAAGATCGTTTTTCGTGAGTAGCCATGGTCGACAAACTGACGCACTTGAAACAGCTGGAGGCGGAGAGCATCCACATCATCCGCGAGGTGGCCGCCGAGTTCGACAACCCGGTGATGCTGTACTCGATCGGC  >  NC_002947/1489226‑1489511

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: