Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,984,043 A→C 8.5% intergenic (+37/+153) ilvE → / ← PP_3512 branched‑chain‑amino acid aminotransferase/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,984,0430AC 8.5% 173.2 / 7.0 59intergenic (+37/+153)ilvE/PP_3512branched‑chain‑amino acid aminotransferase/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/23);  new base C (2/3);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

CCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGG  >  NC_002947/3983895‑3984186
                                                                                                                                                    |                                                                                                                                               
ccACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCaa                                                                                                                                                 >  5:234405/1‑149 (MQ=255)
  aCGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATg                                                                                                                                               <  4:276351/149‑1 (MQ=255)
  aCGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATg                                                                                                                                               <  1:249968/149‑1 (MQ=255)
    gACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATgg                                                                                                                                              >  5:407783/1‑148 (MQ=255)
       cTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATg                                                                                                                                          >  7:22307/1‑149 (MQ=255)
        tGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGc                                                                                                                                         >  1:416319/1‑149 (MQ=255)
         ggAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGcc                                                                                                                                        >  2:226602/1‑149 (MQ=255)
          ttccaGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGc                                                                                                                                         <  3:241243/143‑1 (MQ=255)
          gAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGcc                                                                                                                                        <  1:412531/148‑1 (MQ=255)
           aaaaGGTCGGCCCGGTTACCCATAAGCTCTACAACGATCTGACCGGTATCCAGAGAGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGcccc                                                                                                                                      <  4:258445/149‑1 (MQ=255)
               ggTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGttt                                                                                                                                  <  4:60840/149‑1 (MQ=255)
               ggTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCAGACCGCTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGACTGAAACGCATCAGATGAATGAAAACGGGGCATGCCAGCAAAGGCATGCACCGttt                                                                                                                                  >  4:90198/1‑149 (MQ=255)
                   ggCCCGGTTAGCCAGAAGGTCTACAACGAGGTGTCGTGTTTGCCGAGCGGCGATGTAGAGGCACCGGCAGGCTGGAGCGTCTAGGTCGCGTGATACGCATCAGCTGTATGAAAACGGGGCATGCCAGCTATGGCATGCCCCGtttttt                                                                                                                               <  3:90198/148‑1 (MQ=255)
                     cccGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGc                                                                                                                                             >  8:238493/1‑132 (MQ=255)
                     cccGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGc                                                                                                                                             <  7:238493/132‑1 (MQ=255)
                       cGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGGGTTGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGGATCAGCTGTATGAAAACGGGGAATGCCAGCAATGGCATGCCCCGTTTTTTCAtgt                                                                                                                          >  5:349721/1‑149 (MQ=255)
                             cccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTac                                                                                                                    <  7:28961/149‑1 (MQ=255)
                              ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTaca                                                                                                                   <  1:178119/149‑1 (MQ=255)
                                 gAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACgg                                                                                                                >  6:34908/1‑149 (MQ=255)
                                         aCAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCg                                                                                                        >  7:348158/1‑149 (MQ=255)
                                          cAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTGATGTGTTTACACGGCCCCCTCGt                                                                                                       <  5:95985/149‑1 (MQ=255)
                                            aCGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAg                                                                                                     <  2:325193/149‑1 (MQ=255)
                                             cGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCACCGTTTTTTCATGTGTTACACGGCCCCCTCGTAg                                                                                                     >  4:397695/1‑147 (MQ=255)
                                               aGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAg                                                                                                  >  5:3640/1‑149 (MQ=255)
                                               aGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAg                                                                                                  <  2:416319/149‑1 (MQ=255)
                                                gCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGc                                                                                                 <  8:22307/149‑1 (MQ=255)
                                                      cGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGtt                                                                                            >  7:235915/1‑148 (MQ=255)
                                                        gTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTc                                                                                                         >  4:148264/1‑133 (MQ=255)
                                                        gTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTc                                                                                                         <  3:148264/133‑1 (MQ=255)
                                                        gTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTAc                                                                                         >  6:306380/1‑149 (MQ=255)
                                                              agagCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGa                                                                                   <  6:407783/149‑1 (MQ=255)
                                                                    gcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                             <  4:410682/149‑1 (MQ=255)
                                                                         gTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCacac                                                                        >  1:149308/1‑149 (MQ=255)
                                                                               gCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCgc                                                                    >  7:45999/1‑147 (MQ=255)
                                                                                cACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCgccag                                                                 >  5:384644/1‑149 (MQ=255)
                                                                                      cAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCaata                                                           <  6:234405/149‑1 (MQ=255)
                                                                                        ggCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCaataat                                                         <  1:226602/149‑1 (MQ=255)
                                                                                                 gTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGga                                                 <  5:306380/148‑1 (MQ=255)
                                                                                                    aaGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCa                                            <  3:397695/149‑1 (MQ=255)
                                                                                                          gCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGtt                                       >  8:209624/1‑149 (MQ=255)
                                                                                                            cTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGga                                                 <  4:236628/137‑1 (MQ=255)
                                                                                                            cTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGga                                                 >  3:236628/1‑137 (MQ=255)
                                                                                                            cTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACACGCGAATAGGACACACCCGCCAGCCAATAATGCCTGGGAGATAACAGTTg                                      >  8:96627/1‑148 (MQ=255)
                                                                                                               aaaCGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAg                                  <  6:3640/149‑1 (MQ=255)
                                                                                                                             tATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAg                    >  1:255891/1‑149 (MQ=255)
                                                                                                                              aTGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGc                   >  4:152046/1‑149 (MQ=255)
                                                                                                                                   aaCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGc               >  5:316308/1‑148 (MQ=255)
                                                                                                                                    aCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGc               >  5:255679/1‑147 (MQ=255)
                                                                                                                                      ggggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCTAATAGGCCACACCCGCCCGCCAATAATGCCTGGGAGATCACCGTTGTCTGTGCAGCCGCATTTGCGGGCAAGc           <  2:149308/149‑1 (MQ=255)
                                                                                                                                       gggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGcc          >  4:344585/1‑149 (MQ=255)
                                                                                                                                        ggCATGCCATTTCTGGCATGACCCGTTTTTTCATGTGTTTACACGGCCCCCTTGTTGGAGTGTGTTTACTCGCGAATAGGCTACACTCGCCAGCTAATAATGTCTGGGAGATCATTGTTGTTAGTGCa                              <  8:64955/128‑1 (MQ=255)
                                                                                                                                        ggCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGccc         >  7:5000/1‑149 (MQ=255)
                                                                                                                                        ggCATGCCATTGCTGGCATGCCCCGTTTTTTAATGTGTTTACAAGGCCCCCTCGTAGGAGCGGGTTTACCAGAGAATAGGCCACACCCGCAAGAAAAAAATGCCTGGGAGAACACCGTTGTCAGTgaa                              >  7:64955/1‑126 (MQ=255)
                                                                                                                                        ggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGccc         >  4:394060/1‑149 (MQ=255)
                                                                                                                                         gCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCcg        >  5:177984/1‑149 (MQ=255)
                                                                                                                                              ccAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCg   >  7:55491/1‑149 (MQ=255)
                                                                                                                                               cAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGTCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTACAGCCGCATTAGCGGGCAAGCCCGCGCCgg  <  8:235915/149‑1 (MQ=255)
                                                                                                                                               cAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCgg  >  3:320511/1‑149 (MQ=255)
                                                                                                                                               cAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCg   >  1:157713/1‑148 (MQ=255)
                                                                                                                                               cAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCATGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGTCAAGCCCGCGCCgg  >  8:180872/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGG  >  NC_002947/3983895‑3984186

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: