Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,070,462 T→C K165R (AAG→AGG)  rutE ← putative malonic semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,070,4620TC98.3% 212.9 / ‑4.3 59K165R (AAG→AGG) rutEputative malonic semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (32/26);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

CGAAAAGTTTGTTGATCGACAATGTTCATAGTTTGCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCTTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGAACTGTTGCGAAAC  >  NC_000913/1070326‑1070593
                                                                                                                                        |                                                                                                                                   
cgAAAAGTTTGTTGATCGACAATGTTCATAGTTTGCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCa                                                                                                                                  >  1:1018599/1‑140 (MQ=255)
cgAAAAGTTTGTTGATCGACAATGTTCATAGTTTGCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCa                                                                                                                                  >  1:1185279/1‑140 (MQ=255)
          gttgATCGACAATGTTCATAGTTTGCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGcc                                                                                                                        <  2:1018599/140‑1 (MQ=255)
                 gACAATGTTCATAGTTTGCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAGAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTaaaa                                                                                                                 >  2:1099603/1‑140 (MQ=255)
                 gACAATGTTCATAGTTTGCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTaaaa                                                                                                                 <  2:509843/140‑1 (MQ=255)
                 gACAATGTTCATAGTTTGCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTaaaa                                                                                                                 <  1:281631/140‑1 (MQ=255)
                       gTTCATAGTTTGCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGcg                                                                                                           <  1:1099603/140‑1 (MQ=255)
                          cATAGTTTGCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGcgtcg                                                                                                        <  1:188917/140‑1 (MQ=255)
                            tAGTTTGCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCgtcgtc                                                                                                      <  1:375509/140‑1 (MQ=255)
                                  gCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACgtg                                                                                                >  2:828672/1‑140 (MQ=255)
                                  gCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACgtg                                                                                                >  1:523716/1‑140 (MQ=255)
                                                cccGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGc                                                                                  >  1:411424/1‑140 (MQ=255)
                                                cccGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGc                                                                                  >  2:262436/1‑140 (MQ=255)
                                                cccGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGc                                                                                  >  2:34269/1‑140 (MQ=255)
                                                cccGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGc                                                                                  >  2:1063088/1‑140 (MQ=255)
                                                cccGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGc                                                                                  >  1:1208026/1‑140 (MQ=255)
                                                cccGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTAAAAGc                                                                                  >  1:359035/1‑140 (MQ=255)
                                                  cGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGccc                                                                                >  2:304811/1‑140 (MQ=255)
                                                  cGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGccc                                                                                >  1:356247/1‑140 (MQ=255)
                                                  cGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAACCTTCCTGCTCTCGCCATAGCCCATCTTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACCTATTGACGGTCAAAGccc                                                                                >  2:983232/1‑140 (MQ=255)
                                                           cttcAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCGATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCGCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCggg                                                                       <  1:983232/140‑1 (MQ=255)
                                                           cttcAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCggg                                                                       >  1:404005/1‑140 (MQ=255)
                                                           cttcAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCggg                                                                       <  2:404005/140‑1 (MQ=255)
                                                            ttcAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACAt                                                                           >  1:438570/1‑135 (MQ=255)
                                                            ttcAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACAt                                                                           <  2:438570/135‑1 (MQ=255)
                                                            ttcAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGc                                                                      <  2:1226273/140‑1 (MQ=255)
                                                             tcAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGcc                                                                     <  2:1185279/140‑1 (MQ=255)
                                                              cAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCg                                                                    <  2:1024904/140‑1 (MQ=255)
                                                                     cAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCa                                                             >  2:896501/1‑140 (MQ=255)
                                                                      aGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAg                                                            >  2:631116/1‑140 (MQ=255)
                                                                            ggCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTAAAAGCCCGACATCGGGCCGGTATCCAGTCCCAg                                                      <  2:359035/140‑1 (MQ=255)
                                                                             gCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTAt                                                                <  2:131929/129‑1 (MQ=255)
                                                                             gCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTAt                                                                >  1:131929/1‑129 (MQ=255)
                                                                             gCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAgc                                                     >  2:927642/1‑140 (MQ=255)
                                                                             gCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAgc                                                     >  2:657384/1‑140 (MQ=255)
                                                                                  cgcgcATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCg                                                >  1:1029734/1‑140 (MQ=255)
                                                                                        tAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGc                                          <  2:523716/140‑1 (MQ=255)
                                                                                        tAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGc                                          <  2:584524/140‑1 (MQ=255)
                                                                                        tAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGc                                          <  2:1029734/140‑1 (MQ=255)
                                                                                        tAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGc                                          <  1:631116/140‑1 (MQ=255)
                                                                                        tAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGc                                          <  2:128912/140‑1 (MQ=255)
                                                                                          aaGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGcga                                        <  1:927642/140‑1 (MQ=255)
                                                                                            gCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGcgacg                                      <  2:498112/140‑1 (MQ=255)
                                                                                                 ctgctATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAg                                 >  1:1263077/1‑140 (MQ=255)
                                                                                                 ctgctATCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAg                                 >  1:734512/1‑140 (MQ=255)
                                                                                                       tCGCCATAGCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGgc                           >  2:803347/1‑140 (MQ=255)
                                                                                                               gCCGATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCa                   <  1:744064/140‑1 (MQ=255)
                                                                                                                  gATATTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATgg                <  1:262436/140‑1 (MQ=255)
                                                                                                                   atatTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGa               >  1:1170127/1‑140 (MQ=255)
                                                                                                                   atatTAATCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGa               >  1:1067948/1‑140 (MQ=255)
                                                                                                                         aTCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGAACTGtt         >  2:348671/1‑140 (MQ=255)
                                                                                                                         aTCAGCAGATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGAACTGtt         >  2:698966/1‑140 (MQ=255)
                                                                                                                            agcagATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGAACTGTTGCg      >  1:1216384/1‑140 (MQ=255)
                                                                                                                            agcagATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGAACTGTTGCg      <  1:348671/140‑1 (MQ=255)
                                                                                                                              cagATTGCTCTTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGAACTGTTGCGaa    >  2:96917/1‑140 (MQ=255)
                                                                                                                              cagATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGaa              <  1:245218/130‑1 (MQ=255)
                                                                                                                              cagATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGaa              >  2:245218/1‑130 (MQ=255)
                                                                                                                              cagATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGAACTg           >  2:636597/1‑133 (MQ=255)
                                                                                                                              cagATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGAACTg           <  1:636597/133‑1 (MQ=255)
                                                                                                                                gATTGCTCCTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGAACTGTTGCGAAAc  >  1:10769/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                   
CGAAAAGTTTGTTGATCGACAATGTTCATAGTTTGCTCCTTACAACAGCCCGCAGGCTTCTTCAAAGGACAGACGTGGCAGGCGCGCATAAAGCTTGCTGCTATCGCCATAGCCGATATTAATCAGCAGATTGCTCTTCAGCGTGCTGCCCGTAAAAAAGGCGTCGTCCACGTGTTGACGGTCAAAGCCCGACATCGGGCCGGTATCCAGTCCCAGCGCCCGGCAGGCGACGATCAGATAGGCCGCCTGCATGGAACTGTTGCGAAAC  >  NC_000913/1070326‑1070593

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