Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,476,116 A→T I82F (ATC→TTC)  yfdT → CPS‑53 (KpLE1) prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,476,1160AT98.3% 214.3 / NA 58I82F (ATC→TTC) yfdTCPS‑53 (KpLE1) prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (29/28);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGATCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAATATCCTGATCTACCATCGCCGTCATAGAGCGTATTT  >  NC_000913/2475971‑2476262
                                                                                                                                                 |                                                                                                                                                  
cGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGttcat                                                                                                                                                <  2:186291/150‑1 (MQ=255)
 gCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCt                                                                                                                                              <  2:337791/151‑1 (MQ=255)
   aCTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTgg                                                                                                                                            >  1:516623/1‑151 (MQ=255)
   aCTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATAGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTgg                                                                                                                                            >  1:483120/1‑151 (MQ=255)
      cGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGgc                                                                                                                                           >  1:504722/1‑149 (MQ=255)
       gTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGgcagc                                                                                                                                        >  1:454348/1‑151 (MQ=255)
        tATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGgcagc                                                                                                                                        >  1:113557/1‑150 (MQ=255)
         aTCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGgcagca                                                                                                                                       >  2:489724/1‑150 (MQ=255)
               ggAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAg                                                                                                                                  <  1:381409/149‑1 (MQ=255)
               ggAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCt                                                                                                                                <  1:615514/151‑1 (MQ=255)
                 aaGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTgg                                                                                                                              <  2:30363/151‑1 (MQ=255)
                 aaGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTgg                                                                                                                              <  2:531611/151‑1 (MQ=255)
                  aGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTgg                                                                                                                              >  2:455641/1‑150 (MQ=255)
                     agagCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTAt                                                                                                                           <  2:80168/150‑1 (MQ=255)
                     agagCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATc                                                                                                                          >  2:65349/1‑151 (MQ=255)
                      gagCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATcc                                                                                                                         <  2:643927/151‑1 (MQ=255)
                        gCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCAt                                                                                                                       >  2:513746/1‑151 (MQ=255)
                          cAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCAttt                                                                                                                     <  1:450637/151‑1 (MQ=255)
                          cAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCAttt                                                                                                                     <  1:489431/151‑1 (MQ=255)
                                  cgcgAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTcaacaa                                                                                                             <  1:425460/151‑1 (MQ=255)
                                   gcgAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAAt                                                                                                            >  1:91988/1‑151 (MQ=255)
                                   gcgAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAAt                                                                                                            >  2:253911/1‑151 (MQ=255)
                                   gcgAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGATCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAAt                                                                                                            <  2:34453/151‑1 (MQ=255)
                                     gAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGCAAACAAAAAACGGCGATATAGTAAAAAGGGAAGAAGATAGAACGAGAGGTAGATACATCTGACAGCAACGAGCTGGTATCCATTTTTCAAAAATGt                                                                                                          <  1:679366/151‑1 (MQ=255)
                                     gAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGt                                                                                                          >  1:157378/1‑151 (MQ=255)
                                         ccTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTcac                                                                                                       >  2:58145/1‑150 (MQ=255)
                                            aTGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCacaac                                                                                                    >  1:59343/1‑150 (MQ=255)
                                              gAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCAcaacaa                                                                                                  <  2:281498/150‑1 (MQ=255)
                                              gAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAag                                                                                                 <  1:127698/151‑1 (MQ=255)
                                              gAATTATTCAGGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCacaaca                                                                                                   >  2:105731/1‑149 (MQ=255)
                                              gAATTAATCAAGGAAAAGCGTTTCGGTGAATTAGGAGATCGGCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCacaa                                                                                                     <  2:101491/147‑1 (MQ=255)
                                                 ttattCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAgg                                                                                              >  1:471063/1‑151 (MQ=255)
                                                 ttattCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTCCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAgg                                                                                              >  2:518944/1‑151 (MQ=255)
                                                  tattCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAgg                                                                                              <  2:349984/150‑1 (MQ=255)
                                                         ggAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGCGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAaata                                                                                       >  2:651441/1‑150 (MQ=255)
                                                         ggAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAaata                                                                                       >  2:601673/1‑150 (MQ=255)
                                                            aaaGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAaataat                                                                                     >  2:393804/1‑149 (MQ=255)
                                                                    tCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATAc                                                                           >  2:656900/1‑151 (MQ=255)
                                                                     cGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATAc                                                                           >  2:565233/1‑150 (MQ=255)
                                                                              agtagATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCacac                                                                 >  1:347398/1‑151 (MQ=255)
                                                                                  gATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGa                                                             <  2:557252/151‑1 (MQ=255)
                                                                                        cgGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGc                                                       >  1:5981/1‑151 (MQ=255)
                                                                                         ggAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGcc                                                      >  1:493772/1‑151 (MQ=255)
                                                                                                   aaaCGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTc                                            >  2:117051/1‑151 (MQ=255)
                                                                                                     aCGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCtt                                          <  2:421839/151‑1 (MQ=255)
                                                                                                      cGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTa                                         <  1:455653/151‑1 (MQ=255)
                                                                                                         cGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTaaa                                       >  1:246644/1‑150 (MQ=255)
                                                                                                             aaCGAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAatat                                   <  2:91294/150‑1 (MQ=255)
                                                                                                                gAATACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAATATCCTg                               <  2:658822/151‑1 (MQ=255)
                                                                                                                   tACATGGCATGGGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAATATCCTGATc                            >  2:85629/1‑151 (MQ=255)
                                                                                                                           atggGATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAATATCCTGATCTACCATCg                    <  2:208392/151‑1 (MQ=255)
                                                                                                                             gggATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAATATCCTGATCTACCATCGcc                  <  2:141356/151‑1 (MQ=255)
                                                                                                                             gggATATGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAATATCCTGATCTACCATCGc                   <  1:85641/150‑1 (MQ=255)
                                                                                                                                atatGACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAATATCCTGATCTACCATCGCCGTc               <  2:269414/151‑1 (MQ=255)
                                                                                                                                    gACTCTCGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAATATCCTGATCTACCATCGCCGTCATag           <  1:269679/151‑1 (MQ=255)
                                                                                                                                         tcGGTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAATATCCTGATCTACCATCGCCGTCATAGAGCGt      >  2:477441/1‑151 (MQ=255)
                                                                                                                                            gTTGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAATATCCTGATCTACCATCGCCGTCATAGAGCGTAtt   <  1:8825/151‑1 (MQ=255)
                                                                                                                                             ttGGTTCATCTGGCAGCAACGAGCTAGTATCCATTTTTCAACAATGTCACAACACGAGGTGAAATAATAGAGTCATACAGCCTCACTCTCGATGAGACCTGTCAGTTTCTTAAAATATCCTGATCTACCATCGCCGTCATAGAGCGTAtt   <  2:289753/150‑1 (MQ=255)
                                                                                                                                              tGGTTCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAATATCCTGATCTACCATCGCCGTCATAGAGCGTAttt  <  1:450736/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGCACTCGTATCTCTGGAAGTAGAGCCAGATGAACGCGAAGCCTATGAATTATTCATGGAAAAGCGTTTCGGTGACTTAGTAGATCGTCGGAGAGCAAAAAACGGCGATAACGAATACATGGCATGGGATATGACTCTCGGTTGGATCATCTGGCAGCAACGAGCTGGTATCCATTTTTCAACAATGTCACAACAAGAGGTGAAATAATAGAGCCATACAGCCTCACACTCGATGAGGCCTGTCAGTTTCTTAAAATATCCTGATCTACCATCGCCGTCATAGAGCGTATTT  >  NC_000913/2475971‑2476262

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