Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,017,566 C→T 6.9% P338L (CCG→CTG)  fliI → flagellum‑specific ATP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,017,5660CT 6.9% 182.0 / 3.4 58P338L (CCG→CTG) fliIflagellum‑specific ATP synthase
Reads supporting (aligned to +/- strand):  ref base C (25/29);  new base T (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGTGAGCAACATTACGCGC  >  NC_000913/2017425‑2017712
                                                                                                                                             |                                                                                                                                                  
cccccGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGTAGGCGATGACCAGCAGGATCCGATTGcc                                                                                                                                             <  7:62833/149‑1 (MQ=255)
   ccGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGAc                                                                                                                                          >  1:179938/1‑149 (MQ=255)
          caAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCcgcgc                                                                                                                                   <  2:179938/149‑1 (MQ=255)
            aaGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCgg                                                                                                                                 <  4:331040/149‑1 (MQ=255)
                  aTCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATc                                                                                                                           >  3:315657/1‑149 (MQ=255)
                      aCCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGAGCTCACTGAAGGCGATGACCAGCAGGATCTGATTGCCGACTCCGCGCGGGCGATCCTCg                                                                                                                       >  3:212386/1‑149 (MQ=255)
                      aCCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGAGCTCACTGAAGGCGATGACCAGCAGGATCTGATTGCCGACTCCGCGCGGGCGATCCTCg                                                                                                                       >  3:217291/1‑149 (MQ=255)
                         gTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGAGGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACg                                                                                                                    <  6:64254/149‑1 (MQ=255)
                         ctaggggtTTGCAAAATTACCGGCAATGGTCGCGAGTGCCGGCAATGTCATTAGCGGCGGCGGCTCGAGCACAGCTTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCTGGCGAGCCTCGACg                                                                                                                    <  8:327026/143‑1 (MQ=255)
                           cGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGt                                                                                                                  <  4:97080/149‑1 (MQ=255)
                              tgtTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTcac                                                                                                               <  3:136697/149‑1 (MQ=255)
                               gtTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTcaca                                                                                                              <  8:357122/149‑1 (MQ=255)
                                     cAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGc                                                                                                        >  1:38364/1‑149 (MQ=255)
                                     cAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGc                                                                                                        >  5:149723/1‑149 (MQ=255)
                                         ttACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGCCCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGtc                                                                                                    <  2:43546/149‑1 (MQ=255)
                                         ttACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGAGCTCACTGAAGGCGATGACCAGCAGGATCTGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGtc                                                                                                    <  4:217291/149‑1 (MQ=255)
                                         ttACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGAGCTCACTGAAGGCGATGACCAGCAGGATCTGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGtc                                                                                                    <  4:212386/149‑1 (MQ=255)
                                             cGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGc                                                                                                >  3:312692/1‑149 (MQ=255)
                                                  cTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGAGGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACt                                                                                           <  5:93220/149‑1 (MQ=255)
                                                     gTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACTGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTgg                                                                                         <  3:193623/148‑1 (MQ=255)
                                                        gAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGa                                                                                     <  8:304852/149‑1 (MQ=255)
                                                         aGCGTGGCGGAAATGGCATTTGCGGCGGCGGCTCGATTAGCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGaa                                                                                    <  4:290688/149‑1 (MQ=255)
                                                         aGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGaa                                                                                    <  4:98215/149‑1 (MQ=255)
                                                              gCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCgg                                                                               >  7:399680/1‑149 (MQ=255)
                                                                 ggAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCa                                                                            >  4:363930/1‑149 (MQ=255)
                                                                       ggCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATcc                                                                      >  5:277291/1‑149 (MQ=255)
                                                                       ggCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGCCGATGACCAGCAGGATCCGATTGCCGACACAGCGCGTCCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATcc                                                                      >  6:63146/1‑149 (MQ=255)
                                                                                   ggcggcTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGatat                                                          <  3:113946/149‑1 (MQ=255)
                                                                                    gcggcTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATAtt                                                         <  3:366866/149‑1 (MQ=255)
                                                                                       gcTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGaa                                                      <  6:277291/149‑1 (MQ=255)
                                                                                           gATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGt                                                  >  4:131538/1‑149 (MQ=255)
                                                                                            aTTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTc                                                 >  7:4363/1‑149 (MQ=255)
                                                                                                ccGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGaa                                                      >  2:290415/1‑140 (MQ=255)
                                                                                                ccGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGaa                                                      <  1:290415/140‑1 (MQ=255)
                                                                                                  gcgTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGc                                                                             >  4:277366/1‑115 (MQ=255)
                                                                                                  gcgTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGc                                                                             <  3:277366/115‑1 (MQ=255)
                                                                                                    gTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGcc                                         <  3:72628/149‑1 (MQ=255)
                                                                                                          taCCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCaa                                   >  5:219637/1‑149 (MQ=255)
                                                                                                            ccGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAAt                                  <  6:14318/148‑1 (MQ=255)
                                                                                                               tGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGAc                               >  4:242098/1‑148 (MQ=255)
                                                                                                               tGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGAc                               <  3:242098/148‑1 (MQ=255)
                                                                                                                 cTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACgg                             <  8:4363/148‑1 (MQ=255)
                                                                                                                    aCTGCAGGCGATGACCAGCAGGATCCGAGTGCCGACTCCGCGCGGGCGATCCTCGACGGGCTCATTGTGCTGTCTGGCCGACTGGCGGAAGCCGGGCACTATCCGGCTTTCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGtt                         <  8:360782/149‑1 (MQ=255)
                                                                                                                    aCTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCggt                         >  5:25316/1‑147 (MQ=255)
                                                                                                                        aaGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGCTCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGAt                      >  1:187014/1‑148 (MQ=255)
                                                                                                                         aGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGACGCGTCGATCAGCCGCGCAATGACGGCGTTGAt                      >  6:384045/1‑147 (MQ=255)
                                                                                                                           gCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGt                  <  6:149723/149‑1 (MQ=255)
                                                                                                                             gATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGTGa                <  8:399680/149‑1 (MQ=255)
                                                                                                                                gACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGTGAGc              >  1:361129/1‑148 (MQ=255)
                                                                                                                                   cagcagGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGTGAGCAACa          >  3:327875/1‑149 (MQ=255)
                                                                                                                                    agcagGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGTGAGCAACAt         >  8:290598/1‑149 (MQ=255)
                                                                                                                                     gcagGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGa                                               >  3:368260/1‑110 (MQ=255)
                                                                                                                                     gcagGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGa                                               <  4:368260/110‑1 (MQ=255)
                                                                                                                                     gcagGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGTGAGCAACAt         >  3:207886/1‑148 (MQ=255)
                                                                                                                                      cagGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGTGAGCAACATTa       >  8:227243/1‑149 (MQ=255)
                                                                                                                                         gATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGTGAGCAACATTAcgc    <  2:38364/149‑1 (MQ=255)
                                                                                                                                           tCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGt                  <  5:130941/133‑1 (MQ=255)
                                                                                                                                           tCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGt                  >  6:130941/1‑133 (MQ=255)
                                                                                                                                            ccGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGTGAGCAACATTAcgcgc  >  6:322023/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CCCCCGCCACCAAAGGTTATCCACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCGCAATGACGGCGTTGATCAGTGAGCAACATTACGCGC  >  NC_000913/2017425‑2017712

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