Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,289,857 A→T 8.6% N375I (AAT→ATT)  PP_1126 → hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,289,8570AT 8.6% 173.4 / 7.1 58N375I (AAT→ATT) PP_1126hydrolase
Reads supporting (aligned to +/- strand):  ref base A (32/21);  new base T (3/2);  total (35/23)
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.49e-01

CCCGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACC  >  NC_002947/1289711‑1290003
                                                                                                                                                  |                                                                                                                                                   
cccGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGc                                                                                                                                                  <  2:258423/150‑1 (MQ=255)
          tCCCTGGGCGTGAACGGCGTGGCAGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAACGCATGTTCATGGCCTTGCCCGCTGAGGCTGATGCACTAGCCGCTTTCTACGGAGAAAGCAACAACGACGATGGCAATGCCTCGTCGCaa                                                                                                                                       <  2:272702/150‑1 (MQ=255)
            ccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtcgttgtt                                                                                                                                     <  2:1589663/150‑6 (MQ=255)
            ccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtc                                                                                                                                           <  1:1117206/144‑1 (MQ=255)
            ccTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtc                                                                                                                                           >  2:1117206/1‑144 (MQ=255)
                      aaCGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCtcgtc                                                                                                                                           >  2:1540599/1‑134 (MQ=255)
                      aaCGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCtcgtc                                                                                                                                           <  1:1540599/134‑1 (MQ=255)
                        cGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCaa                                                                                                                         <  1:675011/150‑1 (MQ=255)
                         gCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAg                                                                                                                        >  1:348125/1‑150 (MQ=255)
                             tGGCGGCAGACGGCGTGTGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAgg                                                                                                                    <  1:491439/150‑1 (MQ=255)
                               gcggcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAgggg                                                                                                                  <  1:46221/150‑1 (MQ=255)
                                cggcAGAGGGCGGGCCTACCTATGTCGCCGGGCGCGAAGGCATGTTGCAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAACGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGgc                                                                                                                 <  2:600286/150‑1 (MQ=255)
                                cggcAGACGGCGGTCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGgc                                                                                                                 >  1:1215311/1‑150 (MQ=255)
                                         gcgggcgTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAg                                                                                                                     >  1:630002/1‑137 (MQ=255)
                                         gcgggcgTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAg                                                                                                                     <  2:630002/137‑1 (MQ=255)
                                           gggcgTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGgcgc                                                                                                               <  1:888779/141‑1 (MQ=255)
                                           gggcgTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGgcgc                                                                                                               >  2:888779/1‑141 (MQ=255)
                                                    tATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTAcgcc                                                                                             <  2:1215311/150‑1 (MQ=255)
                                                                   gAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGt                                                                              >  2:1590887/1‑150 (MQ=255)
                                                                         aTGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTAt                                                                        >  2:1145791/1‑150 (MQ=255)
                                                                             tCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCACTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCa                                                                    >  2:1186715/1‑150 (MQ=255)
                                                                               aaGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCaaa                                                                  >  2:1431663/1‑150 (MQ=255)
                                                                                aGGTCTTGCCCGCTGAGGGCGATGCACTGGCCGATATGCCCGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAAc                                                                 <  1:499004/150‑1 (MQ=255)
                                                                                 ggCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAAcc                                                                <  1:1145791/150‑1 (MQ=255)
                                                                                  gCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAAccc                                                               >  1:1522507/1‑150 (MQ=255)
                                                                                    cTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTAAGCCGCCAACCTGCAGCGTACCTATGGCCAATCCTg                                                             >  2:92850/1‑150 (MQ=255)
                                                                                                gggCGATGCACTGACCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGt                                                 >  2:1573174/1‑150 (MQ=255)
                                                                                                   cGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCCGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGa                                              >  1:192926/1‑150 (MQ=255)
                                                                                                    gATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAg                                             >  1:686570/1‑150 (MQ=255)
                                                                                                     aTGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCATGCGCGGCTGCTGGTTCAACTCGGACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAgg                                            <  1:1019734/150‑1 (MQ=255)
                                                                                                           tGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgccagc                                      >  1:173956/1‑150 (MQ=255)
                                                                                                           tGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgccagc                                      >  1:173817/1‑150 (MQ=255)
                                                                                                               cGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCcagccagc                                  >  2:1039001/1‑150 (MQ=255)
                                                                                                               cGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCcagccagc                                  <  1:1431663/150‑1 (MQ=255)
                                                                                                               cGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCcagccagc                                  >  1:1398811/1‑150 (MQ=255)
                                                                                                                  tttCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGa                               >  2:499919/1‑150 (MQ=255)
                                                                                                                        tGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTgcgc                         >  2:202058/1‑150 (MQ=255)
                                                                                                                        tGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGACCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTgcgc                         >  1:377503/1‑150 (MQ=255)
                                                                                                                            gAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTc                                                                                                         >  2:980831/1‑66 (MQ=255)
                                                                                                                            gAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTc                                                                                                         <  1:980831/66‑1 (MQ=255)
                                                                                                                            gAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGcgccg                     <  1:202058/150‑1 (MQ=255)
                                                                                                                             aaaGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGcgccgc                    >  1:1048133/1‑150 (MQ=255)
                                                                                                                              aaGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCgccgcc                   >  1:787313/1‑150 (MQ=255)
                                                                                                                                gCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGccgccgc                 >  1:865351/1‑149 (MQ=255)
                                                                                                                                  ccttgcgacgaGGCATTGCCATCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTg               >  1:183329/6‑150 (MQ=255)
                                                                                                                                  aacaacGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTg               >  2:1211681/1‑150 (MQ=255)
                                                                                                                                    ttgcgacgaGGCATTGCCATCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCt             >  2:1492516/4‑150 (MQ=255)
                                                                                                                                    caacGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCt             >  2:17854/1‑150 (MQ=255)
                                                                                                                                     aacGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGc                                                                                   >  2:1322352/1‑79 (MQ=255)
                                                                                                                                     aacGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGc                                                                                   <  1:1322352/79‑1 (MQ=255)
                                                                                                                                     aacGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTg            <  2:1398811/150‑1 (MQ=255)
                                                                                                                                      acgacgATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTgg           >  2:1266019/1‑150 (MQ=255)
                                                                                                                                        gacgaTGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGa                                              >  1:596602/1‑113 (MQ=255)
                                                                                                                                        gacgaTGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGa                                              <  2:596602/113‑1 (MQ=255)
                                                                                                                                           gaTGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAAACCTACCCGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGcc      >  1:1601484/1‑150 (MQ=255)
                                                                                                                                             tGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGa    <  1:17854/150‑1 (MQ=255)
                                                                                                                                                cAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGAcc  <  2:1055072/149‑1 (MQ=255)
                                                                                                                                                cAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGAcc  >  1:1055072/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCCGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACC  >  NC_002947/1289711‑1290003

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