Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,656,349 C→A 10.5% A53S (GCA→TCA)  USA300HOU_2509 ← possible lactoylglutatione lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,656,3490CA10.5% 181.5 / 2.4 63A53S (GCA→TCA) USA300HOU_2509possible lactoylglutatione lyase
Reads supporting (aligned to +/- strand):  ref base C (25/30);  new base A (2/6);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCATAATATGTTAATGCGTCATCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCTGTCATTGCAGTAACATGGTGAATACCTAGTA  >  CP000730/2656211‑2656486
                                                                                                                                          |                                                                                                                                         
gTTCATAATATGTTAATGCGTCATCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCaa                                                                                                                                         <  2:356211/141‑1 (MQ=255)
gTTCATAATATGTTAATGCGTCATCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCaa                                                                                                                                         <  1:267582/141‑1 (MQ=255)
 ttCATAATATGTTAATGCGTCATCGTTAGGCACTCTAAAAGCCGGTCTTGTAATGGAATTTGTTCATGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCaaa                                                                                                                                        <  1:228161/141‑1 (MQ=255)
   cATAATATGTTAATGCGTCATCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCaaaaa                                                                                                                                      >  1:536503/1‑141 (MQ=255)
       atatGTTAATGCGTCATCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGt                                                                                                                                  <  1:315909/141‑1 (MQ=255)
            ttAATGCGTCATCGTTAGGCAAAATAAAAGACGGTCTTGTAATGGAATTTGTACCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAATAACGTAATGTTTGTCCCTGCCGAACTTAAAAAAAAAAAAAAAAAAGTATGAt                                                                                                                             <  2:493595/141‑1 (MQ=255)
               aTGCGTCATCGTTAGGCACTCTAAAAGACGAACTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCAGCTGCAAAAAAAGTATGATACg                                                                                                                          >  1:136677/1‑141 (MQ=255)
                tGCGTCATCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGt                                                                                                                         >  1:553440/1‑141 (MQ=255)
                agCGTCAACGATAAGCAATAAAAAAGACGGTTTTGTAATAGAATTTGTTCCCGAAAGCCCTCTTGTAATATTTGGAAAATCAAAGAACATAATGTCTGTACCTGCCGAACCTAAAAAAAAAAAAAAAAAAGTATGATACGt                                                                                                                         <  2:339527/140‑1 (MQ=255)
                 gCGTCATCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAATGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGta                                                                                                                        <  1:349858/141‑1 (MQ=255)
                 gCGTCATCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGta                                                                                                                        >  1:233435/1‑141 (MQ=255)
                  cGTCATCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGtat                                                                                                                       <  1:330484/141‑1 (MQ=255)
                  cGTCATCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGtat                                                                                                                       >  2:567137/1‑141 (MQ=255)
                    tCATCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATaa                                                                                                                     >  2:466396/1‑141 (MQ=255)
                       tCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTCCATCATCTGCAAAAAAAGTATGATACGTATAAata                                                                                                                  >  1:205592/1‑141 (MQ=255)
                       tCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAata                                                                                                                  <  1:32406/141‑1 (MQ=255)
                        cGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAatat                                                                                                                 >  2:309767/1‑141 (MQ=255)
                           tAGGCACGCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAAGCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATATGCAAAAAAAGTATGATACGTATAAATatcat                                                                                                              <  2:486776/141‑1 (MQ=255)
                                tttctAAAAGACGGTTTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCGTCTGCAAAAAAAGTATGATACGTATAAATATCATCttga                                                                                                         <  1:378614/139‑1 (MQ=255)
                                       aataCGTTCTTGTAATCGAATTTGTTCCTACCTGCCCTTTTGTAATATTTGAAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCAACTAAAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTg                                                                                                  <  2:254231/138‑1 (MQ=255)
                                       aaGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTg                                                                                                  >  1:125577/1‑141 (MQ=255)
                                       aaGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTg                                                                                                  <  1:569677/141‑1 (MQ=255)
                                         gACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTc                                                                                                <  1:201741/141‑1 (MQ=255)
                                            ggTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATACTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTAAAAAAAAAGAAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCttt                                                                                             <  2:188005/141‑1 (MQ=255)
                                             gTCTTGTAATGGCATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCtttt                                                                                            >  2:539195/1‑141 (MQ=255)
                                                    gagggAAGTGGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCCGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTaa                                                                                     <  1:444060/138‑1 (MQ=255)
                                                    aagggAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTTTACCTGCCGAACCTACAGCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTaa                                                                                     <  2:232981/138‑1 (MQ=255)
                                                    aaTGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTaa                                                                                     <  2:484229/141‑1 (MQ=255)
                                                      tGGAATTTGTACCTGCCTGCCCTTGTGTAATATTTTGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAAc                                                                                   <  2:9095/141‑1 (MQ=255)
                                                       ggAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACg                                                                                  <  1:4646/141‑1 (MQ=255)
                                                           tttGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATg                                                                              >  1:86350/1‑141 (MQ=255)
                                                               ttCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGACTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTa                                                                          >  2:370160/1‑141 (MQ=255)
                                                                tCCTGCCTGCCCTTTTGTAATCTTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATGCGTATAAATATCATCTTGATTGGCTGTCTTTTTAACTAAACGCATGCCTAg                                                                         >  1:225349/1‑141 (MQ=255)
                                                                        gccCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCCTCTGCAAAAAAAGTATGATACGTATACATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACttttt                                                                 >  1:328709/1‑138 (MQ=255)
                                                                        gccCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTg                                                                 <  2:35804/141‑1 (MQ=255)
                                                                         cccTTTTGTAATATTTGGAAAATCAAAGAAAGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGt                                                                <  1:295883/141‑1 (MQ=255)
                                                                         cccTTTGGTAATATTTGGAAAATCAAAGAAGGTCATGTCTGTACCTGACGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGt                                                                <  2:155725/141‑1 (MQ=255)
                                                                           cTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACGGTCTTTTTAACTAAAAGCATGCCTAGTACTTCTGTaa                                                              >  1:425334/1‑141 (MQ=255)
                                                                             tttGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTaaaa                                                            <  2:233435/141‑1 (MQ=255)
                                                                              ttGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCTTAGTACTTCTGTaaaaa                                                           >  2:488279/1‑141 (MQ=255)
                                                                              ttGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTaaaaa                                                           >  1:558084/1‑141 (MQ=255)
                                                                                 tAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAAtt                                                        >  2:292240/1‑141 (MQ=255)
                                                                                 tAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAAtt                                                        <  1:198070/141‑1 (MQ=255)
                                                                                   atatTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGAAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAAAGACTACCTAGAACTTCTGTCAAAAAtttt                                                      >  1:199261/1‑140 (MQ=255)
                                                                                       ttGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAAt                                                  <  2:239715/141‑1 (MQ=255)
                                                                                                 aaaGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGc                                        <  1:287927/141‑1 (MQ=255)
                                                                                                 aaaGAACGTCATGACTGTACCTGCCCAACCTACATCATATGCAAAAAAAGTATGATACGTCTAACTATCATATAGAGTGCCTGTCTTTTTATCTAAACGCCTGCCTAGTACTTCTGTAAAAAATGTATAATTACGTTCTGc                                        >  2:79316/1‑141 (MQ=255)
                                                                                                   aGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGcat                                      >  2:277898/1‑141 (MQ=255)
                                                                                                    gAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGcatc                                     >  2:376819/1‑141 (MQ=255)
                                                                                                     aaCGTCATGTCTGTACCTGCCGAACCTACATCATCTGAAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTCGTACTTCTGTAACCACGTTAGAATTACGTTATAcatcc                                    >  2:116747/1‑140 (MQ=255)
                                                                                                        ttCATGTCTGTACCTGCCGAACCTAAATAATATGAAAATAAAGTATGATACGTATGAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCt                                 <  2:126596/140‑1 (MQ=255)
                                                                                                        ctCATGTCTGTACCTGCCGCTCCGACAGCATCTTCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCt                                 <  1:528490/140‑1 (MQ=255)
                                                                                                         tCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCAGGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCTg                                <  1:190956/141‑1 (MQ=255)
                                                                                                          cATGTCTGTACCTGCCGAACCTACATCCTCTGCAAAAAAAGTATGATCCGTATAAATATCATCTTGATTTACTGTATTTTTAACTAAACGCATGCCTAGTACTTCTCTAAGAATTTCATAATTACGTTCTggacactctgt                               >  2:129216/1‑131 (MQ=255)
                                                                                                           aTGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCTGTc                              >  2:302964/1‑141 (MQ=255)
                                                                                                                  tACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCTGTCATTGCAg                       <  1:126703/141‑1 (MQ=255)
                                                                                                                      tGCCGAACCTACATCATCTGCAAAAAAAGTAGGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCTGTCATTGCAGTAAc                   <  2:220201/141‑1 (MQ=255)
                                                                                                                                   tcatcTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCTGTCATTGCAGTAACATGGTGAATACCt      >  1:171997/1‑141 (MQ=255)
                                                                                                                                   tcatcTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCTGTCATTGCAGTAAAATGGTGAATACCt      >  1:254163/1‑141 (MQ=255)
                                                                                                                                     ctcTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCTGTCATTGCAGTAACATGGTGAATACCTAg    <  1:570270/140‑1 (MQ=255)
                                                                                                                                     atcTGAAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCCTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCTGTCATTGCAGTAACATGGTGAATACCTAg    <  2:536503/141‑1 (MQ=255)
                                                                                                                                      tcTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAACCGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCTGTCATTGCAGTAACATGGTGAATACCTAGt   <  2:205592/141‑1 (MQ=255)
                                                                                                                                       cTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCTGTCATTGCAGTAACATGGTGAATACCTAGTa  <  2:86350/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
GTTCATAATATGTTAATGCGTCATCGTTAGGCACTCTAAAAGACGGTCTTGTAATGGAATTTGTTCCTGCCTGCCCTTTTGTAATATTTGGAAAATCAAAGAACGTCATGTCTGTACCTGCCGAACCTACATCATCTGCAAAAAAAGTATGATACGTATAAATATCATCTTGATTGACTGTCTTTTTAACTAAACGCATGCCTAGTACTTCTGTAAAAAATTTATAATTACGTTCTGCATCATCTGTCATTGCAGTAACATGGTGAATACCTAGTA  >  CP000730/2656211‑2656486

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