Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,983,800 C→T 6.1% S257L (TCG→TTG)  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,8000CT 6.1% 192.8 / 2.8 65S257L (TCG→TTG) ilvEbranched‑chain‑amino acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (29/32);  new base T (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 = 7.83e-01

ACCTCGACCCGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGT  >  NC_002947/3983653‑3983937
                                                                                                                                                   |                                                                                                                                         
aCCTCGACCCGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATcgc                                                                                                                                         >  2:1071675/1‑150 (MQ=255)
  cTCGACCCGCTGACCCACACCAAGATCGAGGAAGTCGGCTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCgcgc                                                                                                                                       >  1:731403/1‑150 (MQ=255)
   tCGACCCGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCAGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGcc                                                                                                                                      >  1:726094/1‑150 (MQ=255)
     gACCCGCTGACCCACACCAAGCTCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTg                                                                                                                                    <  1:573029/150‑1 (MQ=255)
       cccGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTggg                                                                                                                                  <  1:488122/150‑1 (MQ=255)
        ccGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCAGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGc                                                                                                                                 <  2:726094/150‑1 (MQ=255)
             gACCCACACCAAGATCGAGGAAGCTGTCTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGAc                                                                                                                            <  1:36/150‑1 (MQ=255)
                ccACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGt                                                                                                                         >  2:780614/1‑150 (MQ=255)
                  acacCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGa                                                                                                                       <  1:780614/150‑1 (MQ=255)
                      cAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGa                                                                                                                   >  1:498325/1‑150 (MQ=255)
                      cAAGATCGAGGAAGTCGGTTCGGCCAACTTATTCGGCATCACCGCAATCGCAGGATTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAAGCGCGCCGGGGCATGACCGTTATCGa                                                                                                                   >  2:470681/1‑150 (MQ=255)
                        aGATCGAGGAAGTCGGCTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAg                                                                                                                 <  2:731403/150‑1 (MQ=255)
                         gATCGAGGAAGTCGGTTCTCCCAACTTCTTCGTCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAgg                                                                                                                <  2:568163/150‑1 (MQ=255)
                             gAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGAt                                                                                                            >  1:436175/1‑150 (MQ=255)
                               ggAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGGCCGTGATCGAAGGCGATGt                                                                                                          >  2:427824/1‑150 (MQ=255)
                                       gTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCa                                                                                                  >  2:971937/1‑150 (MQ=255)
                                       gTTCGACCAAACTCCTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATCGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCa                                                                                                  <  2:335820/150‑1 (MQ=255)
                                           ggCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATATCGAGATCAGCaa                                                                                              <  2:812866/150‑1 (MQ=255)
                                               aaCTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTg                                                                                          >  1:224341/1‑150 (MQ=255)
                                                aCTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCGATTGCGccagtt                                                                                                                                  >  1:190195/1‑105 (MQ=255)
                                                aCTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCGATTGCGccagtt                                                                                                                                  <  2:190195/109‑5 (MQ=255)
                                                 cttcttCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGa                                                                                        <  1:776219/150‑1 (MQ=255)
                                                  ttcttcGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGAc                                                                                       <  2:909730/150‑1 (MQ=255)
                                                              aCCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCGATTGCGcca                                                                                                                                     >  1:815611/1‑91 (MQ=255)
                                                              aCCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCGATTGCGcca                                                                                                                                     <  2:815611/92‑2 (MQ=255)
                                                               ccGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAg                                                                          >  2:487475/1‑150 (MQ=255)
                                                                cGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGc                                                                         >  2:898276/1‑150 (MQ=255)
                                                                 gCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGcc                                                                        >  1:1003941/1‑150 (MQ=255)
                                                                    aacaacGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGt                                                                     >  1:1146305/1‑150 (MQ=255)
                                                                      caacGAGTTCGTAACCCCGCAACCGGCCTAGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGc                                                                   <  2:498325/150‑1 (MQ=255)
                                                                       aacGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGa                                                                            <  1:810975/140‑1 (MQ=255)
                                                                       aacGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGa                                                                            >  2:810975/1‑140 (MQ=255)
                                                                         cGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCgg                                                                      >  2:539503/1‑144 (MQ=255)
                                                                         cGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCgg                                                                      <  1:539503/144‑1 (MQ=255)
                                                                            gTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCgg                                                             <  2:224341/150‑1 (MQ=255)
                                                                                     cccGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGc                                                                                                                                 <  2:817244/73‑1 (MQ=255)
                                                                                     cccGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGc                                                                                                                                 >  1:817244/1‑73 (MQ=255)
                                                                                         aaaTCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGtt                                                <  1:898276/150‑1 (MQ=255)
                                                                                              ggCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCAc                                           <  1:1071675/150‑1 (MQ=255)
                                                                                              ggCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCCGCCGCTATCAc                                           >  2:305678/1‑150 (MQ=255)
                                                                                                     gTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCGTGCGGCACCGCGGCCGTTATCACCCCGATc                                    >  1:319814/1‑150 (MQ=255)
                                                                                                         tGCCAGCCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATcggcg                                <  1:470681/150‑1 (MQ=255)
                                                                                                          gCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGAt                                                                                                            >  2:462091/1‑73 (MQ=255)
                                                                                                          gCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGAt                                                                                                            <  1:462091/73‑1 (MQ=255)
                                                                                                           ccAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCggcggc                              >  1:746570/1‑150 (MQ=255)
                                                                                                            cAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCa                                            >  2:202099/1‑135 (MQ=255)
                                                                                                            cAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCa                                            <  1:202099/135‑1 (MQ=255)
                                                                                                                cATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGa                         <  2:1003941/150‑1 (MQ=255)
                                                                                                                  tCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGGCCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGt                       <  1:427824/150‑1 (MQ=255)
                                                                                                                     cccGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTtcgt                                                                               >  2:517081/1‑91 (MQ=255)
                                                                                                                     cccGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTtcgt                                                                               <  1:517081/91‑1 (MQ=255)
                                                                                                                     cccGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACa                    >  2:855235/1‑150 (MQ=255)
                                                                                                                      ccGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCgaga                                                                                                     <  1:825941/68‑1 (MQ=255)
                                                                                                                      ccGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCgaga                                                                                                     >  2:825941/1‑68 (MQ=255)
                                                                                                                        gCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACg                 <  1:971937/150‑1 (MQ=255)
                                                                                                                          ctctCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGc                                                                                                <  2:781988/69‑1 (MQ=255)
                                                                                                                          ctctCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGc                                                                                                >  1:781988/1‑69 (MQ=255)
                                                                                                                           tctcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGaa                                                                           <  2:1035354/89‑1 (MQ=255)
                                                                                                                           tctcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGaa                                                                           >  1:1035354/1‑89 (MQ=255)
                                                                                                                           tctcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCa              >  2:33405/1‑150 (MQ=255)
                                                                                                                           tctcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCa              >  1:293027/1‑150 (MQ=255)
                                                                                                                            ctcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCaa             <  2:1146305/150‑1 (MQ=255)
                                                                                                                             tcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAg            <  1:33405/150‑1 (MQ=255)
                                                                                                                                       aaCTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACgtgt  <  1:855235/150‑1 (MQ=255)
                                                                                                                                            gcgcAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTg                                                                                          >  2:253707/1‑57 (MQ=255)
                                                                                                                                            gcgcAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTg                                                                                          <  1:253707/57‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
ACCTCGACCCGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGT  >  NC_002947/3983653‑3983937

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: