Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,362,216 C→G 12.6% Q48H (CAG→CAC tpiA ← triose phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,362,2160CG12.6% 148.5 / 11.1 55Q48H (CAG→CACtpiAtriose phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base G (4/3);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

GCTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATGCTTT  >  NC_002947/5362068‑5362364
                                                                                                                                                     |                                                                                                                                                    
gCTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAcc                                                                                                                                                      >  5:153465/1‑149 (MQ=255)
 cTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCt                                                                                                                                                     <  8:47033/149‑1 (MQ=255)
        tGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCtgattgat                                                                                                                                              <  6:153465/149‑1 (MQ=255)
          gCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatgg                                                                                                                                            <  4:81780/149‑2 (MQ=255)
            cAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGaa                                                                                                                                           >  2:113478/1‑148 (MQ=255)
             aaTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAAc                                                                                                                                          >  6:25349/1‑148 (MQ=255)
              aTCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAAc                                                                                                                                          >  6:198797/1‑147 (MQ=255)
               tCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAg                                                                                                                                       <  5:25349/149‑1 (MQ=255)
                  acacCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGcc                                                                                                                                    >  3:281495/1‑149 (MQ=255)
                    acCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCgg                                                                                                                                  <  5:245175/149‑1 (MQ=255)
                            gCAACCTGCTTCAAACAGCTTCTAGGGAGCAACTTCCCCGGTCAGCGCACACTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACATGAGTGATGAAAAAGGCCGGTGGAAATa                                                                                                                          >  2:45516/1‑149 (MQ=255)
                             cAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATAc                                                                                                                         >  7:12086/1‑149 (MQ=255)
                              aaCCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTTTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACg                                                                                                                        <  5:198797/149‑1 (MQ=255)
                                ccTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGc                                                                                                                      <  4:281495/149‑1 (MQ=255)
                                      tCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATtgatga                                                                                                                                            >  1:120559/1‑121 (MQ=255)
                                      tCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATtgatga                                                                                                                                            <  2:120559/121‑1 (MQ=255)
                                        aaCCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatg                                                                                                                                             <  6:145985/118‑1 (MQ=255)
                                        aaCCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCAGGTGATtgatg                                                                                                                                             >  5:145985/1‑118 (MQ=255)
                                        aaCCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCt                                                                                                                                                     <  3:286847/110‑1 (MQ=255)
                                        aaCCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCt                                                                                                                                                     >  4:286847/1‑110 (MQ=255)
                                                 ctCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCg                                                                                                     >  4:151548/1‑149 (MQ=255)
                                                    gggAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGc                                                                                                   <  5:27930/148‑1 (MQ=255)
                                                     ggAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCaa                                                                                                 <  6:186669/149‑1 (MQ=255)
                                                        gCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGc                                                                                              >  1:298030/1‑149 (MQ=255)
                                                         cAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGc                                                                                              >  6:175806/1‑148 (MQ=255)
                                                         cAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGc                                                                                              >  6:289285/1‑148 (MQ=255)
                                                            cTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGAAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTTCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATAAGGCGACTTCCAATCCGCTCGGCTAGGc                                                                                              >  8:286020/1‑145 (MQ=255)
                                                               ccccGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTa                                                                                       >  7:247328/1‑149 (MQ=255)
                                                                          gcACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCtt                                                                             <  3:151548/148‑1 (MQ=255)
                                                                           cACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg                                                                           >  4:313037/1‑149 (MQ=255)
                                                                           cACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATATCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTg                                                                           >  1:233085/1‑149 (MQ=255)
                                                                             cccTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGGACAGTAATTAATTTACATGCCAGGCCATCAATCACGTGATTGATGAACAAGGCCGTTGGAAATACGGCGAGTACCAGTCCGATCGGCAAGGCGAGATTAATCAAGCGTTTgg                                                                          >  1:214774/1‑149 (MQ=255)
                                                                             cccTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGt                                                                         <  1:113478/149‑1 (MQ=255)
                                                                              ccTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGt                                                                         >  8:189155/1‑148 (MQ=255)
                                                                                  ttCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGtcagct                                                                    >  3:85998/1‑149 (MQ=255)
                                                                                      ggTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTCCCAGGCCAACAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGtcagctcagc                                                                >  6:230995/1‑149 (MQ=255)
                                                                                        ttGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGa                                                              >  6:73441/1‑149 (MQ=255)
                                                                                           tACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACg                                                            <  5:289285/148‑1 (MQ=255)
                                                                                           tACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACg                                                            <  5:175806/148‑1 (MQ=255)
                                                                                                  gAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGcgcg                                                    <  2:214774/149‑1 (MQ=255)
                                                                                                        tgtgCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTAcc                                              >  3:294120/1‑149 (MQ=255)
                                                                                                          tgCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTAc                                               <  8:12086/146‑1 (MQ=255)
                                                                                                          tgCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGt                                            >  7:97372/1‑149 (MQ=255)
                                                                                                             gcCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCa                                         <  5:73441/149‑1 (MQ=255)
                                                                                                                 gACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCAttt                                      >  6:80241/1‑149 (MQ=255)
                                                                                                                     gTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCa                                  <  8:247328/149‑1 (MQ=255)
                                                                                                                        aTTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGt                                >  4:161466/1‑148 (MQ=255)
                                                                                                                               tACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAg                          <  7:189155/147‑1 (MQ=255)
                                                                                                                                 ccTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTa                       >  5:85765/1‑148 (MQ=255)
                                                                                                                                         ttcATCAATCAGGTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTc                                                                                    >  3:185021/3‑79 (MQ=39)
                                                                                                                                         ttcATCAATCAGGTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTc                                                                                    <  4:185021/77‑1 (MQ=39)
                                                                                                                                         gCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGc              >  8:263461/1‑149 (MQ=255)
                                                                                                                                            atcaatcaGGTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGAc           >  7:250428/1‑149 (MQ=255)
                                                                                                                                            atcaatcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGAc           <  5:80241/149‑1 (MQ=255)
                                                                                                                                                atcaCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCAt       >  7:255805/1‑149 (MQ=255)
                                                                                                                                                   aCCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATGc     <  2:298030/148‑1 (MQ=255)
                                                                                                                                                     cTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTgg                                                                          <  1:137577/77‑1 (MQ=255)
                                                                                                                                                     cTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTgg                                                                          >  2:137577/1‑77 (MQ=255)
                                                                                                                                                     cTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATGCttt  <  7:82960/149‑1 (MQ=255)
                                                                                                                                                     cTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATGCtt   >  1:23750/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTCCGAATGGCCAATCAACACCAACTTGCAACCTGCTTCAACCAGCTGCTCGGGAGCAACTTCCCCGGTCAGCGCACCCTGTTCGGGTTGTACAGCAGAATTCTGTGCGCCGACAGTAATTTCTTTACCTGCCAGGCCATCAATCACCTGATTGATGAACAAGGCCGGTGGAAATACGGCGACTTCCACTCCGCTCGGCAAGGCGAGATTACTCAAGCCTTTGGTCAGCTCAGCGACGCTAGCGCGGGTACCGTGCATTTTCCAGTTACCAGCTACCATAGGGCGACGCATGCTTT  >  NC_002947/5362068‑5362364

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