Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,228,274 A→C 7.9% T143P (ACG→CCG) ‡ holB → DNA polymerase III subunit delta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,228,2740AC 7.9% 184.0 / 5.5 63T143P (ACG→CCG) ‡holBDNA polymerase III subunit delta'
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base C (2/3);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GGTGGCCTTTGTGGTGCAGACCGCGCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGCGACGAGGTCGAGCGCGATGAGC  >  NC_002947/2228129‑2228412
                                                                                                                                                 |                                                                                                                                          
ggtggCCTTTGTGGTGCAGACCGCGCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACg                                                                                                                                          >  8:319801/1‑148 (MQ=255)
  tggCCTTTGTGGTGCAGACCGCGCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCAc                                                                                                                                       <  5:153582/149‑1 (MQ=255)
     ccTTTGTGGTGCAGACCGCGCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCAccag                                                                                                                                    <  8:156463/149‑1 (MQ=255)
      cTTTGTGGTGCAGACCGCGCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCAccagc                                                                                                                                   >  1:230815/1‑149 (MQ=255)
       tttGTGGTGCAGACCGCGCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCAccagca                                                                                                                                  >  7:261295/1‑148 (MQ=255)
         tgtgGTGCAGACCGCGCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGCGTCAccagcaa                                                                                                                                 <  5:256265/148‑3 (MQ=255)
         tgtgGTGCAGACCGCGCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCAccagcc                                                                                                                                  <  1:36951/147‑1 (MQ=255)
            ggtgCAGACCGCGCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACcagcccagcc                                                                                                                             >  1:262681/1‑149 (MQ=255)
              tgCAGACCGCGCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGGGACCAACGCCCTGCTCTAAAGCATTGAAGAACCTTCCGGTGATACAGTGTTGTTGCTGGTGACGCACAAGCCCAGCCg                                                                                                                            >  1:304089/1‑148 (MQ=255)
                   aCCGCGCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGtt                                                                                                                      <  2:67755/149‑1 (MQ=255)
                     cgcgCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGCGTCAccagcaac                                                                                                                                >  2:322717/1‑134 (MQ=255)
                     cgcgCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGCGTCAccagcaac                                                                                                                                <  1:322717/137‑4 (MQ=255)
                           gTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCcgacga                                                                                                              <  2:230815/149‑1 (MQ=255)
                            ttAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCcgacga                                                                                                              >  4:234082/1‑148 (MQ=255)
                               ggcgggcgCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCaa                                                                                                          <  2:254781/149‑1 (MQ=255)
                                    gcgcAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCACTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGcc                                                                                                     <  7:136060/149‑1 (MQ=255)
                                       cAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCt                                                                                                  >  3:121971/1‑149 (MQ=255)
                                       cAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCt                                                                                                  >  5:42025/1‑149 (MQ=255)
                                          ggtggtGTTGATCGAGCCGGTGGAGGCGTTGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGc                                                                                                >  1:209378/1‑148 (MQ=255)
                                             ggtgTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGACCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATGCCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGcacgc                                                                                            >  7:352497/1‑145 (MQ=255)
                                                gtTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGCTCAAGAGCCGCTGCCAGCAgg                                                                                         >  6:78828/1‑149 (MQ=255)
                                                 ttGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGt                                                                                        >  4:39434/1‑149 (MQ=255)
                                                    aTCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGc                                                                                     >  2:132045/1‑149 (MQ=255)
                                                     tCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGcc                                                                                    >  6:306935/1‑149 (MQ=255)
                                                         gCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTc                                                                                <  3:39434/149‑1 (MQ=255)
                                                            gttGGAGGCGTTGACATTCATCGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCACTGCCAGCAGGTGGCCTGTCCg                                                                              >  8:259009/3‑148 (MQ=255)
                                                             ttGGAGGCGTTGACATTCATCGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCg                                                                              >  6:127084/2‑147 (MQ=255)
                                                                gAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGcc                                                                         >  8:97866/1‑149 (MQ=255)
                                                                   gCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAg                                                                      >  2:59518/1‑149 (MQ=255)
                                                                     gATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTc                                                                    >  7:75454/1‑149 (MQ=255)
                                                                      aTGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCt                                                                   <  3:234082/149‑1 (MQ=255)
                                                                         aaTGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGc                                                                <  7:97866/149‑1 (MQ=255)
                                                                          aTGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGc                                                                <  4:121971/148‑1 (MQ=255)
                                                                               aaCGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCagag                                                          >  8:334421/1‑149 (MQ=255)
                                                                               aaCGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCagag                                                          >  6:337021/1‑149 (MQ=255)
                                                                               aaCGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCaga                                                           <  2:262681/148‑1 (MQ=255)
                                                                                aCGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGc                                                         <  5:323689/149‑1 (MQ=255)
                                                                                aCGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGc                                                         >  3:244254/1‑149 (MQ=255)
                                                                                     tCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCt                                                                                                  >  2:90780/1‑103 (MQ=255)
                                                                                     tCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCt                                                                                                  <  1:90780/103‑1 (MQ=255)
                                                                                        aaCGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCAtg                                                 <  7:319801/149‑1 (MQ=255)
                                                                                                 cTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCt                                                                                                  <  3:134574/91‑1 (MQ=255)
                                                                                                 cTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCt                                                                                                  >  4:134574/1‑91 (MQ=255)
                                                                                                    aaaaGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCAtggctggctggc                                       <  6:160897/147‑1 (MQ=255)
                                                                                                    aaaaGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCAtggctggctggc                                       <  5:306935/147‑1 (MQ=255)
                                                                                                    aaaaGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGgcgc                                     <  5:127084/149‑1 (MQ=255)
                                                                                                    aaaaGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGCCAtgcctggctctgg                                      <  5:78828/148‑5 (MQ=255)
                                                                                                          cTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCAAGCGCGTTAcc                               >  1:316777/1‑149 (MQ=255)
                                                                                                             gaagaaCCTTCCGGTGATCCCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTAccc                              <  8:261295/147‑1 (MQ=255)
                                                                                                             gaagaaCCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTAccc                              <  1:59518/147‑1 (MQ=255)
                                                                                                                  aCCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGc                        <  8:75454/148‑1 (MQ=255)
                                                                                                                   ccTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGc                                                                                                   <  3:19941/72‑1 (MQ=255)
                                                                                                                   ccTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGc                                                                                                   >  4:19941/1‑72 (MQ=255)
                                                                                                                             aTACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGCGACGAGGTCGa             <  2:316777/148‑1 (MQ=255)
                                                                                                                              tACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGCGACGAGGTGGAgc           >  7:112916/1‑149 (MQ=255)
                                                                                                                              tACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGCGACGAGGTCGAgc           <  7:334421/149‑1 (MQ=255)
                                                                                                                                 cGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGCGACGAGGTCGAgcgc         >  2:125580/1‑148 (MQ=255)
                                                                                                                                 cGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGCGACGAGGTCGAgcgc         >  1:24341/1‑148 (MQ=255)
                                                                                                                                  gtgtTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGCGACGAGGTCGAGCGCGa       <  1:194831/149‑1 (MQ=255)
                                                                                                                                       gttgCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGCGACGAGGTCGAGCGCGATGAgc  >  1:76322/1‑149 (MQ=255)
                                                                                                                                       gttgCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGCGACGAGGTCGAGCGCGATGAgc  >  2:253676/1‑149 (MQ=255)
                                                                                                                                         ggCTGGTGCGTCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGCGACGAg                  <  4:158546/130‑1 (MQ=255)
                                                                                                                                         ggCTGGTGCGTCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGCGACGAg                  >  3:158546/2‑131 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
GGTGGCCTTTGTGGTGCAGACCGCGCAGTTAGGCGGGCGCAAGGTGGTGTTGATCGAGCCGGTGGAGGCGATGAATGTCAACGCCTCCAACGCCCTGCTCAAAAGCCTTGAAGAACCTTCCGGTGATACCGTGTTGTTGCTGGTGACGCACCAGCCCAGCCGTCTGTTGCCGACGATCAAGAGCCGCTGCCAGCAGGTGGCCTGTCCGCAGCCCAGTCTGGCCCAGAGCCAGGCATGGCTGGCTGGCGCGTTACCCGACTGCGACGAGGTCGAGCGCGATGAGC  >  NC_002947/2228129‑2228412

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