Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,897,050 +G coding (277/666 nt) yieH → phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,897,0461.G96.4% 196.4 / ‑0.2 56coding (273/666 nt)yieHphosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (28/26);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCA  >  NC_000913/3896911‑3897173
                                                                                                                                        |                                                                                                                                 
cAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAggg                                                                                                                               <  2:922674/139‑1 (MQ=255)
   aGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAgggggc                                                                                                                           <  2:902857/139‑3 (MQ=255)
        tAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAg                                                                                                                       >  2:973587/1‑139 (MQ=255)
             cTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTc                                                                                                                  <  1:959572/139‑1 (MQ=255)
               gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCt                                                                                                                >  2:739405/1‑139 (MQ=255)
               gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCt                                                                                                                >  1:288156/1‑139 (MQ=255)
                        tATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGAt                                                                                                       >  1:101928/1‑139 (MQ=255)
                          tCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                     >  1:583510/1‑139 (MQ=255)
                          tCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                     >  2:620136/1‑139 (MQ=255)
                            gATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACt                                                                                                   <  1:227840/139‑1 (MQ=255)
                                 tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                    <  2:162179/103‑1 (MQ=255)
                                 tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                    >  1:162179/1‑103 (MQ=255)
                                 tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGcc                                                                                              >  1:938384/1‑139 (MQ=255)
                                  gTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                             <  1:739405/139‑1 (MQ=255)
                                     tCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtg                                                                                          >  1:609405/1‑139 (MQ=255)
                                     tCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtg                                                                                          <  2:583510/139‑1 (MQ=255)
                                        cTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgt                                                                                       >  1:172183/1‑139 (MQ=255)
                                         ttGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                      <  2:288156/139‑1 (MQ=255)
                                         ttGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                      <  1:620136/139‑1 (MQ=255)
                                           gAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGa                                                                                                                        >  1:474544/1‑103 (MQ=255)
                                           gAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGa                                                                                                                        <  2:474544/103‑1 (MQ=255)
                                                      tACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                    <  2:920560/82‑1 (MQ=255)
                                                      tACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                    >  1:920560/1‑82 (MQ=255)
                                                      tACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACgg                                                                         >  1:940194/1‑139 (MQ=255)
                                                      tACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACgg                                                                         >  2:875512/1‑139 (MQ=255)
                                                       aCGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGc                                                                        <  2:885132/139‑1 (MQ=255)
                                                           tAGCGTAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                >  1:696439/1‑111 (MQ=255)
                                                           tAGCGTAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                <  2:696439/111‑1 (MQ=255)
                                                           tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTg                                                                                                  >  1:328297/1‑109 (MQ=255)
                                                           tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTg                                                                                                  <  2:328297/109‑1 (MQ=255)
                                                           tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCaa                                                                    >  1:660847/1‑139 (MQ=255)
                                                            aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCaaa                                                                   >  2:594058/1‑139 (MQ=255)
                                                             gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgc                                                                                                 >  2:614033/1‑108 (MQ=255)
                                                             gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgc                                                                                                 <  1:614033/108‑1 (MQ=255)
                                                             gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGc                                                                        <  1:973891/133‑1 (MQ=255)
                                                             gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGc                                                                        >  2:973891/1‑133 (MQ=255)
                                                             gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAAt                                                                  >  2:687735/1‑139 (MQ=255)
                                                                         gCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcag                                                      <  2:495727/139‑1 (MQ=255)
                                                                             aaCACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGcc                                                                                              <  2:786938/95‑1 (MQ=255)
                                                                             aaCACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGcc                                                                                              >  1:786938/1‑95 (MQ=255)
                                                                              acacGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTc                                                 <  1:594058/139‑1 (MQ=255)
                                                                               cacGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCt                                                <  2:132670/139‑1 (MQ=255)
                                                                                   tttACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATgg                                            <  1:248932/139‑1 (MQ=255)
                                                                                    ttACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGAGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATggg                                           <  1:560862/139‑1 (MQ=255)
                                                                                     tACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                          <  2:322414/139‑1 (MQ=255)
                                                                                         gTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGc                                      >  1:914270/1‑139 (MQ=255)
                                                                                          tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAg                                                                                                                                  >  1:145934/1‑47 (MQ=255)
                                                                                          tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAg                                                                                                                                  <  2:145934/47‑1 (MQ=255)
                                                                                            cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGa                                   >  2:97487/1‑139 (MQ=255)
                                                                                            cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGa                                   <  1:97487/139‑1 (MQ=255)
                                                                                            cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGa                                   <  2:266377/139‑1 (MQ=255)
                                                                                            cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGa                                   >  2:241679/1‑139 (MQ=255)
                                                                                            cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGa                                   <  2:830534/139‑1 (MQ=255)
                                                                                            cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGa                                   >  2:863422/1‑139 (MQ=255)
                                                                                                    gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                    <  2:387779/36‑1 (MQ=255)
                                                                                                    gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                    >  1:387779/1‑36 (MQ=255)
                                                                                                    gCTCGGCTGTTCGAGTCAGAACTGGAAGCCATCGAAGAGGGCTGGAGCCCTCCTGTCAGCGATCACCGCGCCAACGTGTGTGGTATCTATCGGCCCAAATAACTAAATGCAGCAATCGATGGGCAAGCTGAATATGTTg                           >  2:316393/1‑139 (MQ=255)
                                                                                                            gTTCGATTGAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGTGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACtt                   <  2:645697/139‑1 (MQ=255)
                                                                                                            gTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcag                                                      <  2:17655/104‑1 (MQ=255)
                                                                                                            gTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcag                                                      >  1:17655/1‑104 (MQ=255)
                                                                                                              tCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTcc                 <  2:940194/139‑1 (MQ=255)
                                                                                                               cGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTccc                <  1:241679/139‑1 (MQ=255)
                                                                                                                  ttCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGa             >  1:912316/1‑139 (MQ=255)
                                                                                                                            gAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTc   >  2:875102/1‑139 (MQ=255)
                                                                                                                            gAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTc   >  2:595166/1‑139 (MQ=255)
                                                                                                                             aaGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCa  <  1:875512/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                 
CAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCA  >  NC_000913/3896911‑3897173

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