Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,845,963 G→A 8.4% intergenic (+35/+56) yfhL → / ← mutM 4Fe‑4S cluster‑containing protein/formamidopyrimidine‑DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,845,9630GA 8.4% 198.9 / 5.9 71intergenic (+35/+56)yfhL/mutM4Fe‑4S cluster‑containing protein/formamidopyrimidine‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (32/33);  new base A (3/3);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

AGCCGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCC  >  NC_002947/5845816‑5846107
                                                                                                                                                     |                                                                                                                                                
aGCCGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCTAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCgg                                                                                                                                                 <  8:274039/149‑1 (MQ=255)
  ccgccgTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCGGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGGCCCTGTGAGAGCGGGc                                                                                                                                               <  1:110737/149‑1 (MQ=255)
  ccGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCTCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGc                                                                                                                                               <  1:115348/149‑1 (MQ=255)
  ccGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCggg                                                                                                                                                >  3:79987/1‑148 (MQ=255)
    gCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCAt                                                                                                                                             <  6:180350/149‑1 (MQ=255)
    gCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCAt                                                                                                                                             <  4:79987/149‑1 (MQ=255)
      aGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGCAAATCGATTAACTGCTTATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGc                                                                                                                                           <  1:252776/149‑1 (MQ=255)
             gcagGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTcc                                                                                                                                    <  2:280984/149‑1 (MQ=255)
               agGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCca                                                                                                                                  <  1:158744/149‑1 (MQ=255)
                 gTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGGGAGCAGGCATGCCCGCTCTca                                                                                                                                  >  5:196230/1‑147 (MQ=255)
                 gTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGCGCGGGCATGCCTGCTcccc                                                                                                                                  >  5:10304/1‑146 (MQ=255)
                 gTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCca                                                                                                                                  >  2:298168/1‑147 (MQ=255)
                    tGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATAATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAggt                                                                                                                             >  6:7699/1‑147 (MQ=255)
                      cccGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGGGAGCAGGCATGCCCGCTCTCACAgggac                                                                                                                           <  6:196230/149‑5 (MQ=255)
                      cccGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTa                                                                                                                           <  4:287237/149‑1 (MQ=255)
                             gACTGCATCCCGCTGGATGAAGCGCATCCGTAAACCGAAGAACAGCTGAGGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGttt                                                                                                                            <  8:195789/141‑1 (MQ=255)
                             gACTGAATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCAGGAT‑‑CACTGGCAAGCGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGttt                                                                                                                            >  7:195789/1‑141 (MQ=255)
                              aCTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGATTGTGCTGGAAAGGGCTGATCATTGCCCTACCCTGAACGGTCCCTGTGAGAGCGGGAATCCCTGATCCCACAGTTTAGGTTTc                                                                                                                     >  1:249697/1‑149 (MQ=255)
                                   aTCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCaa                                                                                                              <  2:89420/149‑1 (MQ=255)
                                         ctcgATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGtt                                                                                                        <  5:285405/146‑1 (MQ=255)
                                             aTGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCa                                                                                                    <  8:291856/149‑1 (MQ=255)
                                              tGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAg                                                                                                   >  2:35855/1‑149 (MQ=255)
                                               gAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGc                                                                                                  >  2:162547/1‑149 (MQ=255)
                                                  gcgcATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGaa                                                                                               <  2:249697/149‑1 (MQ=255)
                                                          ggAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCaa                                                                                                              >  6:258217/1‑126 (MQ=255)
                                                          ggAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCaa                                                                                                              <  5:258217/126‑1 (MQ=255)
                                                          ggAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCAt                                                                                       >  5:49963/1‑149 (MQ=255)
                                                                cGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                 <  6:49963/149‑1 (MQ=255)
                                                                cGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                 >  6:131073/1‑149 (MQ=255)
                                                                      aCAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAgcgc                                                                           <  1:162547/149‑1 (MQ=255)
                                                                       cAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGcgcg                                                                          >  7:3018/1‑149 (MQ=255)
                                                                         gCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGAAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCggg                                                                        <  8:4983/149‑1 (MQ=255)
                                                                             aTGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAAc                                                                                                             <  4:260470/108‑1 (MQ=255)
                                                                             aTGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAAc                                                                                                             >  3:260470/1‑108 (MQ=255)
                                                                             aTGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGt                                                                    <  8:3018/149‑1 (MQ=255)
                                                                                 ccAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTc                                                                                                       <  6:311614/110‑1 (MQ=255)
                                                                                 ccAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTc                                                                                                       >  5:311614/1‑110 (MQ=255)
                                                                                   aaGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGc                                                              >  5:177584/1‑149 (MQ=255)
                                                                                      tATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGATCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCg                                                           >  7:150214/1‑149 (MQ=255)
                                                                                          gCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCct                                                       >  1:47431/1‑149 (MQ=255)
                                                                                                t‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGccc                                                 >  4:46634/1‑149 (MQ=255)
                                                                                                     cTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAg                                              <  7:128872/149‑1 (MQ=255)
                                                                                                     cTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGGTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAg                                              >  5:251104/1‑149 (MQ=255)
                                                                                                        gCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTc                                                                                         >  4:91525/1‑103 (MQ=255)
                                                                                                        gCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTc                                                                                         <  3:91525/103‑1 (MQ=255)
                                                                                                        gCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCtt                                           >  6:148908/1‑149 (MQ=255)
                                                                                                                tGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCACCCTCTGCCCCAGCTTCACTTcgc                                   >  1:43220/1‑149 (MQ=255)
                                                                                                                  aTCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                 <  5:27378/101‑1 (MQ=255)
                                                                                                                  aTCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTgg                                                                                 >  6:27378/1‑101 (MQ=255)
                                                                                                                     aTTGCGCTTCCCTGCACGGTCCCTGTGAGAGCGGGGATGCCTGGTCGCACAGTTTAGGTTTCGCTTCTACAGGTTCGCAGCGAACGTTTGATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGAAAc                              <  6:251104/149‑1 (MQ=255)
                                                                                                                      ttGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACt                             >  8:168105/1‑149 (MQ=255)
                                                                                                                         cccTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCa                           <  6:177584/148‑1 (MQ=255)
                                                                                                                         cccTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGCTTCGCTTCAACAGGTTCGCAGCTAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCa                           <  5:131073/148‑1 (MQ=255)
                                                                                                                          ccTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGt                         <  3:46634/149‑1 (MQ=255)
                                                                                                                            ttCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcgc                                   <  3:104103/137‑1 (MQ=255)
                                                                                                                            ttCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTcgc                                   >  4:104103/1‑137 (MQ=255)
                                                                                                                            ttCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGc                       >  8:15935/1‑149 (MQ=255)
                                                                                                                              cccTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCg                     >  4:48815/1‑149 (MQ=255)
                                                                                                                               ccTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCg                     <  1:298168/148‑1 (MQ=255)
                                                                                                                                cTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCACCTCAGTGCCGCa                   >  6:230607/1‑149 (MQ=255)
                                                                                                                                 tGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCAc                                        <  7:131073/127‑1 (MQ=255)
                                                                                                                                 tGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCAc                                        >  8:131073/1‑127 (MQ=255)
                                                                                                                                    aCGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATAcc               <  3:48815/149‑1 (MQ=255)
                                                                                                                                         aacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGc           >  5:235886/3‑148 (MQ=255)
                                                                                                                                         aacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGc           <  6:235886/146‑1 (MQ=255)
                                                                                                                                          ccTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCAc         >  5:247611/1‑149 (MQ=255)
                                                                                                                                          acTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCa                   >  2:167522/2‑139 (MQ=255)
                                                                                                                                          acTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCa                   <  1:167522/138‑1 (MQ=255)
                                                                                                                                           cTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCatg        >  8:64392/1‑147 (MQ=255)
                                                                                                                                            tgtgAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACgg       >  7:21776/1‑149 (MQ=255)
                                                                                                                                                agagCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGc   <  6:10304/149‑1 (MQ=255)
                                                                                                                                                  agCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGcc  <  4:90133/148‑1 (MQ=255)
                                                                                                                                                  agCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGcc  >  3:90133/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AGCCGCAGTGCCAGCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCC  >  NC_002947/5845816‑5846107

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: