Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 681,694 A→G 9.0% intergenic (+31/+29) djlC → / ← hscC co‑chaperone DjlC/chaperone protein HscC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913681,6940AG 9.0% 127.3 / 3.5 57intergenic (+31/+29)djlC/hscCco‑chaperone DjlC/chaperone protein HscC
Reads supporting (aligned to +/- strand):  ref base A (27/23);  new base G (4/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

AAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCC  >  NC_000913/681554‑681830
                                                                                                                                            |                                                                                                                                        
aaaaTGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGc                                                                                                                                        <  1:430766/143‑1 (MQ=255)
 aaaTGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATgcggcg                                                                                                                                       >  1:123353/1‑143 (MQ=255)
        gaGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTAGGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGAGGTGGTTCTTGCGACGATGCCGGATGCCGCGGTTAcgg                                                                                                                                >  1:1013938/1‑142 (MQ=255)
        gaGCGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgc                                                                                                                                <  1:706069/143‑1 (MQ=255)
         ggggTGAAAAGCCACAGGGCATTACCAGCATTATCAACATCAGTTTAACGTTGCTGGTGATAGTTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgcc                                                                                                                               <  2:1034850/142‑1 (MQ=255)
                aagaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCg                                                                                                                        <  1:339523/140‑1 (MQ=255)
                 aaaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCgg                                                                                                                       >  2:13021/1‑143 (MQ=255)
                   tccgcaGAGCTTTCCTAGCATCACCCAAGTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAACAGCGCGATTTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCa                                                                                                                     <  2:38695/139‑1 (MQ=255)
                        agagCATTACCGGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTGTCGACAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGCTGCAGTGTTTACGCCGCATCCGGGAGCCgt                                                                                                                >  2:614195/1‑143 (MQ=255)
                              tgaCCCGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTGATGCTCTCGGCAAAATCGCGATGTTGGTTAGCGGACTAGGTGTTGTTTGCTCCGATTCCGGAGGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTAt                                                                                                          <  1:840530/141‑1 (MQ=255)
                                 ccAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATggg                                                                                                       >  1:840922/1‑143 (MQ=255)
                                 ccAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGAGTTTATggg                                                                                                       >  2:306731/1‑143 (MQ=255)
                                  cAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCAACCGCAATCCGTGTTTATGGGa                                                                                                      >  1:146882/1‑143 (MQ=255)
                                     cATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGCGTTTATGGGAtct                                                                                                   >  2:98061/1‑143 (MQ=255)
                                        tATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGAGTACTATGCGTCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCaa                                                                                                <  1:303792/143‑1 (MQ=255)
                                             aaaTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGCCGATGTGTTTTTTGCGCCGATCCCGGATGCGGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCACTGGCt                                                                                           >  1:578097/1‑143 (MQ=255)
                                                    tttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCGGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaataat                                                                                    <  1:991545/143‑1 (MQ=255)
                                                    tttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCGGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaataat                                                                                    >  2:766681/1‑143 (MQ=255)
                                                    tttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTCGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGCGTTTATGGGACCTCAATGGCTCaataat                                                                                    >  2:456544/1‑143 (MQ=255)
                                                     ttAACGCTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGGTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtc                                                                                   <  1:50034/143‑1 (MQ=255)
                                                       aaCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTATGCCGCATCCGCCAGCCGGGGTCATGGGATCTCTATTGCTAAATAAtcgc                                                                                 >  2:290916/1‑143 (MQ=255)
                                                        aCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTTTTTTGTTTCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTAGGGGATCTCAATGCCTAAATAAtcgct                                                                                >  2:929811/1‑143 (MQ=255)
                                                          ggTGCTGTTGATAGGTTATGCGCTCGGCAAAAGCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctcg                                                                              <  2:299798/143‑1 (MQ=255)
                                                             gggtgTGATGGGTTATGCGTTCGGCAAAATCGCGAGGTGGACGAGCTGACTATGTGTTTTGTGCGCCGATGGCGGATGCAGCGTTTGCGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgccc                                                                           <  2:779630/139‑1 (MQ=255)
                                                             gCTGGTGATAGGTTATGCTCTCTGCAAAATCGCGATGTTGTTTATCTGACTATGTTTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCAGCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgccc                                                                           <  1:185413/143‑1 (MQ=255)
                                                              cTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTCTGGCACCTCACCGGCTTATTACTAGctcgctcg                                                                          >  1:744190/1‑143 (MQ=255)
                                                              cTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGACGTGATAATCGCGCTGATGCCCGTAGCAGCGGTTACGCCGCATCCGGCAGCCGTGTTTATGGCATCTCAACGGCTAAACAATCgcccgcccc                                                                          >  2:492912/1‑142 (MQ=255)
                                                               tGGTGATAGGTTATGCTCTCGGCACAATCGCGATGTGGTTTAGCTGACGATGTGGTGTCTGCGCCGATGCCGGATGCAGCGTTCACGCCGCCTCCGGCAGCCGTGTTTCTGGGATCTCACTGGCGAAATCATCGCTcgcccgc                                                                         >  1:137488/1‑143 (MQ=255)
                                                               tGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccgc                                                                         >  1:964807/1‑143 (MQ=255)
                                                               tGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTGTTTAGACGATGCTGGATGCAGCTGTGACGCCGCATCCGGCAGGCGTGTGTCTGGGATCTCAATGGCGAAATAAGCGTCcgcccgg                                                                         >  2:142267/1‑142 (MQ=255)
                                                                  tGATAGGTTATGCGCTCGGCAAAATCGGGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCccgccgc                                                                      <  1:231061/143‑1 (MQ=255)
                                                                   gATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGTCTAAATAATCGCTCGCccgccgcc                                                                     >  1:855863/1‑143 (MQ=255)
                                                                    aTAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCAGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCccgccgccg                                                                    >  1:630951/1‑143 (MQ=255)
                                                                       ggTTATGCTCTCGGCTAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGGTGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAAc                                                                 <  2:958802/143‑1 (MQ=255)
                                                                         ttATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTGTTTGCGCCGATGCCGGATGCAGCGGGTTCGCCGCATCCGGCAGCCGTGTTTATGGGGCCGCAACTGCGAACTAGTCGCTCGCCCGCCGCCGAACtt                                                               >  2:282726/1‑143 (MQ=255)
                                                                           aTGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGATTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATCGCTAAGTAATCGCTCGCCCGCCGCCGAACttaa                                                             >  2:901845/1‑141 (MQ=255)
                                                                             gCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGACCGCCGCCGAACTTCAtg                                                           >  2:352599/1‑142 (MQ=255)
                                                                                     cAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCACCGCggg                                                   >  1:365365/1‑143 (MQ=255)
                                                                                       aaaTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAACAAACACACTTTTGTAGCGCTAGCCCGCCGCCGAACTTCATCCACGCGGGcc                                                 >  1:1093208/1‑143 (MQ=255)
                                                                                            gcgATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCAGAACTTCATTCACGAGGGCCAGATc                                            >  2:912847/1‑143 (MQ=255)
                                                                                              gATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCCTGTTTATGGGATCTCACTGGCTAAATAATCGCTCGCCCGCCGCCGACCTTCCTCCACCCGGGccccctccg                                          >  2:196874/1‑136 (MQ=255)
                                                                                               aTGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGtt                                         <  2:600392/143‑1 (MQ=255)
                                                                                                  ttgttTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGa                                      <  2:855388/143‑1 (MQ=255)
                                                                                                   tgttGAGCGGACGATGTGTTTTTTGCGCCGATGCCGGAGGCAGCGCTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGAc                                     <  1:925120/143‑1 (MQ=255)
                                                                                                     tttAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTg                                   <  2:851757/143‑1 (MQ=255)
                                                                                                               gaggtgtTGTGTGCGCCGATGGTGGATGGAGCGTTTACGCCGCATGCGGCAGCCGTGTTTATGGGATCTGAATGGCTAAATAATGGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCa                         <  1:948003/140‑1 (MQ=255)
                                                                                                                  gtgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCCTATCGTTTGACTGCACCGCCGCATCa                      >  1:495805/1‑143 (MQ=255)
                                                                                                                  gtgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCAACGCGGGCCATAAAGTTTGACTACAACGCCGCATCa                      >  1:159605/1‑143 (MQ=255)
                                                                                                                      ttttgGCGCCGATGCCGGATGCGGCGTTGACGCCGCCGCAGGCAGCTGTGTTTATTGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAAt                  <  2:230896/143‑1 (MQ=255)
                                                                                                                      ttttgGCGCCGATGCCGGAGGCAGCGTTTCCGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAAt                  <  1:349502/143‑1 (MQ=255)
                                                                                                                      gttgtGCGCGGAGGCCGCATGCAGCGTTGAGGCCGCATCGGGCAGCAGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAAt                  <  2:159912/139‑1 (MQ=255)
                                                                                                                        tctGCCCCGATGCCGGATGCAGCTTTTACGCCCCACCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATCATCGCTCGCCCGCCGCCCAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCg                <  2:375396/141‑1 (MQ=255)
                                                                                                                          cgcgcCGATGTCGGATGCAGCGTTTACGCCGCAGCCGGCAGCCGGGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGt              <  2:677207/142‑1 (MQ=255)
                                                                                                                             gcCGATGCCGGATGCAGCGTTTACGCCGCCTCCGGCAGCCGTGTGTATGGGATCTCAATGGCGAAATGATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCCTATCGTTTCGCCGCAACCCCACCTCAAAACCGGTGAt           >  2:244287/1‑143 (MQ=255)
                                                                                                                                  tGCCGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGCCTGCAACGCCGCACCAAATCGGTGATCACCCg     >  2:799810/1‑143 (MQ=255)
                                                                                                                                   gCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCg     <  1:1009389/143‑1 (MQ=255)
                                                                                                                                     cGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCCCCCGGc   >  2:347402/1‑143 (MQ=255)
                                                                                                                                      ggATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGCACTTCATCCCCGCGGGCCATATCGTTTGACTGCAACACAACATAAAAATCGGTGATGACCCGGcc  >  1:1038272/1‑143 (MQ=255)
                                                                                                                                            aGCGTTTACGCCGCATCCGGCAGCCGTGTCTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAAcgc                             >  2:697517/1‑110 (MQ=255)
                                                                                                                                            aGCGTTTACGCCGCATCCGGCAGCCGTGTCTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAAcgc                             <  1:697517/110‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                        
AAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCC  >  NC_000913/681554‑681830

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: