Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 83,590 G→T 7.6% intergenic (‑61/+32) leuA ← / ← leuL 2‑isopropylmalate synthase/leu operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091383,5900GT 7.6% 126.3 / 2.7 56intergenic (‑61/+32)leuA/leuL2‑isopropylmalate synthase/leu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (31/18);  new base T (2/2);  total (35/21)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGC  >  NC_000913/83448‑83730
                                                                                                                                              |                                                                                                                                            
tttctttcACACTCTAGCTTGGCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATg                                                                                                                                              <  1:382660/143‑1 (MQ=255)
       cacaCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTCTAGAGCGCCTCGCTTCGGGCCTAAAAAACCCCGCGCAACGCCCCggg                                                                                                                                       >  2:340457/1‑143 (MQ=255)
          acttacgCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGttt                                                                                                                                    <  1:645039/137‑1 (MQ=255)
          acTCGAGCTCGCCTGCAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGTGCATAAAAAAACCCGCGCCATTGCGCGGGttt                                                                                                                                    <  2:15889/143‑1 (MQ=255)
          acTCACGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGttt                                                                                                                                    <  1:497287/143‑1 (MQ=255)
          acTCAAGCTTGCCTGTAGCGCCTGTTCGCCGTCGCGCAGTGTGGTATCGAGAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGGCTCTGTTAGAGCGCCTCGCTTCGGGCATAAAAAAACTCGCGCAATGGCGCGGGttt                                                                                                                                    <  2:251567/143‑1 (MQ=255)
          acTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGttt                                                                                                                                    >  1:782674/1‑143 (MQ=255)
           cTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGtttt                                                                                                                                   >  2:337876/1‑143 (MQ=255)
           cTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGtttt                                                                                                                                   <  1:39546/143‑1 (MQ=255)
           cTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAAGAGCGCgggctgt                                                                                                                                   >  2:690427/1‑139 (MQ=255)
            tccggCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAACTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGGttttt                                                                                                                                  <  2:23885/139‑1 (MQ=255)
            tCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGgtttt                                                                                                                                  >  1:713655/1‑143 (MQ=255)
                       tGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATTAAAAAACCCGCGCAATGGCGCCGGGTTTTTGGTTGACTGc                                                                                                                       >  2:493798/1‑143 (MQ=255)
                             gCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTg                                                                                                                 >  1:474989/1‑143 (MQ=255)
                             cccTGTTCACCGTCGCGCAATGCGGTACGGAAAATAATGACGTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTg                                                                                                                 <  1:632766/142‑1 (MQ=255)
                                   tCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTaa                                                                                                           <  1:260998/143‑1 (MQ=255)
                                   tCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCCTac                                                                                                           >  2:20080/1‑142 (MQ=255)
                                     aCCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATg                                                                                                         >  1:475158/1‑143 (MQ=255)
                                     aCCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATg                                                                                                         <  2:214058/143‑1 (MQ=255)
                                        gTCGCGTAATGTTGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTTTCTCTCATAGAGTGCCTCCCTTCGGTCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTg                                                                                                      <  1:376204/143‑1 (MQ=255)
                                           gcgcAATGTGGTATCGAACATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGAt                                                                                                   <  2:558400/143‑1 (MQ=255)
                                                   tgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTGGGGCATAAAAAAACCCGCGCAATGGCGGCGCTTTTTTATTTGAGTGCCTGCTGGCTTACTGCCGGATGCCGCtaa                                                                                           >  2:407633/1‑141 (MQ=255)
                                                   tgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGATGGATGCCGCTcc                                                                                           >  2:89363/1‑142 (MQ=255)
                                                   ggggaTCGAAAATAGTGACTTGCCGGCTCATGGTTTGGGTCCTTGTCGCGGTTAGGGCGCCTCGCTTCGTGCATAAAAATCCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCa                                                                                           <  2:239672/139‑1 (MQ=255)
                                                                    aCTTGCTGTCTCATGGTCTGGGTACGTGTCGCTGGTAGAGCGCCTCGCGTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCaaa                                                                          <  1:293430/143‑1 (MQ=255)
                                                                        gCTTGCTCATTGTTTGGTTCCCTGTCTCCTGTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAGTGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaag                                                                      <  2:464446/143‑1 (MQ=255)
                                                                         cTGGCTCATGGTTTGGGTCCTTGTCTCCTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGCGCTGGCTTAATGCTGGATGCCGCTCACACGTCTCACGCGCCacgaaca                                                                     >  2:788765/1‑141 (MQ=255)
                                                                              tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGTCTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAGAGAAGCTGCgc                                                                >  2:541892/1‑142 (MQ=255)
                                                                              tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGTCTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCACAGAAGATGCGt                                                                >  2:580450/1‑143 (MQ=255)
                                                                              tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGCTGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGt                                                                >  1:46945/1‑143 (MQ=255)
                                                                              tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGt                                                                >  2:217515/1‑143 (MQ=255)
                                                                              tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGt                                                                >  1:742641/1‑143 (MQ=255)
                                                                              tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCCCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAGTGCTGGATGCCGCTCACTCGTCTACCGCGCGACGAAGATGCGt                                                                >  2:453062/1‑143 (MQ=255)
                                                                              tCATGGTTTGGGTCCTTGTCTCTCTTAGAGCGCCTCGCTGCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTCGTCTGCGTGCTGGCTTAATGCTGGATGCCGCCCACTCGTCTACCGCGCCAAGGCGATGCGt                                                                >  2:614605/1‑143 (MQ=255)
                                                                                    tttGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagta                                                          <  1:249354/143‑1 (MQ=255)
                                                                                     ttGGGCCATTGGCGCTGTGAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTGCTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtag                                                         <  1:652409/143‑1 (MQ=255)
                                                                                       gggTCCTTGTCTCTGTTAGAGCCCCTCGCTTCGGGCATAAAAATACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagta                                                       <  2:687590/143‑1 (MQ=255)
                                                                                       cgcgccTTGTCCCTTCTACAGCGCCTCGCTTCGGGCACAAAAAAACCCGCGCAATCGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagta                                                       <  1:193500/139‑1 (MQ=255)
                                                                                         gTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtcgtagt                                                     >  2:675180/1‑143 (MQ=255)
                                                                                          tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagta                                                    <  2:782674/143‑1 (MQ=255)
                                                                                              tGTCTCTTCTCGCGCGCTCCGCTTTGGTCATAAAAAATCCCGCGCAATGGCGCGTGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGTCGCTCACTCGTCTACCGCGCATAGAAGATTCGTTTAGTAGTAGTAGAcc                                                <  2:272846/143‑1 (MQ=255)
                                                                                                  tctTTTAGAGCGCCTCGCTACGGCCATAAACAAACCCGCGCAATGGCGCTGGTTTATTGTTTCACTGCGTGCTAGCTTAACGCTGGATGCCGCCCACTCATCTACCGCGCAGAGAAGAAGCGTTCACTAGACGTAGACCGATa                                            >  2:641467/1‑143 (MQ=255)
                                                                                                    ttttAGAGCGCCTCGCTTCGGGCATCACCCAACCCGCTCAATGGCGCGGGTGTTTTGTTTGACTGCGTGCTGGCTTCATGCTGGATGCCGCTCACTCGTCTACCGGGCACAGAAGAGGCGTGTAGTAGTACTAGACCGAtgag                                          >  2:449585/1‑140 (MQ=255)
                                                                                                    ttttAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATaaa                                          >  2:413874/1‑143 (MQ=255)
                                                                                                     tttAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTTTCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGACGATGCGTTTCGTAGTAGTAGACCGATAAAg                                         >  2:615100/1‑143 (MQ=255)
                                                                                                           gcgcCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACg                                   >  2:582392/1‑143 (MQ=255)
                                                                                                           gcgcCTCGCTTCGGGCATAAAAAAACCAGGGCCATCACGCGTGGTGTTTGTTTGACTGCGTGCTCGCACACGGCTGGCGGCCGCTCACTCGGCTACCGCGCAAAGAAGATTCTCCTGGTAGTAGGAGACCGATGAAGCGAACg                                   >  2:15545/1‑143 (MQ=255)
                                                                                                                tCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtgtg                              >  1:46805/1‑143 (MQ=255)
                                                                                                                           aTAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAAt                   >  2:545298/1‑143 (MQ=255)
                                                                                                                             aaaaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCGCAGAAGATGCGCTTAGGCGTCGTAGCCCGCTAAAGCGAACGAGGTGAGTCCTTCACTCa                 >  2:568052/1‑143 (MQ=255)
                                                                                                                             aaaaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCa                 <  1:690427/143‑1 (MQ=255)
                                                                                                                                    cccGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGCTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAg          >  2:238126/1‑143 (MQ=255)
                                                                                                                                     ccGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGa         >  1:521323/1‑143 (MQ=255)
                                                                                                                                     ccGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGCGTGCCGCTCACTCGTCTACCGCGCAAAGAAGCTGCGTTTAGTGGTAGTAGACTGATAACGAGGACGATGTGCGACATTAAGTCAGCTCCCGa         >  2:665936/1‑143 (MQ=255)
                                                                                                                                      cgcgCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGAt        >  1:172228/1‑143 (MQ=255)
                                                                                                                                       gcgcAATGGCGCGGGTTTTTTGTTTGACTGAGTGCTGGCTTAATTCTGGATTCCGCTCACTCGTCTACCGCGCCAAGAAGATCGGTTTAGTAGTAGTAGACCGATAAAGCGCACGATGTGATTCATTAAATCAACTCCAGATg       >  2:660609/1‑143 (MQ=255)
                                                                                                                                         gcAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGaa     >  2:444729/1‑143 (MQ=255)
                                                                                                                                            agggCGCGTGTTTTTTGTTTGACTGCGTGCTGGTTTAAGGCTGGATGCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGc  <  2:224684/141‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
TTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGC  >  NC_000913/83448‑83730

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: ‑