Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 877,298 G→A 7.0% D246D (GAC→GAT rimO ← ribosomal protein S12 methylthiotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913877,2980GA 7.0% 182.8 / 3.1 57D246D (GAC→GATrimOribosomal protein S12 methylthiotransferase
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base A (2/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCGCGCCACTGTTTGATGCGCGCCAGTTGGCGATCTACAGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCC  >  NC_000913/877155‑877438
                                                                                                                                               |                                                                                                                                            
aTTTCGCGCCACTGTTTGATGCGCGCCAGTTGGCGATCTACAGTACCCGGAGGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGAcgtcgt                                                                                                                                          <  3:163907/148‑1 (MQ=255)
 tttCGCGCCACTGTTTGATGCGCGCCAGTTGGCGAGCTACAGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACgtcgtc                                                                                                                                         <  5:45259/148‑1 (MQ=255)
  ttCGCGCCACTGTTTGATGCGCGCCAGTTGGCGATCTACAGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGACTTTGCCTTCTGCCATCAGTGGGATGACgtcgtc                                                                                                                                         >  3:287752/1‑147 (MQ=255)
      cgcCACTGTTTGATGCGCGCCAGTTGGCGATCTACAGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACAt                                                                                                                                    >  7:233675/1‑148 (MQ=255)
              tttGATGCGCGCCAGTTGGCGATCTACAGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGACCACATGCGGAt                                                                                                                              >  2:59015/1‑146 (MQ=255)
               ttGATGCGCGCCAGTTGGCGATCTACAGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAg                                                                                                                           >  6:89662/1‑148 (MQ=255)
                  aTGCGCGCCAGTTGGCGATCTACAGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAggg                                                                                                                         >  6:31684/1‑147 (MQ=255)
                   tGCGCGCCAGTTGGCGATCTACAGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGt                                                                                                                        >  7:90274/1‑147 (MQ=255)
                             ttGGCGATCTACAGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGc                                                                                                             <  1:59015/148‑1 (MQ=255)
                              tGGCGATCTACAGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTCCCTTCTGCCATCAGTGGGATGACGTCGTCAACATGCGGATAAGGGTAAACTTAGAGCa                                                                                                            >  3:131348/1‑148 (MQ=255)
                                        cAGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTcc                                                                                                   <  8:90274/147‑1 (MQ=255)
                                         aGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGACTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAg                                                                                                 <  4:287752/148‑1 (MQ=255)
                                                  aCGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGt                                                                                        >  7:299738/1‑148 (MQ=255)
                                                   cGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGtt                                                                                       <  5:89662/148‑1 (MQ=255)
                                                       tcatcaGTTTTAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGtt                                                                                       <  7:200558/144‑1 (MQ=255)
                                                       tcatcaGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGtt                                                                                       >  8:200558/1‑144 (MQ=255)
                                                        catcaGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGa                                                                                   <  5:31684/147‑1 (MQ=255)
                                                         atcaGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGAt                                                                                  >  8:163371/1‑147 (MQ=255)
                                                         atcaGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATa                                                                                 <  3:173017/148‑1 (MQ=255)
                                                               ttGAGAATGCGCGGGCAGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGAGATGACGTCGTCCACAAGAGGATAAGCGTAAACGTAGAGCAAACTTGTCAACATCCCAAGTTTAGATAACTGt                                                                            >  2:251164/1‑147 (MQ=255)
                                                                tGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATa                                                                                 >  8:100765/1‑141 (MQ=255)
                                                                tGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATa                                                                                 <  7:100765/141‑1 (MQ=255)
                                                                tGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGtt                                                                           <  7:163371/147‑1 (MQ=255)
                                                                 gagaATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTc                                                                          <  3:242624/147‑1 (MQ=255)
                                                                    aaTGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGca                                                                       >  2:173898/1‑147 (MQ=255)
                                                                      tGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAg                                                                    <  2:119064/148‑1 (MQ=255)
                                                                       gcgcGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAgg                                                                   >  8:334902/1‑148 (MQ=255)
                                                                         gcgGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCt                                                                 >  2:80291/1‑148 (MQ=255)
                                                                             gCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGAcc                                                             <  7:334902/148‑1 (MQ=255)
                                                                                 gCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATg                                                          <  2:314452/147‑1 (MQ=255)
                                                                                          aaCGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTa                                                >  8:12260/1‑148 (MQ=255)
                                                                                             ggAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCg                                                                                    >  8:328276/1‑109 (MQ=255)
                                                                                             ggAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCg                                                                                    <  7:328276/109‑1 (MQ=255)
                                                                                              gAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGCTCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGGTTTTACCgg                                            >  8:103049/1‑148 (MQ=255)
                                                                                               aaTGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGcac                                                                      >  4:74328/1‑121 (MQ=255)
                                                                                               aaTGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGcac                                                                      <  3:74328/121‑1 (MQ=255)
                                                                                                aTGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGtt                                                                           >  4:193237/1‑115 (MQ=255)
                                                                                                aTGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGtt                                                                           <  3:193237/115‑1 (MQ=255)
                                                                                                aTGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCt                                          >  1:189512/1‑148 (MQ=255)
                                                                                                aTGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCt                                          <  8:233675/148‑1 (MQ=255)
                                                                                                  gTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCg                                        >  5:147880/1‑148 (MQ=255)
                                                                                                  gTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTATTTACCGGCTCg                                        >  8:105169/1‑148 (MQ=255)
                                                                                                    ccAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGAcc                                                             >  6:100692/1‑125 (MQ=255)
                                                                                                    ccAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGAcc                                                             <  5:100692/125‑1 (MQ=255)
                                                                                                    ccAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGcc                                      >  4:170623/1‑148 (MQ=255)
                                                                                                      aGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGt                                    <  1:80291/148‑1 (MQ=255)
                                                                                                       gATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACATCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGta                                   >  8:219230/1‑147 (MQ=255)
                                                                                                       gATACCGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACATCGTCCACATGCGGATAAGCGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGATCGCACAGGCTCACCAAGCTCGTTTTTACCGGCTCGCCg                                     >  8:224482/1‑146 (MQ=255)
                                                                                                              cAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGTGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAg                            <  2:189512/148‑1 (MQ=255)
                                                                                                                ggATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGcc                          <  7:12260/148‑1 (MQ=255)
                                                                                                                 gATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCa                         <  8:299738/148‑1 (MQ=255)
                                                                                                                        ccTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGa                   <  7:105169/147‑1 (MQ=255)
                                                                                                                           tCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTgg                               <  5:201276/132‑1 (MQ=255)
                                                                                                                           tCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGtt               >  1:181464/1‑148 (MQ=255)
                                                                                                                           tCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTAACCGGCTCGCCGTTGTgg                               >  6:201276/1‑132 (MQ=255)
                                                                                                                             tGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGttt              >  3:146724/1‑147 (MQ=255)
                                                                                                                                    aGTGGGATGACATCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCa       <  7:224482/147‑1 (MQ=255)
                                                                                                                                    aGTGGGATGACATCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCa       <  7:219230/147‑1 (MQ=255)
                                                                                                                                        ggATTACGTCGTCCACATGCGGATAAGGGTAAAAGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAATCCAGTACGATGTTTAACATCAACGcc  <  1:251164/148‑1 (MQ=255)
                                                                                                                                        ggATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGcc  <  3:170623/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
ATTTCGCGCCACTGTTTGATGCGCGCCAGTTGGCGATCTACAGAACCCGGACGCTTCATCAGTTTGAGAATGCGCGGGCTGGCGTGCTGCAACGGAATGTCCAGATACGGCAGGATTTTGCCTTCTGCCATCAGTGGGATGACGTCGTCCACATGCGGATAAGGGTAAACGTAGTGCAGACGTGTCCAGATCCCCAGTTTCGATAACTGTTCGCACAGGCTGACCATGCTGGTTTTTACCGGCTCGCCGTTGTGGAAGCCAGTACGATGTTTAACATCAACGCC  >  NC_000913/877155‑877438

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: