Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,792,242:1 +A 17.5% intergenic (‑222/+25) ydiV ← / ← nlpC anti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,792,2421.A17.5% 241.1 / 58.0 80intergenic (‑222/+25)ydiV/nlpCanti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (29/37);  new base A (8/6);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

TTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCT  >  NC_000913/1792147‑1792341
                                                                                                |                                                                                                   
tttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatcatc                                                                                                 >  2:1727789/1‑100 (MQ=255)
tttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatcatc                                                                                                 >  2:548826/1‑100 (MQ=255)
 ttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTGTGTATGCAAAATACTACGGGAGTGAGGTGTATATGGACGGCGATGTATCATCg                                                                                                <  2:1365356/100‑1 (MQ=255)
  tttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGc                                                                                               >  1:1695813/1‑100 (MQ=255)
  tttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGc                                                                                               >  1:264292/1‑100 (MQ=255)
   ttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGcc                                                                                              <  1:1074427/100‑1 (MQ=255)
   ttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGcc                                                                                              >  1:2255049/1‑100 (MQ=255)
    tAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCg                                                                                             <  2:1932698/100‑1 (MQ=255)
    tAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCg                                                                                             >  2:502121/1‑100 (MQ=255)
     aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGt                                                                                            <  2:330929/100‑1 (MQ=39)
     aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGt                                                                                            <  1:1757939/100‑1 (MQ=39)
      gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                           <  1:31147/100‑1 (MQ=255)
      gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                           <  1:1585326/100‑1 (MQ=39)
       aTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTcc                                                                                          <  2:1863235/100‑1 (MQ=39)
       aTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTcc                                                                                          <  1:1101720/100‑1 (MQ=39)
        ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCACCGTCCa                                                                                         <  1:1715673/100‑1 (MQ=38)
         tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAt                                                                                        <  1:295541/100‑1 (MQ=255)
         tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAt                                                                                        <  2:1408100/100‑1 (MQ=255)
         tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAt                                                                                        >  2:1724742/1‑100 (MQ=255)
             cTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAttttt                                                                                    >  2:2165128/1‑100 (MQ=255)
             cTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAttttt                                                                                    <  2:1571455/100‑1 (MQ=255)
              tGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTc                                                                                   <  1:1913327/100‑1 (MQ=255)
                 ttATAAACGAATAGTCCATTTACGGAATTTTTGTCTGCAAAAGACTACTGTATTCAGTAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGAc                                                                                <  2:2294958/100‑1 (MQ=255)
                   ataAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACAt                                                                              >  1:1732980/1‑100 (MQ=255)
                    taAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATa                                                                             >  2:2379422/1‑100 (MQ=255)
                     aaaCGAATAATCCATTTACGGACTTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGGCGGCGATGTATCATCCCCGTCCATTTTTCGACATAc                                                                            >  2:617149/1‑100 (MQ=255)
                     aaaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATAc                                                                            <  2:2134762/100‑1 (MQ=255)
                      aaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACt                                                                           >  1:1821036/1‑100 (MQ=255)
                      aaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGAAGATGTATTATCGACGTCCATTTTTCGACATACt                                                                           >  2:2181058/1‑100 (MQ=255)
                       aCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTa                                                                          <  2:2319108/100‑1 (MQ=255)
                          aataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGAt                                                                       <  1:25534/100‑1 (MQ=255)
                                 cATTTACGGAATTTTTGTTTGTAAAATACTACTGTATTCAGGAGTAAATGTACGTCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGtctc                                                                >  2:2280470/1‑100 (MQ=255)
                                    ttACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGcc                                                             <  2:2001188/100‑1 (MQ=255)
                                     tACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCt                                                            >  2:112233/1‑100 (MQ=255)
                                        ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGcc                                                         <  2:1353109/100‑1 (MQ=255)
                                               ttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAt                                                   <  2:1257640/99‑1 (MQ=255)
                                                 gTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAtttt                                                <  1:1817040/100‑1 (MQ=255)
                                                  tCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAttttt                                               <  2:1236438/100‑1 (MQ=255)
                                                   cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                              <  2:1360811/100‑1 (MQ=255)
                                                     gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                            <  1:2269723/100‑1 (MQ=255)
                                                     gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                            <  2:658675/100‑1 (MQ=255)
                                                          aTACTCCTGTATTCAGGAGTAAATGGCAGGCTATGTATCATCGTCGTCCATTTTTCGAGCTACGACATTCGTCCCGCTTGCCAGAAAATTTTCTGCCa                                         <  2:2372344/98‑1 (MQ=255)
                                                             ctactGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGa                                    >  1:2074422/1‑100 (MQ=255)
                                                               actGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                  >  2:399689/1‑100 (MQ=255)
                                                               actGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                  >  2:1221818/1‑100 (MQ=255)
                                                                ctGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAt                                 >  1:382813/1‑100 (MQ=255)
                                                                ctGTATCCAGGAGTAAATGGATGGCGCTGTATCAGCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCGGAAATTTTTCTGCCAATAg                                     <  2:1463531/96‑1 (MQ=255)
                                                                 tGTATTCAGGAGTAAATGGACAGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAtt                                <  2:1719990/100‑1 (MQ=255)
                                                                  gTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTa                               >  1:900602/1‑100 (MQ=255)
                                                                   tATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTAt                              <  1:502123/100‑1 (MQ=255)
                                                                   tATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTAt                              <  2:13334/100‑1 (MQ=255)
                                                                    aTTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATc                             >  1:447685/1‑100 (MQ=255)
                                                                    aTTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATc                             <  1:2036171/100‑1 (MQ=255)
                                                                     ttCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCa                            >  1:339211/1‑100 (MQ=255)
                                                                     ttCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCa                            >  2:354302/1‑100 (MQ=255)
                                                                       cAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                          >  2:1710112/1‑100 (MQ=255)
                                                                       cAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                          >  1:810957/1‑100 (MQ=255)
                                                                        aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGt                         <  1:2058393/100‑1 (MQ=255)
                                                                         ggAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                        <  1:730974/100‑1 (MQ=255)
                                                                         ggAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                        <  1:1033928/100‑1 (MQ=255)
                                                                           aGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAg                      >  2:1397472/1‑100 (MQ=255)
                                                                            gTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAgg                     <  2:1085900/100‑1 (MQ=255)
                                                                             tAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGa                    <  2:311976/100‑1 (MQ=255)
                                                                             tAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGa                    <  2:197534/100‑1 (MQ=255)
                                                                               aaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc                  <  1:1068702/100‑1 (MQ=255)
                                                                               aaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc                  >  2:1218487/1‑100 (MQ=255)
                                                                                 tGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGc                <  2:852339/100‑1 (MQ=255)
                                                                                  ggACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa               <  1:1201931/100‑1 (MQ=255)
                                                                                  ggACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa               <  1:194395/100‑1 (MQ=255)
                                                                                  ggACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATCTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa               >  1:177080/1‑100 (MQ=255)
                                                                                  ggACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa               >  1:382075/1‑100 (MQ=39)
                                                                                   gACGGCGATGTATCATCGCCGTCCATTTTTTTATATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  2:5389/1‑100 (MQ=255)
                                                                                   gACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:2372957/1‑100 (MQ=255)
                                                                                   gACGGCGATGATACATCGCCGTCCATTTTTCTACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:715899/1‑100 (MQ=38)
                                                                                   gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:2236596/1‑100 (MQ=39)
                                                                                   gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:1222050/1‑100 (MQ=39)
                                                                                   gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:1848807/1‑100 (MQ=39)
                                                                                   gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:2000857/1‑100 (MQ=39)
                                                                                   gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCAAGAAATTTTTTTGAGAATAGACATTATCAAGTGAGGAACGAAt              >  2:1468068/1‑100 (MQ=37)
                                                                                    aCGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATc             <  2:2255029/100‑1 (MQ=255)
                                                                                    aCGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATc             >  1:33839/1‑100 (MQ=39)
                                                                                             tATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTg    >  2:2268486/1‑100 (MQ=255)
                                                                                              atcatcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGc   <  1:112234/100‑1 (MQ=255)
                                                                                              atcatcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACACTATCAAGTGAGGAACGCATCACTCCCTTGc   >  2:2294534/1‑100 (MQ=255)
                                                                                               tcatcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCt  <  2:1732966/100‑1 (MQ=255)
                                                                                                |                                                                                                   
TTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCT  >  NC_000913/1792147‑1792341

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: