Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,814,273 G→A 6.8% G136D (GGT→GAT)  katE → catalase HPII, heme d‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,814,2730GA 6.8% 179.9 / 2.7 59G136D (GGT→GAT) katEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base G (25/30);  new base A (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

CCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTT  >  NC_000913/1814132‑1814421
                                                                                                                                             |                                                                                                                                                    
ccGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTc                                                                                                                                               <  3:110165/149‑1 (MQ=255)
 cgacgaTCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTcc                                                                                                                                              <  8:91725/149‑1 (MQ=255)
  gacgaTCAAAACTCACGGCGTGCCGGTAGCCGTGGTCCAACGCTGCGGGAAGATTTTATGCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCa                                                                                                                                             <  3:196137/149‑1 (MQ=255)
     gaTCAAAACTCACTGAGTGCCGGTAGCAGTCGTCCAACTCTGCTCGAAGATTATATTCTGCGCGAGAAAATCACCCACTTTGACAATGAGCGCATTCCGGAACGTAATCTTAATGCACGCGCATCACCCGCTCACGGTTATTTCCAGc                                                                                                                                           >  8:34157/1‑148 (MQ=255)
      aTCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCa                                                                                                                                         <  7:48293/149‑1 (MQ=255)
             cTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCTGATCAGCCGCTCACGGTTATTTCCAGCCATATaaaa                                                                                                                                  <  2:34097/149‑1 (MQ=255)
             cTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATaaaa                                                                                                                                  <  3:173825/149‑1 (MQ=255)
                         cGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACTGTTATTTCCAGCCATATAAAAGCTTAAGCGata                                                                                                                      <  2:207136/149‑1 (MQ=255)
                          ggTAGCCGTGGTCCAACTCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGatat                                                                                                                     >  4:189524/1‑149 (MQ=255)
                              gCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCTGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCCGCCAGATAAAAGCTTAAGCGATATTAc                                                                                                                  >  5:129169/1‑148 (MQ=255)
                              gCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTAc                                                                                                                  >  5:83876/1‑148 (MQ=255)
                               ccGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTAcc                                                                                                                 >  4:626/1‑148 (MQ=255)
                                    gTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGc                                                                                                            <  8:33553/148‑1 (MQ=255)
                                      ccAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGGATTGTTCATGCACGCGGATCAGCCGCTCACGGTTAGTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGa                                                                                                         >  5:228345/1‑149 (MQ=255)
                                                tgGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCctct                                                                                                 >  7:132641/1‑147 (MQ=255)
                                                  gAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAg                                                                                                                           >  6:11967/1‑119 (MQ=255)
                                                  gAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAg                                                                                                                           <  5:11967/119‑1 (MQ=255)
                                                   aaGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGATTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCctct                                                                                                 >  7:50983/1‑144 (MQ=255)
                                                   aaGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGATTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCctct                                                                                                 >  3:4154/1‑144 (MQ=255)
                                                           aTTCTGCGCGAGAAAATCACCCACTTTGCCCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATATCCTCTCAGATCCGAACa                                                                                     >  7:96202/1‑148 (MQ=255)
                                                           aTTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACa                                                                                     >  1:171695/1‑148 (MQ=255)
                                                            ttCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACaaa                                                                                   <  3:626/149‑1 (MQ=255)
                                                                 cgcgAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCa                                                                               <  8:206802/148‑1 (MQ=255)
                                                                  gcgAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAcc                                                                             >  6:80047/1‑149 (MQ=255)
                                                                    gagaAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAccc                                                                            >  6:72865/1‑148 (MQ=255)
                                                                    gagaAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAccc                                                                            >  6:62586/1‑148 (MQ=255)
                                                                     agaAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCa                                                                          <  7:17728/149‑1 (MQ=255)
                                                                         aaTCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTAt                                                                      <  4:165285/149‑1 (MQ=255)
                                                                            cACCCACTTTGACCATGATCGCATTCCGTAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCATCCATATAAAAGCTTAAGCGTTATTACCAAAGCGTATTTCCTCTCAGATCCGAACAAAATCACCCCAGTAttt                                                                    <  4:68551/148‑1 (MQ=255)
                                                                                     tGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGttt                                                            <  6:228345/147‑1 (MQ=255)
                                                                                         cATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTAttt                                                                    <  7:116116/135‑1 (MQ=255)
                                                                                         cATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTAttt                                                                    >  8:116116/1‑135 (MQ=255)
                                                                                         cATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTAc                                                      >  8:144331/1‑149 (MQ=255)
                                                                                         cATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTAc                                                      >  2:66865/1‑149 (MQ=255)
                                                                                          aTGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTAcc                                                     >  4:199263/1‑149 (MQ=255)
                                                                                                gcATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAg                                               <  7:34157/149‑1 (MQ=255)
                                                                                                 cATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATc                                                                                           >  7:171745/1‑104 (MQ=255)
                                                                                                 cATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATc                                                                                           <  8:171745/104‑1 (MQ=255)
                                                                                                 cATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAgg                                              <  6:206674/149‑1 (MQ=255)
                                                                                                     ccGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGggtg                                           <  8:132641/148‑1 (MQ=255)
                                                                                                      cGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGATTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGggtgg                                          <  4:4154/148‑1 (MQ=255)
                                                                                                      cGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGATTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGggtgg                                          <  8:50983/148‑1 (MQ=255)
                                                                                                             tATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCATTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGct                                  <  5:72865/149‑1 (MQ=255)
                                                                                                             tATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGct                                  >  8:158142/1‑149 (MQ=255)
                                                                                                                tGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTctg                               >  5:54855/1‑149 (MQ=255)
                                                                                                                  ttCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTctgct                             <  3:189524/149‑1 (MQ=255)
                                                                                                                      tGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATa                         <  7:144331/149‑1 (MQ=255)
                                                                                                                       gCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATAc                        >  7:26446/1‑149 (MQ=255)
                                                                                                                         aCGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATATCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACcg                      >  6:210276/1‑149 (MQ=255)
                                                                                                                          cgcgGATCAGCCGCTCAAGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACcgt                     <  8:183521/149‑1 (MQ=255)
                                                                                                                              gATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACcgtgcgt                 <  8:26446/149‑1 (MQ=255)
                                                                                                                               aTCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACcgtgcgt                 >  1:153277/1‑148 (MQ=255)
                                                                                                                                tCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGTCTCTGCTGATACCGTGCGTGa               <  5:62586/149‑1 (MQ=255)
                                                                                                                                tCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGa               <  6:54855/149‑1 (MQ=255)
                                                                                                                                  aGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGata             >  1:59845/1‑149 (MQ=255)
                                                                                                                                  aGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGat              <  6:83876/148‑1 (MQ=255)
                                                                                                                                  aGCCGCTCACGGTTAGTTCCCGCCATATAAACGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGat              <  6:129169/148‑1 (MQ=255)
                                                                                                                                   gCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGatat            <  7:122721/149‑1 (MQ=255)
                                                                                                                                          aCGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGc     >  8:42675/1‑149 (MQ=255)
                                                                                                                                             gTTATTTCCAGCCATATAAAAGCTTAAGCGTTATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCttt  <  5:80047/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTT  >  NC_000913/1814132‑1814421

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑