Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,983,831 G→C 6.1% C282S (TGC→TCC) ‡ ilvE → branched‑chain‑amino acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,983,8310GC 6.1% 209.1 / 3.7 65C282S (TGC→TCC) ‡ilvEbranched‑chain‑amino acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (29/32);  new base C (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACC  >  NC_002947/3983683‑3983978
                                                                                                                                                    |                                                                                                                                                   
gAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTg                                                                                                                                                     >  8:246238/1‑149 (MQ=255)
           cccGAAATCGGCCTCGGTACTACCAGGCATCAGCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCgg                                                                                                                                          >  5:240336/1‑149 (MQ=255)
             ggAAATCGGCATCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTTGCGCAATCGCGCCTGGGCATGACCTTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGcc                                                                                                                                        <  8:83523/148‑1 (MQ=255)
             cGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCACTCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGcc                                                                                                                                        <  1:321175/149‑1 (MQ=255)
             cGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGcc                                                                                                                                        >  5:242917/1‑149 (MQ=255)
             cGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGcc                                                                                                                                        <  4:406582/149‑1 (MQ=255)
             cGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCGCAGGCACcggcttc                                                                                                                                        <  3:12700/149‑6 (MQ=255)
                     gCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCAc                                                                                                                                <  5:174945/149‑1 (MQ=255)
                       cTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGCCCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCAccc                                                                                                                              >  7:268504/1‑149 (MQ=255)
                               tGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATcggc                                                                                                                      >  7:103660/1‑149 (MQ=255)
                               tGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCGCAGGCACcggcttc                                                                                                                                        >  5:185736/1‑126 (MQ=255)
                               tGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCGCAGGCACcggcttc                                                                                                                                        <  6:185736/131‑6 (MQ=255)
                                 ccAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATcggcgg                                                                                                                    >  5:134467/1‑149 (MQ=255)
                                 ccAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATcggcgg                                                                                                                    >  1:329356/1‑149 (MQ=255)
                                     gCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATc                                                                                                                <  6:402261/149‑1 (MQ=255)
                                      cATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATc                                                                                                                >  8:189072/1‑148 (MQ=255)
                                       aTCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGa                                                                                                              <  7:246238/149‑1 (MQ=255)
                                         cACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAg                                                                                                             >  4:246157/1‑148 (MQ=255)
                                            ccGCCTCTCGCTGATGGAACTGGCGCGATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTAc                                                                                                          >  1:232974/1‑148 (MQ=255)
                                               ccTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGGGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACg                                                                                                      <  6:240336/149‑1 (MQ=255)
                                                   tcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGc                                                                                                    >  6:85658/1‑147 (MQ=255)
                                                       tGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCaa                                                                                                  <  5:210276/145‑1 (MQ=255)
                                                       tGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCaa                                                                                                  >  6:210276/1‑145 (MQ=255)
                                                              aCTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACgtgt                                                                                       >  2:322118/1‑149 (MQ=255)
                                                                  gcgcAATCGGGCCTGGGCATGACGGTGATGGAAGGCGATGTCGAGTTCAGCATGCTGGACGGCTTCGTGGAAGCCGGTGGCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTcg                                                                                    <  2:412491/148‑2 (MQ=255)
                                                                  gcgcAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCa                                                                                   >  8:334983/1‑149 (MQ=255)
                                                                   cgcAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATc                                                                                                                <  2:129891/119‑1 (MQ=255)
                                                                   cgcAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATc                                                                                                                >  1:129891/1‑119 (MQ=255)
                                                                   cgcAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACAGCGGCCGTTATCACCCCGATCAGCGGCATCGAGTACAACGGCAAGCTGCACCTGTTacac                                                                                  >  3:258445/1‑149 (MQ=255)
                                                                    gcAATCGCGCCTGCGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACAGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCAc                                                                                  >  2:85786/1‑148 (MQ=255)
                                                                     cAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGa                                                                                <  4:282809/149‑1 (MQ=255)
                                                                     cAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGa                                                                                <  7:39690/149‑1 (MQ=255)
                                                                      aaTCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGAc                                                                               <  8:103660/149‑1 (MQ=255)
                                                                       aTCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGAcc                                                                              >  5:51650/1‑149 (MQ=255)
                                                                        tCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCt                                                                             <  6:134467/149‑1 (MQ=255)
                                                                        tCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCt                                                                             <  2:329356/149‑1 (MQ=255)
                                                                        tCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCt                                                                             >  8:352973/1‑149 (MQ=255)
                                                                          gcgcCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTg                                                                            <  8:111962/148‑1 (MQ=255)
                                                                                  gCATGACCGTGATCGAAGGGGATGGCGAGATCAGCAAGGTGGACGGCTTCGTCGAAGCCGGTGCCTGCGGCGCCGCGGCGGTTGTGACGCCGAGCGGGGGCATCGAGTACAACGGGCAGCGGGAGGTGTTCCACGAGCTGGATAAGGTc                                                                   >  4:237966/1‑149 (MQ=255)
                                                                                   cATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCGGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCg                                                                  <  2:212147/149‑1 (MQ=255)
                                                                                   cATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCg                                                                  <  1:85786/149‑1 (MQ=255)
                                                                                   cATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCg                                                                  <  1:322118/149‑1 (MQ=255)
                                                                                           tGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGt                                                                    >  4:420418/1‑139 (MQ=255)
                                                                                           tGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGt                                                                    <  3:420418/139‑1 (MQ=255)
                                                                                                  aGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGa                                                   <  2:232974/149‑1 (MQ=255)
                                                                                                   ggCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGaa                                                  <  4:13119/149‑1 (MQ=255)
                                                                                                       aTGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGctc                                              >  6:168132/1‑149 (MQ=255)
                                                                                                        tGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGctct                                             <  7:189072/149‑1 (MQ=255)
                                                                                                                  aGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGc                                    <  6:242917/148‑1 (MQ=255)
                                                                                                                   gCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCt                                   >  2:315415/1‑148 (MQ=255)
                                                                                                                    cAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGa                                 >  1:293057/1‑149 (MQ=255)
                                                                                                                      aGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGAcc                               >  8:162807/1‑149 (MQ=255)
                                                                                                                             ccGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATcc                        <  1:315415/149‑1 (MQ=255)
                                                                                                                              cGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAAc                                        >  5:337489/1‑132 (MQ=255)
                                                                                                                              cGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAAc                                        <  6:337489/132‑1 (MQ=255)
                                                                                                                              cGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCa                       <  3:237966/149‑1 (MQ=255)
                                                                                                                               gCTTCGTCGAAGCCGTTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAATCTGCACGTGTACCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCag                      <  8:268504/149‑1 (MQ=255)
                                                                                                                                 ttCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCagag                    >  5:189928/1‑149 (MQ=255)
                                                                                                                                   cgtcgAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCagag                    <  6:189928/147‑1 (MQ=255)
                                                                                                                                     tcgAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAgcggc                <  4:131205/149‑1 (MQ=255)
                                                                                                                                         ccgcGGTGCCGCAGGCACCGCGGCCGTTATCACCCCGATCGGCTGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATg            >  3:276351/4‑149 (MQ=255)
                                                                                                                                             ggTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTaga        <  7:162807/149‑1 (MQ=255)
                                                                                                                                                  cTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCAc   <  7:334983/149‑1 (MQ=255)
                                                                                                                                                   tGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCAcc  >  2:178119/1‑149 (MQ=255)
                                                                                                                                                   tGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCAcc  >  1:6971/1‑149 (MQ=255)
                                                                                                                                                   tGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCAcc  >  1:14238/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACC  >  NC_002947/3983683‑3983978

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