Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 107,553 (T)9→8 100% intergenic (+79/‑152) lpxC → / → secM UDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913107,5450T.95.8% 309.6 / 6.6 72intergenic (+71/‑160)lpxC/secMUDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base . (37/32);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

TCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTG  >  NC_000913/107446‑107637
                                                                                                   |                                                                                            
tCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGt                                                                                              <  2:1646927/100‑1 (MQ=255)
   aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttt                                                                                           >  1:1247677/1‑100 (MQ=255)
   aaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttt                                                                                           >  1:1905767/1‑100 (MQ=255)
    aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttt                                                                                          >  2:36708/1‑100 (MQ=255)
      cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                        <  2:2222968/100‑1 (MQ=255)
       gcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttttt                                                                                       <  2:1186670/100‑1 (MQ=255)
       gcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttttt                                                                                       <  2:1568822/100‑1 (MQ=255)
        ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                      <  1:813759/100‑2 (MQ=255)
         cTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttta                                                                                     >  2:1962881/1‑99 (MQ=255)
           tCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttgaat                                                                                   <  1:922317/100‑5 (MQ=255)
            cAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAtt                                                                                 <  2:2310118/100‑1 (MQ=255)
             aGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAttt                                                                                >  1:389933/1‑100 (MQ=255)
              gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAtttt                                                                               >  2:1595621/1‑100 (MQ=255)
              gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAtttt                                                                               >  2:1617823/1‑100 (MQ=255)
                tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTtat                                                                             <  1:928761/100‑1 (MQ=255)
                tGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTtat                                                                             <  2:356979/100‑1 (MQ=255)
                 gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTtata                                                                            >  1:721863/1‑100 (MQ=255)
                  tACTGGCATAACGACATTTATACTGTCGTATAAAATGCGACTGGCAAATCTGGCTCTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTtata                                                                            <  2:2033990/99‑1 (MQ=255)
                   aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAag                                                                          >  1:268145/1‑100 (MQ=255)
                   aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAag                                                                          <  1:614184/100‑1 (MQ=255)
                    cTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAaga                                                                         >  2:2144073/1‑100 (MQ=255)
                     tGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAagag                                                                        <  1:899300/100‑1 (MQ=255)
                     tGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAagag                                                                        <  1:2068565/100‑1 (MQ=255)
                      ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                       <  2:2304434/100‑1 (MQ=255)
                      ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                       >  1:883190/1‑100 (MQ=255)
                       gCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCt                                                                      >  1:2029199/1‑100 (MQ=255)
                        cATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCta                                                                     >  2:14418/1‑100 (MQ=255)
                         aTAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCtat                                                                    >  2:667086/1‑100 (MQ=255)
                          tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCtata                                                                   <  2:2120833/100‑1 (MQ=255)
                          tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCtata                                                                   >  1:1745843/1‑100 (MQ=255)
                          tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCtata                                                                   >  1:329465/1‑100 (MQ=255)
                           aaCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaa                                                                  >  1:1888388/1‑100 (MQ=255)
                               aCATTTATACTGTCGTATAAAATTCGACTGGAAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaaaaa                                                               <  2:695430/99‑1 (MQ=255)
                                cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                             <  2:2056806/100‑1 (MQ=255)
                                 aTTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACg                                                            >  1:240990/1‑100 (MQ=255)
                                   ttATACTTTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGt                                                          <  1:867201/100‑1 (MQ=255)
                                   ttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGt                                                          <  2:389930/100‑1 (MQ=255)
                                    tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                         >  1:1365343/1‑100 (MQ=255)
                                    tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                         <  1:2272566/100‑1 (MQ=255)
                                    tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                         >  1:445211/1‑100 (MQ=255)
                                    tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTACGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                         >  2:255824/1‑100 (MQ=255)
                                        cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                     >  2:176732/1‑100 (MQ=255)
                                        cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGa                                                      >  2:1784786/1‑99 (MQ=255)
                                          gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                   <  1:480745/100‑1 (MQ=255)
                                             gTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTg                                                >  1:640774/1‑100 (MQ=255)
                                               ataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtt                                              >  1:844736/1‑100 (MQ=255)
                                                taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                             <  1:355525/100‑1 (MQ=255)
                                                 aaaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGCGCTATAAAAAACGGTGCGAACGCTGtttt                                            >  2:1414417/1‑100 (MQ=255)
                                                 aaaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtttt                                            <  1:1056276/100‑1 (MQ=255)
                                                    aTTCGACTGGCAAATCTGGCACTCACTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                         <  1:327005/100‑1 (MQ=255)
                                                     ttCGGCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACAAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                        <  2:378025/100‑1 (MQ=255)
                                                      tCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                       <  1:1937955/100‑1 (MQ=255)
                                                       tctcTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                      <  1:2217405/97‑1 (MQ=255)
                                                        gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                       <  1:361248/98‑1 (MQ=255)
                                                          cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAgaga                                   >  2:2217386/1‑97 (MQ=255)
                                                          cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                   >  1:1498385/1‑100 (MQ=255)
                                                          cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                   <  2:951961/100‑1 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                 >  2:1383163/1‑100 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                 >  2:127609/1‑100 (MQ=255)
                                                            ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                 >  1:843995/1‑100 (MQ=255)
                                                             gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                >  1:1901253/1‑100 (MQ=255)
                                                                 aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                            >  2:531724/1‑100 (MQ=255)
                                                                    tGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                         <  2:422533/100‑1 (MQ=255)
                                                                    tGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                         >  1:1249759/1‑100 (MQ=255)
                                                                    tGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                         >  2:748495/1‑100 (MQ=255)
                                                                       cACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                      <  2:1017330/100‑1 (MQ=255)
                                                                        aCTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACt                     <  1:947373/100‑1 (MQ=255)
                                                                          tctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTa                   <  1:2162047/100‑1 (MQ=255)
                                                                           ctctcCGGCCAGGTGAACCAGTCGTTTTTTTTTTTTTTTTTAAGAGCTATAAAAGACAGTGCGATCTCTGTTTTCTTAAGCACTTTTCCGTACACCttct                  >  2:1036722/1‑98 (MQ=255)
                                                                           ctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGTATTTTATAAGAGCTCTAACAAACGGTGCGACCGATGTTTTCTTAAGCACTTTTCTGCACAACTTAt                  >  2:628490/1‑100 (MQ=255)
                                                                           ctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                  >  2:1000237/1‑100 (MQ=255)
                                                                           ctctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                  >  1:1068033/1‑100 (MQ=255)
                                                                           ctctcCAGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGATATAAAAAACAGTGCGAACGCTGATTTCTTAAGCACTGTTCCGCACAACTTAt                  >  2:1036650/1‑100 (MQ=255)
                                                                             ctcCGGCCAGGTGAACCAGTCGTTTTTTTTTTTTTTTTTAAGAGCTATAAAAAACGGTGCTAATTCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                  >  2:1589930/1‑98 (MQ=255)
                                                                             ctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                <  2:844732/100‑1 (MQ=255)
                                                                                 ggCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAt            >  2:821960/1‑100 (MQ=255)
                                                                                    cAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg         <  2:388491/100‑1 (MQ=255)
                                                                                      ggTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg       <  2:1607337/100‑1 (MQ=255)
                                                                                         gAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtg    <  1:303917/100‑1 (MQ=255)
                                                                                          aaCCAGTCGTTTTTTTTTTTTTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgt   >  2:2099897/1‑100 (MQ=255)
                                                                                          aaCCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgt   >  2:680853/1‑100 (MQ=255)
                                                                                           aCCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACTACTTATCTTCATTCGTGCtgtg  <  2:1365335/100‑1 (MQ=255)
                                                                                                   |                                                                                            
TCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTG  >  NC_000913/107446‑107637

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