Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,761,204 A→T 7.9% K86* (AAG→TAG) ‡ gpmI → 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,761,2040AT 7.9% 182.4 / 5.4 63K86* (AAG→TAG) ‡gpmI2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base T (3/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CCTGTGCGCCACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGG  >  NC_002947/5761059‑5761351
                                                                                                                                                 |                                                                                                                                                   
ccTGTGCGCCACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAgg                                                                                                                                                  <  6:118574/149‑1 (MQ=255)
        gccaccAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCcgc                                                                                                                                          >  2:197138/5‑149 (MQ=255)
        ccACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCcgc                                                                                                                                          >  7:133957/1‑149 (MQ=255)
          aCCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGa                                                                                                                                        >  7:276480/1‑149 (MQ=255)
                ccGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGc                                                                                                                                    >  4:117812/1‑147 (MQ=255)
                       gCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttctt                                                                                                                            <  7:134353/148‑1 (MQ=255)
                                 cGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCgg                                                                                                                 <  7:3872/149‑1 (MQ=255)
                                 cGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCgg                                                                                                                 <  6:98095/149‑1 (MQ=255)
                                    cTCGGGCATGGATGTCGGCCTGCCGGACGGGAAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGGCGTGTATCAGTACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTAc                                                                                                              >  5:289016/1‑149 (MQ=255)
                                       gggCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACt                                                                                                             <  3:117812/147‑1 (MQ=255)
                                       gggCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGa                                                                                                           >  6:181136/1‑149 (MQ=255)
                                        ggCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAc                                                                                                          >  7:153742/1‑149 (MQ=255)
                                            tGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGt                                                                                                      <  5:181136/149‑1 (MQ=255)
                                               aTGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAc                                                                                                          <  7:243493/142‑1 (MQ=255)
                                               aTGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAc                                                                                                          >  8:243493/1‑142 (MQ=255)
                                                tGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCtg                                                                                                  >  3:224943/1‑149 (MQ=255)
                                                  tCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCtgtg                                                                                                <  8:276480/149‑1 (MQ=255)
                                                   ccatcTGCCCGTCCGGCAGATGGGCAACTCCGAAGTCGGTCATATGAGCCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGGGACGGCGAGTTCTTCGAAAACCCGGTACTGAGCGGTGCTGTgg                                                                                               >  4:210964/5‑149 (MQ=255)
                                                    ggCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTgg                                                                                               <  2:171285/148‑1 (MQ=255)
                                                      ccTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGAc                                                                                             >  1:155803/1‑148 (MQ=255)
                                                             gACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAATACCCGGTACTGACCGGTGCTGTGGACAAGGCAGt                                                                                     <  1:197138/149‑1 (MQ=255)
                                                              aCGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCtgtg                                                                                                <  6:18127/137‑1 (MQ=255)
                                                              aCGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCtgtg                                                                                                >  5:18127/1‑137 (MQ=255)
                                                               cGGGCATATGGGCAACTCCGAAGTCGGTCATATGAAACTAGGCGCAGGTCGCGTCGTGTATCAGGACTTCACTCGGGTGACCAAGGCCATCCGCGACGGCTAGTTCTTCGAAAACCCGGTACTGACCTGTGCTGTGTAAAAGGCAGTCg                                                                                   <  6:289016/149‑1 (MQ=255)
                                                                   cAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCTGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGc                                                                               <  5:114152/149‑1 (MQ=255)
                                                                    aGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCt                                                                              >  4:260544/1‑149 (MQ=255)
                                                                     gATGGGCACCGCCGAAGTCGGTCATATGCACCTCGGCGCAGGTCGCGTCGTGTATCAGGCCTTCACCCGGGTGACCAAGCCCATCCGCGACGGCGACTTCTTCGAAAACCCGGTACTGACCGGTGCTGCGGACAAGGCAGTCGGTGCTg                                                                             <  8:133957/149‑1 (MQ=255)
                                                                     cctGGGCCACTCCTCAGTCGGCCATAGTAACCTCGGCGCAGGTCGCGTCGTCTATCTGGACTTCACGCGTGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTg                                                                             <  4:224943/147‑1 (MQ=255)
                                                                           cAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAgg                                                                       >  7:275833/1‑149 (MQ=255)
                                                                                 cGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTg                                                                  >  1:144970/1‑148 (MQ=255)
                                                                                 cGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGc                                                                 <  8:275833/149‑1 (MQ=255)
                                                                                            atatGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGa                                                      <  2:155803/149‑1 (MQ=255)
                                                                                            atatGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGa                                                      >  6:66643/1‑149 (MQ=255)
                                                                                              atcaaCCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAACGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGAGCACATCCTCGGAct                                                    >  5:126999/4‑149 (MQ=255)
                                                                                               tGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGATAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGActg                                                   >  4:15254/1‑149 (MQ=255)
                                                                                                gAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTg                                                                  <  5:6228/133‑1 (MQ=255)
                                                                                                gAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTg                                                                  >  6:6228/1‑133 (MQ=255)
                                                                                                   ccTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGActgctg                                                >  4:94507/1‑148 (MQ=255)
                                                                                                    cTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTc                                              >  1:229675/1‑149 (MQ=255)
                                                                                                    cTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTc                                              >  5:755/1‑149 (MQ=255)
                                                                                                       ggCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACAGCAGTCCGa                                           <  6:755/149‑1 (MQ=255)
                                                                                                            aGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                      >  5:182860/1‑149 (MQ=255)
                                                                                                            aGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                      <  5:66643/149‑1 (MQ=255)
                                                                                                              gTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGATAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATggtggt                                    <  3:15254/149‑1 (MQ=255)
                                                                                                                       tgtATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTTCTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGc                           >  5:239530/1‑149 (MQ=255)
                                                                                                                           tCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccac                        >  8:195312/1‑148 (MQ=255)
                                                                                                                                 cTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGACAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGAc                  >  1:94533/1‑148 (MQ=255)
                                                                                                                                  ttCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCa                >  3:294738/1‑149 (MQ=255)
                                                                                                                                    gtcgcgGATGGCCTTGGTCACCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTAcaca                             >  5:104243/3‑134 (MQ=255)
                                                                                                                                    gtcgcgGATGGCCTTGGTCACCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTAcaca                             <  6:104243/132‑1 (MQ=255)
                                                                                                                                    gtcgcgGATGGCCTTGGTCACCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACCCAGCCCCCAGGACCCCt              >  7:285939/3‑149 (MQ=255)
                                                                                                                                     tcgcgGATGGCCTTGGTCACCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                      >  3:158697/2‑124 (MQ=255)
                                                                                                                                     tcgcgGATGGCCTTGGTCACCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                      <  4:158697/123‑1 (MQ=255)
                                                                                                                                     aCGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTacac                              <  7:212451/132‑1 (MQ=255)
                                                                                                                                     aCGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTacac                              >  8:212451/1‑132 (MQ=255)
                                                                                                                                       gcgGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTac                                >  4:258794/1‑128 (MQ=255)
                                                                                                                                       gcgGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTac                                <  3:258794/128‑1 (MQ=255)
                                                                                                                                        cgGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtggt          <  3:260544/149‑1 (MQ=255)
                                                                                                                                          ggTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtggtgg        <  2:257112/149‑1 (MQ=255)
                                                                                                                                          ggTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtggtgg        <  8:153742/149‑1 (MQ=255)
                                                                                                                                              aCCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACATAGCCACCAGGACCACTTGGTGGCCAt    <  3:210964/149‑1 (MQ=255)
                                                                                                                                               ccAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATg   >  6:251678/1‑149 (MQ=255)
                                                                                                                                                cAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGTAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATgg  >  5:23241/1‑149 (MQ=255)
                                                                                                                                                 aaGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTgg                                                                            <  3:98585/74‑1 (MQ=255)
                                                                                                                                                 aaGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTgg                                                                            >  4:98585/1‑74 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCTGTGCGCCACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTCACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGGCCATGG  >  NC_002947/5761059‑5761351

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: ‑