Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,543 C→A P1100Q (CCG→CAG)  rpoB → RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5430CA98.4% 178.6 / NA 62P1100Q (CCG→CAG) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (33/28);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGA  >  NC_000913/4184403‑4184672
                                                                                                                                            |                                                                                                                                 
tCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAGCCCGATCGAAGATATGCCTTACGATGAAAAGGGGAGGCCGGTAGACACCGTACTGAACCAGc                                                                                                                                 >  1:465980/1‑143 (MQ=255)
                   aTCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                              <  1:728488/143‑1 (MQ=255)
                     ccAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGa                                                                                                            >  1:643669/1‑143 (MQ=255)
                      cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                           >  2:468579/1‑143 (MQ=255)
                       aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAc                                                                                                          >  1:33667/1‑143 (MQ=255)
                           tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGCGTGTAAGTGCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCg                                                                                                      <  1:294719/143‑1 (MQ=255)
                             gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCCAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                    <  2:763040/143‑1 (MQ=255)
                               gACAAGATAACAGGGCAGCACTGTAACAAGGGTGTAATTTCTAAGATCAAGCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCAGGAAGACATCGTACGGAACCAGCAGGGGGGAACGGACCGTATGAGCAACGGACa                                                                                                  >  2:510679/1‑143 (MQ=255)
                                aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAg                                                                                                 >  1:453999/1‑143 (MQ=255)
                                aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCCCGTATGAACATCGGTccg                                                                                                 >  1:29568/1‑141 (MQ=255)
                                  aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTAAGCCGCTAGGTTTGGTTGTGGACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                                               >  1:209199/1‑143 (MQ=255)
                                   aGATGGCAGGTCGTCACGGTCACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGTCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                                              <  1:466149/143‑1 (MQ=255)
                                     gtGGCAGGTTGTCACGGTAACACGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGACGACCACGGTCCGCCGGTAGACTTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCt                                                                                            <  2:104195/142‑1 (MQ=255)
                                      gggCAGGTCGGCACGGTAACCACGTTGCAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGTTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                           <  2:101134/142‑1 (MQ=255)
                                           gggcgccACGGTAACGAGGGTGGAAGTTCTAGGGTCACCCCGATCGAAGATATGCCTTACGATGAAGACGGTACGCCGGTAGACATCGTACTGAACCAGCTGTGCGGACCGTCTCGTATGAACAGCGGTCAGATCCTCGAAAc                                                                                      <  1:656599/137‑1 (MQ=255)
                                            gtcgtgACGGTAACAAGGGTGTAATTTCTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAcc                                                                                     >  1:860577/1‑143 (MQ=255)
                                            gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTAAGATAAAGATGGTACGCCGGTAGCAATCGTACTGAACCAGCTGGGCGTACCGTCGCGTATGAAGATCGATAAGAGCCACTaaagc                                                                                     >  2:448198/1‑141 (MQ=255)
                                               gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGCGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                  >  1:275133/1‑143 (MQ=255)
                                                    ggTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGGACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCAACCTggg                                                                             >  2:836445/1‑143 (MQ=255)
                                                      tAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGACACCCACCTGGGTa                                                                           >  2:21599/1‑143 (MQ=255)
                                                       aaCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTAt                                                                          <  1:565442/143‑1 (MQ=255)
                                                       aaCAAGGGTGGACCTTCGAAGAGCGACCCGATCGCAGATACGCCTTACGATGAAAACGGGACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTAt                                                                          <  1:155982/143‑1 (MQ=255)
                                                          aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGc                                                                       >  1:760301/1‑143 (MQ=255)
                                                             ggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGGTGc                                                                    >  1:64413/1‑143 (MQ=255)
                                                             ggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGc                                                                    <  1:516255/143‑1 (MQ=255)
                                                             gccttAAGTTCGAAGATCAACCCGATGTAAGAGACTCCTTACGATGAACACGGTACGCCGGTAGACATCGTACTGGACCAGCTGGGTGTACCGTCTGGGATGAACATCGGTCAGATTCTCGAAACCCTCCTGGGTATGGCTGc                                                                    <  2:336417/139‑1 (MQ=255)
                                                              gtgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCACCTCCGCCAAACGCAACCGGGGATGCCTggg                                                                   >  2:186999/1‑141 (MQ=255)
                                                                gtAGTTTCGAAGATCGCCCCGAGCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                 <  2:616417/143‑1 (MQ=255)
                                                                gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                 <  1:468579/143‑1 (MQ=255)
                                                                gtAATTTCGAAGATCAACCCGATCGAAGAAATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTCAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGACCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                 <  1:718276/143‑1 (MQ=255)
                                                                 tAATTTTTAAAATCAACCCGATCGAACATATGCCTTACGATGAGAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                <  2:411490/143‑1 (MQ=255)
                                                                 tAATTTCTAAGATCATCCCGATCGAAGATATGCCGTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAGCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                <  2:903707/143‑1 (MQ=255)
                                                                  actttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                               <  1:564147/141‑1 (MQ=255)
                                                                  aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                               >  2:57948/1‑143 (MQ=255)
                                                                  aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                               <  1:351568/143‑1 (MQ=255)
                                                                   aTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCGGCGCAAgg                                                              >  2:292663/1‑143 (MQ=255)
                                                                        tAAGATCAACCCGATCGAAGATATGCCTTACGATGAACACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAACTCGGTCCGATCCTCGCAACCCACGTGGGTATGGCTGCGAAAGGTAGCg                                                         >  2:602819/1‑143 (MQ=255)
                                                                        tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTAAACCAGCTGGGCGTACCGTCGCGTATGAACATCGGCCAGATACTCGACACCCACCCGGGTATGGCTACGAAAgaaagtg                                                         >  2:218642/1‑137 (MQ=255)
                                                                            aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTCTGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                     <  1:301726/143‑1 (MQ=255)
                                                                            aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                     >  1:912168/1‑143 (MQ=255)
                                                                                ccccGATCGAAGATATGCCTTACGATGAAAGCGGCACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAg                                                 <  2:290382/142‑1 (MQ=255)
                                                                                 cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCCGGCTATTGCTGCGCAACGTAGCGGCGACAAGa                                                >  2:648062/1‑143 (MQ=255)
                                                                                  ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                               <  2:849478/143‑1 (MQ=255)
                                                                                         aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                        >  2:231782/1‑143 (MQ=255)
                                                                                           gATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGAGCGTACCGTCTCATATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGCAGCGGCGCCAAGGCCCACGCCag                                      >  2:104251/1‑142 (MQ=255)
                                                                                            gtatGCCTTGCGATGATGTCGGTACCGAGGTAGACATCGCACTGAACGAGTTGGGCGTACCGTCTCGTATGATCATCGGTCAGATCCTCGTACCCCCTCTGGGTATGGCGGCGAAAGGTATCGGCGACAAGATCAACGCCATg                                     <  2:281220/142‑1 (MQ=255)
                                                                                            atatGCCTTACGATGAAAACTGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATg                                     <  1:554795/143‑1 (MQ=255)
                                                                                                  cgtACGATGAATACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                               <  2:638205/141‑1 (MQ=255)
                                                                                                  cgtACGATGAAGACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                               <  2:587541/141‑1 (MQ=255)
                                                                                                  cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGCTCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCGACGCCATGCTGaaa                               >  2:711256/1‑143 (MQ=255)
                                                                                                     aCGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag                            >  1:900303/1‑143 (MQ=255)
                                                                                                           aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaa                      >  1:698397/1‑143 (MQ=255)
                                                                                                            aaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                     <  1:71287/143‑1 (MQ=255)
                                                                                                            aaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCCCGAAACCCACCTTGGTATGGCTCCGACAGGTCTCGCCGACAAGGCCAACGCCATCCTGACACAgcaacaca                     >  2:885152/1‑141 (MQ=255)
                                                                                                              aCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGACATGCCTGCGACAGGTCTCGCCGACACCCCCAACGCCATCCTGCAACAccagcacgcc                   >  2:705582/1‑139 (MQ=255)
                                                                                                                     gCCGGTAGACATCGTACTGAACCAGCTGGGCGTTCCGTCTCGTATGAACATCGGTCCGATCCTCGAAACCCACCTGGGTATGCCAGCGAGAGGTATCGGCGACAAGATCCACGCCATGCGGAAACACCAGCCAGACGTcgccg            >  2:719303/1‑141 (MQ=255)
                                                                                                                     gCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGa            >  1:163151/1‑143 (MQ=255)
                                                                                                                     gCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTGAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGa            <  1:109026/143‑1 (MQ=255)
                                                                                                                       cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCCGCAAGAAGTcgcggag          >  1:565536/1‑140 (MQ=255)
                                                                                                                            gagaTCGTACGGAACCAGCTGGGCGTGCCGTCTCGTATGAACATCGGTCAGATCCTCGAACCCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcg     <  2:643669/140‑1 (MQ=255)
                                                                                                                             aCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACACGATCAACGCCCTGCTGACACAGCAGCAAGAAGTCGCGAAACTgcgc    >  1:299516/1‑143 (MQ=255)
                                                                                                                               aTCGTACTGAACCAGCGGGGCGTACCGTCGCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGa  <  1:868866/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                 
TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGA  >  NC_000913/4184403‑4184672

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