Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,845,977 C→T 9.4% intergenic (+49/+42) 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,9770CT 9.4% 174.5 / 6.7 65intergenic (+49/+42)yfhL/mutM4Fe‑4S cluster‑containing protein/formamidopyrimidine‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base C (32/26);  new base T (3/3);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-01

GCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACA  >  NC_002947/5845829‑5846122
                                                                                                                                                      |                                                                                                                                                 
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                                                                                                        <  3:91525/103‑1 (MQ=255)
                                                                                           gCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTc                                                                                                        >  4:91525/1‑103 (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                          <  6:235886/146‑1 (MQ=255)
                                                                                                                            aacTGTGGGAGCAGGCATGCCCGCTCTCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGc                          >  5:235886/3‑148 (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)
                                                                                                                                                    tCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTaca  >  1:247040/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GCAGGTCTGCCCGGTCGACTGCATCCCGCTGGATGAAGCGCATCCGGAAACCGAAGAACAGCTGATGGCCAAGTATCGCCGGAT‑‑CACTGGCAAGGGCTGATCATTGCCCTTCCCTGCACGGTCCCTGTGAGAGCGGGCATGCCTGCTCCCACAGTTTAGGTTTCGCTTCAACAGGTTCGCAGCGAACGTTTCATCGTTGGCAGCGCGGGCAGTACACGCTCGCCCTCTGCCCCAGCTTCACTTCGCGCAACTCAGTGCCGCATACCTTGCACGGCTGCCCGCCACGCCCGTACA  >  NC_002947/5845829‑5846122

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: