Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,972,637 2 bp→GC coding (430‑431/546 nt) wecD → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,972,6370CG97.7% 324.0 / NA 130R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base G (64/61);  total (66/64)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,972,6380GC99.2% 338.1 / NA 130R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (65/64);  total (66/64)
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

Only 100 of 131 total aligned reads displayed.

GCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATAC  >  NC_000913/3972567‑3972706
                                                                      ||                                                                    
gcgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTcgc                                                                  >  1:5117582/1‑76 (MQ=255)
  gAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCgg                                                                <  1:3130640/76‑1 (MQ=255)
    aTTGGCGTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTa                                                              <  2:256139/76‑1 (MQ=255)
     ttGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTaa                                                             <  2:2977865/76‑1 (MQ=255)
     ttGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCCTATGCTCGCGGTaa                                                             <  2:4736266/76‑1 (MQ=255)
     tgggCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTaa                                                             <  1:696792/74‑1 (MQ=255)
      tGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTaaa                                                            <  1:808381/76‑1 (MQ=255)
      tGGCCTCCTGGCTGGACGAGGTGCAGGGGCTGAGCTGATGCAAACGGCGCTGAACTGCGTGTAGTGTCGAGGTata                                                            >  1:4454822/1‑74 (MQ=255)
       ggCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTaaaa                                                           >  1:3040394/1‑76 (MQ=255)
        gCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaac                                                          <  2:3532822/76‑1 (MQ=255)
         ccTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaaca                                                         >  1:4518316/1‑76 (MQ=255)
          ctgctgGCTGGACGCGGTGCAGGTGCTGCGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaacaa                                                        <  2:4829553/76‑1 (MQ=255)
          ctgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAcacaa                                                        >  1:5731044/1‑76 (MQ=255)
           tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaacaac                                                       >  1:4567501/1‑76 (MQ=255)
            gctggctgGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACt                                                      <  2:2537708/76‑1 (MQ=255)
             ctggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGGAAAACAACtt                                                     >  2:4912518/1‑76 (MQ=255)
              tggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACttt                                                    <  2:4742466/76‑1 (MQ=255)
              tggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACttt                                                    <  1:5362438/76‑1 (MQ=255)
                gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGc                                                  <  1:4597212/76‑1 (MQ=255)
                gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGc                                                  <  1:2389794/76‑1 (MQ=255)
                 ctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCg                                                 <  1:883784/76‑1 (MQ=255)
                    gACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAGCTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGt                                              <  2:4873534/76‑1 (MQ=255)
                    gACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGt                                              <  1:5469788/76‑1 (MQ=255)
                    gACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGt                                              <  1:1143251/76‑1 (MQ=255)
                    gACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGt                                              >  1:3216990/1‑76 (MQ=255)
                     gcgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTg                                             <  2:3767620/75‑1 (MQ=255)
                     gcgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTg                                             <  2:1211543/75‑1 (MQ=255)
                     aCGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTg                                             >  2:5525582/1‑76 (MQ=255)
                     aCGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTg                                             >  2:3486567/1‑76 (MQ=255)
                     aCGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTg                                             <  2:1381414/76‑1 (MQ=255)
                      cgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTgg                                            <  2:78926/76‑1 (MQ=255)
                      cgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTgg                                            <  2:2805616/76‑1 (MQ=255)
                       gcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGc                                           >  2:4259920/1‑76 (MQ=255)
                       gcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGc                                           >  1:1145746/1‑76 (MQ=255)
                         ggTGGAGGGGCGGAGCTGGTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGa                                         <  2:233962/76‑1 (MQ=255)
                         ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAACACAACTTTGCGGGGGGCgc                                         >  2:5819699/1‑75 (MQ=255)
                         ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGa                                         >  2:2837023/1‑76 (MQ=255)
                         ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGa                                         >  1:3201533/1‑76 (MQ=255)
                         ggTGCAGGGGCGGAGCTGAGGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGa                                         <  2:5443874/76‑1 (MQ=255)
                          gTGCAGGTGCTGCGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGAc                                        <  1:5127167/76‑1 (MQ=255)
                           tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGAcc                                       >  2:1858779/1‑76 (MQ=255)
                           tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGCGGAGAcc                                       >  1:2300007/1‑76 (MQ=255)
                           tGCAGGTGCGGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGAca                                       <  2:4889208/76‑2 (MQ=255)
                            gCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGGGGCGAccc                                      >  2:142773/1‑76 (MQ=255)
                             cAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCa                                     >  2:2385138/1‑76 (MQ=255)
                              aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAGCTTTGCGGGTGGCGTCCCaa                                    <  2:3420676/76‑1 (MQ=255)
                              aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCaa                                    >  2:2277872/1‑76 (MQ=255)
                              aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCaa                                    >  1:4730019/1‑76 (MQ=255)
                              aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCaa                                    >  2:4257131/1‑76 (MQ=255)
                               ggTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCaca                                   >  1:2618272/1‑74 (MQ=255)
                               ggTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCaaa                                   >  1:2792163/1‑76 (MQ=255)
                               ggTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCaaa                                   >  2:4435956/1‑76 (MQ=255)
                                gTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAAt                                  <  1:5648692/76‑1 (MQ=255)
                                  gctgagctgaTGCAAACGGCGCTAAACTGGGCGCATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAACgg                                >  1:4022595/1‑76 (MQ=255)
                                    tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGc                              <  1:1017148/76‑1 (MQ=255)
                                     gagctgaTGCAAACGTCGCTAAACGAGTCTTATGCTCGCGGTAAACCAACGTTGCGGGTTGCGACCCAACTGTGCa                             <  1:2664159/76‑1 (MQ=255)
                                     gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCa                             <  2:927955/76‑1 (MQ=255)
                                      ggctaaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCaa                            <  2:5731042/75‑1 (MQ=255)
                                       gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAac                           <  1:5002378/76‑1 (MQ=255)
                                        ctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAaca                          >  2:3874461/1‑76 (MQ=255)
                                        ctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAaca                          >  1:5736975/1‑76 (MQ=255)
                                           gtGCAGACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccg                       <  2:215507/75‑1 (MQ=255)
                                           aTGCAAACGGCGCTAACCTGGGCTTATTCTCGCGTTAAAACAACTTTGCGGGGGTAGACCACAATGTGCACGAccg                       >  2:2727275/1‑76 (MQ=255)
                                           aTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccg                       >  1:608573/1‑76 (MQ=255)
                                           aTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccg                       >  2:4927557/1‑76 (MQ=255)
                                           aTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAACGGGCAACAccg                       >  2:2699638/1‑76 (MQ=255)
                                             gCAACCGGCGCGAGACTGGGCGTATGCTCGCGGTAAAACCGCTTTGCGGGTGGCGGCCCAAATGGGCAACAgcgcc                     <  1:267100/76‑1 (MQ=255)
                                             gCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgcc                     <  1:5525584/76‑1 (MQ=255)
                                             gCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgcc                     <  2:2956394/76‑1 (MQ=255)
                                              aaaaCGGCGCGAAGCGGGGCGTAGGCTCGCGGTAAAACACCTTGGGGGGTGGCGCCCAAAAGGGGCAACAccgccg                    <  2:5117581/75‑1 (MQ=255)
                                               aaaCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                   >  1:3780863/1‑76 (MQ=255)
                                                aggggCGCTAAGGGGGGCGTATGCTGGGGGTAAAACAACTTGGGGGGGGGCGACCCAAAGGGGCAACACCGCcgcg                  <  2:4287551/73‑1 (MQ=255)
                                                accGGCGCAAGACTGGGCGGAGGCTGGCGGTAAAACAATTTGGCGGGTGGCGCCCAAAAGGGGCAACACCGCcgcg                  <  2:3123634/74‑1 (MQ=255)
                                                aaCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCcgcg                  >  2:3985863/1‑76 (MQ=255)
                                                 aCGGCGCTAAACTGGGCGTATGCTAGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCgcgc                 >  2:3326321/1‑76 (MQ=255)
                                                  cGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCt                >  2:823452/1‑76 (MQ=255)
                                                  cGGCGCTAAACGGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCt                <  1:1691591/76‑1 (MQ=255)
                                                   ggCGCTAACGTGGGCGTATGCTCGCGGTAAAACAACGTTGCGGTGGGCGACCCAAATGGGCAACACCGCCGCGCtt               <  1:3626562/76‑1 (MQ=255)
                                                   ggCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt               <  1:4883464/76‑1 (MQ=255)
                                                   ggCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt               >  1:1734646/1‑76 (MQ=255)
                                                   ggCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt               >  2:5304131/1‑76 (MQ=255)
                                                   gccgcAAAACGGGGCGTATCCCCGCGGAAAAACACCTTGGCGGGGGGCGACCCAAAGGGGCAACACCGCCGCGCtt               <  1:1989849/74‑1 (MQ=255)
                                                    gcgcTAAACTGGGCGTATGCTCGCGGTACAACAACTTTGCGGGTGGATACCCAAATTTGCAACACCGCCGAGCCTa              >  1:1485666/1‑76 (MQ=255)
                                                    gcgcTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa              >  1:1441239/1‑76 (MQ=255)
                                                    gcgcGAAACGGGGGGTAGCCGCGCGGTAAAACAACTTTGCGGGTGGCGACCAAAAGGGGCAACACCGCCGCGCTTa              <  2:1145746/76‑1 (MQ=255)
                                                     ggcTAAACGGGGCGTAGGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaa             <  2:3190285/75‑1 (MQ=255)
                                                     ggcTAAACGGGGCGTAGGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaa             <  2:3159511/75‑1 (MQ=255)
                                                     cgcTAACCGGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGCCCCAAATGGGCAACACCGCCGCGCTTaa             <  1:3983023/76‑1 (MQ=255)
                                                      gcTAAACTGGGCGTATGCTTGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa            >  1:4072583/1‑76 (MQ=255)
                                                       cTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAc           >  2:1209590/1‑76 (MQ=255)
                                                        tAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAcn          <  1:3591325/76‑2 (MQ=255)
                                                        tAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACg          >  1:3062561/1‑76 (MQ=255)
                                                        tAAACGGGGCGTAGCCTCGGGGAAAACCAACTTGGGGGGTGGCGACCAAAGTGGGCAACCCCGCCGCGCTTAAACg          <  2:1610880/76‑1 (MQ=255)
                                                         aaaCTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGa         >  1:1855574/1‑76 (MQ=255)
                                                           aCTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGata       >  1:1750037/1‑76 (MQ=255)
                                                            cTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac      >  1:354706/1‑76 (MQ=255)
                                                              gggCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacat    >  1:4078167/1‑76 (MQ=255)
                                                              gggCGTAGGGTGGCGGTAAAACAACTTTGCGGGGGGCGACCAAAAGGGGAAACACCGCCGCGCTTAAACGatacat    <  2:5404472/76‑1 (MQ=255)
                                                               ggCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacata   >  2:4839664/1‑76 (MQ=255)
                                                                gCGTATGCTCGCGGAAAAACAACTTGGCGGGTGGCGCCCAAAATGGGCAACACCGCCGCGCTTAAACGatacatac  <  1:2450955/76‑1 (MQ=255)
                                                                      ||                                                                    
GCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATAC  >  NC_000913/3972567‑3972706

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