Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,165,628 T→A L223Q (CTG→CAG)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,165,6280TA99.2% 1100.9 / ‑1.4 261L223Q (CTG→CAG) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (121/138);  total (121/140)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 265 total aligned reads displayed.

GAACAGTTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAG  >  REL606/4165555‑4165702
                                                                         |                                                                          
gAACAGTTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGt                                                                          <  2:721005/76‑1 (MQ=255)
gAACAGTTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGt                                                                          <  1:4097447/76‑1 (MQ=255)
 aaCAGTTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGtt                                                                         <  1:2548946/76‑1 (MQ=255)
  aCAGTTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTg                                                                        <  2:4233913/76‑1 (MQ=255)
   cAGTTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTgg                                                                       <  1:633875/76‑1 (MQ=255)
    aGTTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGa                                                                      >  2:2687857/1‑76 (MQ=255)
    aGTTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGa                                                                      <  1:1123314/76‑1 (MQ=255)
    aGTTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGa                                                                      >  1:3443901/1‑76 (MQ=255)
     gTTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGaa                                                                     <  2:565039/76‑1 (MQ=255)
       tGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGc                                                                   >  2:2737061/1‑76 (MQ=255)
        gCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGcg                                                                  >  1:2279742/1‑76 (MQ=255)
        gCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGcg                                                                  <  2:6444542/76‑1 (MQ=255)
         cGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGcgt                                                                 <  1:6960191/76‑1 (MQ=255)
          gTGAAGAGCTGAACGGAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGcgtt                                                                >  1:3912769/1‑76 (MQ=255)
          gTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGcgtt                                                                <  2:638032/76‑1 (MQ=255)
           tGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGcgttc                                                               >  1:5750118/1‑76 (MQ=255)
            gAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGcgttcg                                                              >  1:475357/1‑76 (MQ=255)
            gAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGcgttcg                                                              <  1:3225220/76‑1 (MQ=255)
             aaGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGcgttcgt                                                             >  1:4896437/1‑76 (MQ=255)
             aaGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGcgttcgt                                                             <  1:283099/76‑1 (MQ=255)
             aaGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGcgttcgt                                                             <  2:6843829/76‑1 (MQ=255)
              agagCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGcgttcgtt                                                            <  2:3472958/76‑1 (MQ=255)
               gagCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCgttcgttc                                                           <  1:824389/76‑1 (MQ=255)
                  cTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGt                                                        <  2:3150803/76‑1 (MQ=255)
                  cTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGt                                                        <  1:2581565/76‑1 (MQ=255)
                    gAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCt                                                      >  1:3001693/1‑76 (MQ=255)
                    gAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCt                                                      <  1:1817072/76‑1 (MQ=255)
                    aaaCGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCt                                                      <  1:3993525/75‑1 (MQ=255)
                     aaCGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTg                                                     >  1:5727273/1‑76 (MQ=255)
                     aaCGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTg                                                     <  1:968901/76‑1 (MQ=255)
                     aaCGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTg                                                     >  1:1172168/1‑76 (MQ=255)
                       cGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGt                                                   <  2:1979862/76‑1 (MQ=255)
                       cGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGt                                                   <  2:5384302/76‑1 (MQ=255)
                        gAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTa                                                  <  1:5135682/76‑1 (MQ=255)
                        gAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTa                                                  <  2:2768782/76‑1 (MQ=255)
                        gAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTa                                                  <  1:420187/76‑1 (MQ=255)
                        gAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTa                                                  <  2:3410884/76‑1 (MQ=255)
                        gAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTa                                                  <  2:5728377/76‑1 (MQ=255)
                        gAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTa                                                  >  1:1751337/1‑76 (MQ=255)
                         aaaCCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTaa                                                 <  2:1175904/76‑1 (MQ=255)
                         aaaCCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTaa                                                 <  1:1200318/76‑1 (MQ=255)
                          aaCCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAAc                                                >  1:776343/1‑76 (MQ=255)
                           aCCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACa                                               <  1:6983699/76‑1 (MQ=255)
                           aCCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACa                                               <  2:5456318/76‑1 (MQ=255)
                           aCCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACa                                               <  2:6819944/76‑1 (MQ=255)
                            ccAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACaa                                              <  1:6270879/76‑1 (MQ=255)
                              aaCTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAAc                                            >  2:490031/1‑76 (MQ=255)
                              aaCTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAAc                                            >  2:5241489/1‑76 (MQ=255)
                                cTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCa                                          <  1:4534917/76‑1 (MQ=255)
                                 tCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCag                                         <  1:6732662/76‑1 (MQ=255)
                                 tCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCag                                         <  2:1518607/76‑1 (MQ=255)
                                 tCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCag                                         <  2:1838683/76‑1 (MQ=255)
                                   cGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCagag                                       >  2:4207472/1‑76 (MQ=255)
                                   cGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCagag                                       >  2:3848293/1‑76 (MQ=255)
                                   cGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCagag                                       >  2:6038664/1‑76 (MQ=255)
                                     aaaCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTg                                     >  1:894943/1‑76 (MQ=255)
                                     aaaCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTg                                     >  2:5043931/1‑76 (MQ=255)
                                     aaaCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTg                                     <  1:3795669/76‑1 (MQ=255)
                                      aaCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTgg                                    >  2:4346977/1‑76 (MQ=255)
                                      aaCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTgg                                    <  2:824102/76‑1 (MQ=255)
                                      aaCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTgg                                    <  2:4626199/76‑1 (MQ=255)
                                         cAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGgatg                                 <  1:250118/76‑1 (MQ=255)
                                         cAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGgatg                                 >  1:1623314/1‑76 (MQ=255)
                                         cAAGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGgatg                                 <  1:5425950/76‑1 (MQ=255)
                                           aGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGgatgat                               <  2:970260/76‑1 (MQ=255)
                                           aGCGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGgatgat                               <  2:3805753/76‑1 (MQ=255)
                                             cGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATcc                             >  2:5091055/1‑76 (MQ=255)
                                             cGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATcc                             >  1:3839505/1‑76 (MQ=255)
                                             cGTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATcc                             >  1:3739299/1‑76 (MQ=255)
                                              gTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCt                            <  2:165341/76‑1 (MQ=255)
                                              gTAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCt                            >  2:3204734/1‑76 (MQ=255)
                                               tAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTg                           >  2:4525364/1‑76 (MQ=255)
                                               tAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTg                           >  1:6102500/1‑76 (MQ=255)
                                               tAAAAAGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTg                           <  2:176052/76‑1 (MQ=255)
                                                aaaaaGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGa                          <  1:6409069/76‑1 (MQ=255)
                                                aaaaaGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGa                          >  2:1946711/1‑76 (MQ=255)
                                                aaaaaGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGa                          <  1:3318752/76‑1 (MQ=255)
                                                 aaaaGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGAc                         >  2:4387667/1‑76 (MQ=255)
                                                  aaaGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGCGTGGATGATCCTGAcc                        >  2:2350993/1‑76 (MQ=255)
                                                  aaaGCTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGAcc                        <  2:2320861/76‑1 (MQ=255)
                                                      cTGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTc                    >  1:6212208/1‑76 (MQ=255)
                                                       tGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCt                   <  1:6609909/76‑1 (MQ=255)
                                                       tGACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCt                   >  1:172685/1‑76 (MQ=255)
                                                        gACCAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTg                  >  2:6895215/1‑76 (MQ=255)
                                                           cAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCg               >  1:1392184/1‑76 (MQ=255)
                                                           cAAGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCg               >  2:1354584/1‑76 (MQ=255)
                                                             aGCGTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGt             <  1:5978427/76‑1 (MQ=255)
                                                                gTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACt          <  2:4822850/76‑1 (MQ=255)
                                                                gTATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACt          <  1:202547/76‑1 (MQ=255)
                                                                 tATCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTg         >  1:5119512/1‑76 (MQ=255)
                                                                  aTCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTgc        >  1:4441841/1‑76 (MQ=255)
                                                                  aTCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTgc        >  2:885314/1‑76 (MQ=255)
                                                                   tCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTgcc       >  1:5243160/1‑76 (MQ=255)
                                                                   tCAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTgcc       >  2:4837593/1‑76 (MQ=255)
                                                                    cAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTgccg      >  2:1294020/1‑76 (MQ=255)
                                                                    cAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTgccg      >  2:228008/1‑76 (MQ=255)
                                                                    cAAACAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTgccg      >  2:315824/1‑76 (MQ=255)
                                                                     aaaCAGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTgccgc     >  2:3688945/1‑76 (MQ=255)
                                                                      aaCGGTTGGAACCTGTCGTTCACCTTGGTAACAAGCCAGAGTGGATAATCCTGACCGTTCTGCCGGTACTgccgcc    <  2:4135986/76‑1 (MQ=255)
                                                                        cTGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAg  <  1:2160981/76‑1 (MQ=255)
                                                                         |                                                                          
GAACAGTTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGTTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAG  >  REL606/4165555‑4165702

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: