Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,698,695 A→G 6.7% V162V (GTT→GTC malI ← transcriptional repressor of Mal regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,698,6950AG 6.7% 180.4 / 3.2 60V162V (GTT→GTCmalItranscriptional repressor of Mal regulon
Reads supporting (aligned to +/- strand):  ref base A (23/33);  new base G (2/2);  total (25/35)
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.63e-01

ATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGATTCAGTAACAG  >  NC_000913/1698550‑1698838
                                                                                                                                                 |                                                                                                                                               
aTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAAcc                                                                                                                                              <  4:15904/149‑1 (MQ=255)
      gtagAGTTGCACAATCGCCCCCAACCCGTACTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATc                                                                                                                                         <  8:94115/148‑1 (MQ=255)
      gtagAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATc                                                                                                                                         <  3:231773/148‑1 (MQ=255)
        agagTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATc                                                                                                                                         <  5:89905/146‑1 (MQ=255)
        agagTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATc                                                                                                                                         >  6:89905/1‑146 (MQ=255)
         gagTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAc                                                                                                                                      <  7:117541/148‑1 (MQ=255)
         gagTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAc                                                                                                                                      <  1:149203/148‑1 (MQ=255)
           gTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTa                                                                                                                                           >  6:284518/1‑141 (MQ=255)
           gTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTa                                                                                                                                           <  5:284518/141‑1 (MQ=255)
           gTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAcatc                                                                                                                                   <  7:204605/149‑1 (MQ=255)
           gTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAcat                                                                                                                                    >  3:291422/1‑148 (MQ=255)
                       tcattaaCCCGTGCAGAACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCCTTGCGAATGAGATGCTCCGTCCACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAGATCTTCAAGATTAcg                                                                                                                       >  6:109499/6‑148 (MQ=255)
                         cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTgg                                                                                                                     <  3:174041/149‑1 (MQ=255)
                          ccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTgg                                                                                                                     >  4:164811/1‑148 (MQ=255)
                          ccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTgg                                                                                                                     >  3:199128/1‑148 (MQ=255)
                            aaCCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACACCTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCAc                                                                                                                  >  1:329210/1‑149 (MQ=255)
                                      gCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGa                                                                                                         >  4:260365/1‑148 (MQ=255)
                                                 tGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGa                                                                                             >  1:180632/1‑149 (MQ=255)
                                                   aGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATa                                                                                           <  6:10880/149‑1 (MQ=255)
                                                   aGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTCGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGAtc                                                                                           <  5:249769/149‑2 (MQ=255)
                                                   aGGAACTTTGCCCTCCCAGCCAGGCGATCAGCTGATGCACATTGCGAATGAGATTCTCCGTCAACACCTGTGCCGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGCGa                                                                                             >  2:78220/1‑147 (MQ=255)
                                                       aCTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCtt                                                                                       <  3:223842/149‑1 (MQ=255)
                                                        cTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCttt                                                                                      <  5:26224/149‑1 (MQ=255)
                                                            gCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTttc                                                                                   <  1:46279/148‑1 (MQ=255)
                                                              ccTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTttcttc                                                                                <  1:283928/149‑1 (MQ=255)
                                                               cTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTtcttct                                                                               <  4:304247/149‑1 (MQ=255)
                                                                tCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTg                                                                              >  5:116704/1‑149 (MQ=255)
                                                                  ccagccagGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGc                                                                             <  5:304331/148‑1 (MQ=255)
                                                                      ccagGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATcc                                                                        >  2:303847/1‑149 (MQ=255)
                                                                      ccagGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATcc                                                                        <  5:249316/149‑1 (MQ=255)
                                                                             aTCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGTACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAg                                                                 <  5:109499/149‑1 (MQ=255)
                                                                                     aTGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATc                                                           <  6:116704/147‑1 (MQ=255)
                                                                                         ccATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGGACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACt                                                     >  7:342722/1‑149 (MQ=255)
                                                                                         ccATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGGACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACt                                                     >  7:347621/1‑149 (MQ=255)
                                                                                         ccATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACt                                                     >  6:168096/1‑149 (MQ=255)
                                                                                           aTTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTc                                                   <  2:329210/149‑1 (MQ=255)
                                                                                              gCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATCACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCt                                                 <  1:200889/148‑1 (MQ=255)
                                                                                                 aaTGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAg                                                                 >  8:322753/1‑129 (MQ=255)
                                                                                                 aaTGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAg                                                                 <  7:322753/129‑1 (MQ=255)
                                                                                                      gaTGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCg                                        >  6:357031/1‑149 (MQ=255)
                                                                                                           tCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGAt                                                                                                                           <  3:121442/61‑1 (MQ=255)
                                                                                                           tCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGAt                                                                                                                           >  4:121442/1‑61 (MQ=255)
                                                                                                             cGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAAtt                                  <  3:164811/148‑1 (MQ=255)
                                                                                                              gTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGGACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTAc                                <  8:347621/149‑1 (MQ=255)
                                                                                                              gTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGGACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTAc                                <  8:342722/149‑1 (MQ=255)
                                                                                                              gTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTAc                                <  3:260365/149‑1 (MQ=255)
                                                                                                                 aacaacTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACca                              >  6:203180/1‑148 (MQ=255)
                                                                                                                 aacaacTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACca                              >  6:208302/1‑148 (MQ=255)
                                                                                                                   caacTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCacac                           >  4:185433/1‑149 (MQ=255)
                                                                                                                    aacTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCaa                                    <  2:358494/139‑1 (MQ=255)
                                                                                                                    aacTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCaa                                    >  1:358494/1‑139 (MQ=255)
                                                                                                                         tgCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATc                       <  3:185433/147‑1 (MQ=255)
                                                                                                                           cAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGaca                   <  4:291422/149‑1 (MQ=255)
                                                                                                                             gCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAAtt                                  >  4:209742/1‑132 (MQ=255)
                                                                                                                             gCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAAtt                                  <  3:209742/132‑1 (MQ=255)
                                                                                                                             gCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACAcc                 <  5:168096/149‑1 (MQ=255)
                                                                                                                                   aTGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGAtt           <  4:199128/149‑1 (MQ=255)
                                                                                                                                         tgcGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAATTCACCGGGATAGCTTTTTCTTCTGCCATCAGTCGCACTTCATCGCTACTTCCTGCAGAGCCGGCAATTACCAAACCATCTACACCCTGCTTCAGTa      >  3:42102/3‑148 (MQ=255)
                                                                                                                                           cGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGATTCAGTAACa   <  1:303847/149‑1 (MQ=255)
                                                                                                                                            gggCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGATTCAGTAACAg  >  8:157096/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACTGGCACGGGAAGCGAAAATCACCGGGATAGCTTTTTCTTCTGCCATCCGTCGCAGGTCATCGCTACTTCCTGCAGCTCCGGCAATTACCACACCATCGACACCCTGATTCAGTAACAG  >  NC_000913/1698550‑1698838

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: