Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,999,387 Δ1 bp coding (1405/2163 nt) uvrD → DNA‑dependent ATPase I and helicase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,999,3870G.100.0% 274.3 / NA 57coding (1405/2163 nt)uvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (32/25);  total (32/25)

CTGGACGTGGTACGTCAGACATCGCGCGATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTAGCGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAACGCG  >  NC_000913/3999249‑3999532
                                                                                                                                          |                                                                                                                                                 
cTGGACGTGGTACGTCAGACATCGCGCGATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGc                                                                                                                                      >  2:170478/1‑151 (MQ=255)
  ggACGTGGTACGTCAGACATCGCGCGATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAGACTGCCg                                                                                                                                    <  2:550997/151‑1 (MQ=255)
  ggACGTGGTACGTCAGACATCGCGCGATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCg                                                                                                                                    >  1:295006/1‑151 (MQ=255)
           aCGTCAGACATCGCGCGATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGc                                                                                                                           >  1:931010/1‑151 (MQ=255)
               cAGACATCGCGCGATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCa                                                                                                                       <  1:647810/151‑1 (MQ=255)
                 gACATCGCGCGATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCATCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCAt                                                                                                                      <  1:921142/150‑1 (MQ=255)
                  aCATCGCGCGATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGt                                                                                                                    <  1:226392/151‑1 (MQ=255)
                  aCATCGCGCGATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGt                                                                                                                    >  1:377946/1‑151 (MQ=255)
                      cgcgcgATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAg                                                                                                                <  2:645055/151‑1 (MQ=255)
                      cgcgcgATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAg                                                                                                                <  2:24354/151‑1 (MQ=255)
                       gcgcgATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTAc                                                                                                                  >  1:733272/1‑148 (MQ=255)
                         gcgATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACt                                                                                                             <  2:813660/151‑1 (MQ=255)
                          cgATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTg                                                                                                            >  1:979219/1‑151 (MQ=255)
                                  aGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGt                                                                                                     >  1:241542/1‑150 (MQ=255)
                                  aGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGAGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGt                                                                                                     >  1:647039/1‑150 (MQ=255)
                                  aGTTAACACTCAGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGt                                                                                                     >  2:889492/1‑150 (MQ=255)
                                       acacTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTaa                                                                                               <  1:1011914/151‑1 (MQ=255)
                                        cacTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGCGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTaaa                                                                                              <  1:285113/151‑1 (MQ=255)
                                                    caTGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTg                                                                                  <  1:995227/151‑1 (MQ=255)
                                                     aTGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGc                                                                                 >  2:283927/1‑151 (MQ=255)
                                                     aTGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGc                                                                                 >  1:96755/1‑151 (MQ=255)
                                                            gAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCAt                                                                          <  2:16298/151‑1 (MQ=255)
                                                                tGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTAt                                                                      >  2:267691/1‑151 (MQ=255)
                                                                   tGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATg                                                                     >  1:45358/1‑149 (MQ=255)
                                                                    gCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAg                                                                   <  2:464533/150‑1 (MQ=255)
                                                                      aGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAg                                                                <  1:146263/151‑1 (MQ=255)
                                                                        gAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGga                                                              <  2:456081/151‑1 (MQ=255)
                                                                         aaaaaGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGgag                                                             <  2:1040748/151‑1 (MQ=255)
                                                                          aaaaGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAgg                                                               >  2:462203/1‑148 (MQ=255)
                                                                              gCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAggg                                                        <  2:547111/151‑1 (MQ=255)
                                                                                   cGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGaa                                                    >  2:319795/1‑150 (MQ=255)
                                                                                   cGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGaa                                                    >  2:305152/1‑150 (MQ=255)
                                                                                       gggCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAgg                                               <  1:85295/151‑1 (MQ=255)
                                                                                        ggCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAgg                                               >  1:546236/1‑150 (MQ=255)
                                                                                         gCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTc                                             >  2:764174/1‑151 (MQ=255)
                                                                                          cGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCa                                            <  2:123900/151‑1 (MQ=255)
                                                                                           gAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTc                                             >  2:868699/1‑149 (MQ=255)
                                                                                                 ccAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGt                                     >  2:316412/1‑151 (MQ=255)
                                                                                                  cAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGt                                     >  2:635847/1‑150 (MQ=255)
                                                                                                     cgcCTTGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCg                                 >  1:504452/1‑151 (MQ=255)
                                                                                                          tGCAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAAc                            >  2:770537/1‑151 (MQ=255)
                                                                                                            cAGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACtt                          <  1:562352/151‑1 (MQ=255)
                                                                                                             aGCGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACtt                          >  2:781299/1‑150 (MQ=255)
                                                                                                               cGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACtt                          >  2:456050/1‑148 (MQ=255)
                                                                                                               cGATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTaga                       <  1:173230/151‑1 (MQ=255)
                                                                                                                gATTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTag                        >  2:542239/1‑149 (MQ=255)
                                                                                                                 aTTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAgg                     <  2:48488/151‑1 (MQ=255)
                                                                                                                 aTTTATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAgg                     <  1:589399/151‑1 (MQ=255)
                                                                                                                  tttATGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAgg                     >  1:361272/1‑150 (MQ=255)
                                                                                                                     aTGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACt                 >  1:28562/1‑151 (MQ=255)
                                                                                                                     aTGGAATTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACt                 >  1:913380/1‑151 (MQ=255)
                                                                                                                         aaTTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTg             <  1:137248/151‑1 (MQ=255)
                                                                                                                          aTTAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTg             >  2:559249/1‑150 (MQ=255)
                                                                                                                            tAATCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACg          >  1:826438/1‑151 (MQ=255)
                                                                                                                              aTCGACGCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGc        <  1:170416/151‑1 (MQ=255)
                                                                                                                                    gCCTTA‑CGCAGGAAACTGCCGATATGCCTCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAAcgcg  >  1:382171/1‑151 (MQ=255)
                                                                                                                                    gCCTTA‑CGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAAcgcg  <  2:707120/151‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
CTGGACGTGGTACGTCAGACATCGCGCGATCGCCAGTTAACACTCTGGCAGGCATGTCGTGAGCTGTTGCAGGAAAAAGCCCTCGCCGGGCGAGCTGCCAGCGCCTTGCAGCGATTTATGGAATTAATCGACGCCTTAGCGCAGGAAACTGCCGATATGCCGCTGCATGTACAGACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAACGCG  >  NC_000913/3999249‑3999532

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑