Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 70,289 G→T intergenic (‑241/‑98) araB ← / → araC L‑ribulokinase/transcriptional regulator for ara operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091370,2890GT98.5% 400.6 / NA 134intergenic (‑241/‑98)araB/araCL‑ribulokinase/transcriptional regulator for ara operon
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (67/65);  total (67/67)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

TTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTT  >  NC_000913/70216‑70360
                                                                          |                                                                       
ttGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCCGGTTgg                                                                        >  1:687926/1‑76 (MQ=255)
 tGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGt                                                                      >  2:1711966/1‑76 (MQ=255)
  gTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGtt                                                                     <  1:5001399/76‑1 (MQ=255)
    tACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTCATAATCGGGGTTACCGGTTTGGTTAg                                                                   <  2:5364853/76‑1 (MQ=255)
    tACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAg                                                                   <  1:5613726/76‑1 (MQ=255)
    tACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAg                                                                   <  1:2186573/76‑1 (MQ=255)
    tACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGGTTGGTTAg                                                                   >  2:844076/1‑76 (MQ=255)
     aCGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGc                                                                  <  1:1370755/76‑1 (MQ=255)
     aCGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGc                                                                  >  1:2098342/1‑76 (MQ=255)
      cgcgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCg                                                                 >  1:5325710/1‑76 (MQ=255)
      cgcgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCg                                                                 <  2:2646939/76‑1 (MQ=255)
      cgcgTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCTTTTAATACGCGGGGTTCCCGGTTTGGTTAGCg                                                                 <  1:5558648/76‑1 (MQ=255)
        cgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCgag                                                               >  2:3433423/1‑76 (MQ=255)
          tttttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaa                                                             >  2:554764/1‑76 (MQ=255)
          tttttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaa                                                             <  1:1464204/76‑1 (MQ=255)
          tttttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaa                                                             <  1:2178791/76‑1 (MQ=255)
           ttttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaag                                                            <  2:5094214/76‑1 (MQ=255)
             ttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAagag                                                          >  2:4380819/1‑76 (MQ=255)
             ttGTCAGGGCTTTGGACCCGATTTGTTACAGAATGATTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAagag                                                          >  1:3405236/1‑76 (MQ=255)
              tGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGc                                                         <  1:5345504/76‑1 (MQ=255)
               gTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGcc                                                        <  1:2097081/76‑1 (MQ=255)
               gTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGcc                                                        <  1:32550/76‑1 (MQ=255)
               gTCATGGCGGAGGTCCCGCTGTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGcc                                                        <  1:46287/76‑1 (MQ=255)
                tCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCa                                                       <  2:4415968/76‑1 (MQ=255)
                tCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCa                                                       <  1:2417979/76‑1 (MQ=255)
                tCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCa                                                       >  2:3452339/1‑76 (MQ=255)
                tCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCa                                                       <  2:2508880/76‑1 (MQ=255)
                 cATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAg                                                      >  2:4909112/1‑76 (MQ=255)
                 cATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAg                                                      >  1:2892032/1‑76 (MQ=255)
                  aTGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaaga                                                           <  2:1133825/70‑1 (MQ=255)
                  aTGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaaga                                                           >  1:1133825/1‑70 (MQ=255)
                  aTGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGt                                                     >  2:4925176/1‑76 (MQ=255)
                  aTGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGt                                                     >  1:3743210/1‑76 (MQ=255)
                    ggCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTaa                                                   <  1:2235270/76‑1 (MQ=255)
                    ggCTTTGGGCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTaa                                                   <  2:1316960/76‑1 (MQ=255)
                     gCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTaaa                                                  >  2:1584155/1‑76 (MQ=255)
                     gCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTaaa                                                  >  1:5511135/1‑76 (MQ=255)
                      cTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTaaaa                                                 >  2:5376925/1‑76 (MQ=255)
                       tttGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAg                                                >  1:4904805/1‑76 (MQ=255)
                         tGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGAc                                              <  1:1073689/76‑1 (MQ=255)
                         tGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGAc                                              <  1:3085299/76‑1 (MQ=255)
                          ggTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACg                                             >  1:3368364/1‑76 (MQ=255)
                          gagcccGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGCTTGGGTAGCGAGAAGAGCCAGTAAAAGACg                                             >  2:5137168/4‑76 (MQ=255)
                           gTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGc                                            <  2:4494485/76‑1 (MQ=255)
                             cccGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAg                                          >  1:5379454/1‑76 (MQ=255)
                              ccGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAACGCAGCAGt                                         >  1:4235097/1‑76 (MQ=255)
                              ccGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGt                                         >  1:3241709/1‑76 (MQ=255)
                                gtgttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGa                                       <  2:2859682/73‑1 (MQ=255)
                                gCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGa                                       >  1:3311184/1‑76 (MQ=255)
                                  tttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACg                                     <  1:2012032/76‑1 (MQ=255)
                                  tttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACg                                     >  1:2170415/1‑76 (MQ=255)
                                  tttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACg                                     >  2:1843125/1‑76 (MQ=255)
                                  tttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACg                                     >  1:1579702/1‑76 (MQ=255)
                                   ttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACgg                                    >  1:5422856/1‑76 (MQ=255)
                                    tGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCGAGTAAAAGACGCAGTGACGGc                                   <  2:4095628/76‑1 (MQ=255)
                                    tGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGc                                   <  1:605915/76‑1 (MQ=255)
                                    tGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGc                                   <  1:4598877/76‑1 (MQ=255)
                                    tGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGc                                   <  2:310140/76‑1 (MQ=255)
                                     gTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCa                                  >  2:4977377/1‑76 (MQ=255)
                                      ttACAGAATGCTTTTAATAAGCGGGGTTAACGGTTTGGCTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCaa                                 >  1:3638042/1‑76 (MQ=255)
                                       tACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAAt                                >  2:4726617/1‑76 (MQ=255)
                                       tACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAAt                                >  2:599344/1‑76 (MQ=255)
                                        aCAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATg                               >  1:2525661/1‑76 (MQ=255)
                                        aCAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATg                               >  1:4732616/1‑76 (MQ=255)
                                         cAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGt                              >  1:5384580/1‑76 (MQ=255)
                                         cAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGt                              >  1:346424/1‑76 (MQ=255)
                                           gAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCt                            <  1:421402/76‑1 (MQ=255)
                                           gAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCt                            >  2:5238435/1‑76 (MQ=255)
                                           gAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCt                            >  2:3428832/1‑76 (MQ=255)
                                           aaaTGTTTTTAAGAAGGGGGGTTACGGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCt                            <  2:5739897/75‑1 (MQ=255)
                                            aaTGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTg                           <  1:5498098/76‑1 (MQ=255)
                                             aTGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGa                          >  1:3219188/1‑76 (MQ=255)
                                              tGCTTTTAATAAGCGGGGTTGCCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGAt                         <  2:2344069/76‑1 (MQ=255)
                                              tGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGAt                         >  1:4911875/1‑76 (MQ=255)
                                              tGCTTGTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGAt                         <  2:1378386/76‑1 (MQ=255)
                                               gCTTTTAATAAGCGGGGTTACGGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATg                        <  2:4426726/76‑1 (MQ=255)
                                                cTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGc                       <  1:3780350/76‑1 (MQ=255)
                                                 ttttAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCa                      <  2:3288533/76‑1 (MQ=255)
                                                    taataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAata                   <  1:1342190/76‑1 (MQ=255)
                                                     aatagGCGGGGTTGCCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAatat                  <  2:2072730/76‑1 (MQ=255)
                                                     aataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAatat                  >  2:855966/1‑76 (MQ=255)
                                                     aataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAatat                  <  1:720199/76‑1 (MQ=255)
                                                      ataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGGGACGGAAATGTCTGATGCAATATg                 >  1:4147139/1‑76 (MQ=255)
                                                         aGCGGGGTTACGGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGAc              <  2:2413974/76‑1 (MQ=255)
                                                         aGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGAc              >  1:2855895/1‑76 (MQ=255)
                                                          gCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACa             <  1:1041528/76‑1 (MQ=255)
                                                           cGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACaa            <  1:2411789/76‑1 (MQ=255)
                                                            ggggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAAt           >  2:2351073/1‑76 (MQ=255)
                                                             gggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAAtt          >  1:179378/1‑76 (MQ=255)
                                                              ggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTg         >  1:4354709/1‑76 (MQ=255)
                                                              ggTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTg         <  1:934139/76‑1 (MQ=255)
                                                                ttACCGGTTTGGTTAGCGAGAGGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGt       <  2:2706215/76‑1 (MQ=255)
                                                                ttACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGCCAATTGGt       >  2:4827406/1‑76 (MQ=255)
                                                                ttACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGt       >  1:2228059/1‑76 (MQ=255)
                                                                 tACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGtt      <  2:2202579/76‑1 (MQ=255)
                                                                  aCCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGc                       <  1:1557191/58‑1 (MQ=255)
                                                                  aCCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGttt     >  2:4850618/1‑76 (MQ=255)
                                                                    cgggTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTttc    <  2:3572421/74‑1 (MQ=255)
                                                                     cGGTTTGGTTAGCGAGGAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTttc  <  2:210868/76‑2 (MQ=255)
                                                                      ggTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACAGCAATGTCTGATGCAATATGGACAATTGGTttctt  >  1:478054/1‑76 (MQ=255)
                                                                          |                                                                       
TTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTT  >  NC_000913/70216‑70360

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: