Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,132,403 G→T 6.3% I51I (ATC→ATA wcaB ← putative acyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,4030GT 6.3% 233.4 / 2.8 63I51I (ATC→ATAwcaBputative acyl transferase
Reads supporting (aligned to +/- strand):  ref base G (23/36);  new base T (2/2);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

AGTGACGCCGTGGCGAATGGTAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGC  >  NC_000913/2132265‑2132539
                                                                                                                                          |                                                                                                                                        
aGTGACGCCGTGGCGAATGGTAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGa                                                                                                                                         >  1:1088915/1‑140 (MQ=255)
aGTGACGCCGTGGCGAATGGTAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGa                                                                                                                                         >  1:1579646/1‑140 (MQ=255)
    aCGCCGTGGCGAATGGTAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAAt                                                                                                                                     >  1:1138251/1‑140 (MQ=255)
    aCGCCGTGGCGAATGGTAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAAt                                                                                                                                     >  1:754815/1‑140 (MQ=255)
     cGCCGTGGCGAATGGTAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATg                                                                                                                                    <  1:340620/140‑1 (MQ=255)
     cGCCGTGGCGAATGGTAAAATCATCCCCCGCCACCCCGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAACAGCATTCGGTGATAATg                                                                                                                                    <  1:340619/140‑1 (MQ=255)
           ggCGAATGGTAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATAc                                                                                                                              <  2:674048/140‑1 (MQ=255)
                 tGGTAAAATCATCCCCCGCTACCTCGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCAcc                                                                                                                        <  1:1738634/140‑1 (MQ=255)
                 tGGTAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCAcc                                                                                                                        <  1:819882/140‑1 (MQ=255)
                 tGGTAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCAcc                                                                                                                        <  1:1370014/140‑1 (MQ=255)
                 tGGTAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCAcc                                                                                                                        <  1:1487544/140‑1 (MQ=255)
                 tGGTAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCAcc                                                                                                                        <  1:1896577/140‑1 (MQ=255)
                    tAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACcagc                                                                                                                     >  2:1651729/1‑140 (MQ=255)
                    tAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACcagc                                                                                                                     >  2:1651675/1‑140 (MQ=255)
                      aaaTCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCTTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACcagcag                                                                                                                   >  1:494025/1‑140 (MQ=255)
                      aaaTCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACcagcag                                                                                                                   >  2:1167776/1‑140 (MQ=255)
                      aaaTCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGAGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATATCAAAACCGAAAAAGCATACGGTGATAATGCGATACAGCACcagcag                                                                                                                   >  1:494013/1‑140 (MQ=255)
                            tcCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGgc                                                                                                             <  1:1773771/140‑1 (MQ=255)
                            tcCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGgc                                                                                                             >  1:1617696/1‑140 (MQ=255)
                            tcCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGgc                                                                                                             >  1:155953/1‑140 (MQ=255)
                            tcCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGgc                                                                                                             >  2:740688/1‑140 (MQ=255)
                                    tACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAACGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCcaca                                                                                                     <  1:1197156/140‑1 (MQ=255)
                                                        cgacgGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCttt                                                                                 <  2:764008/140‑1 (MQ=255)
                                                                tAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGGAAAAGCAATCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCca                                                                         <  2:928842/140‑1 (MQ=255)
                                                                  aCCGTGTTGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCcaca                                                                       <  1:155809/140‑1 (MQ=255)
                                                                  aCCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCcaca                                                                       <  1:722703/140‑1 (MQ=255)
                                                                      tGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGa                                                                   >  1:1329581/1‑140 (MQ=255)
                                                                        aTGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAAc                                                                 >  1:763577/1‑140 (MQ=255)
                                                                             tAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACaaaaa                                                            <  2:1721700/140‑1 (MQ=255)
                                                                                    gcggcgGCCAATGGTCGCGGCAGCGTGGATTTCATAACCGAAATAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCaa                                                     <  2:1590461/140‑1 (MQ=255)
                                                                                    gcggcgGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCaa                                                     <  2:318094/140‑1 (MQ=255)
                                                                                       gcggcCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTTATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACaaaa                                                             <  2:1067436/129‑1 (MQ=255)
                                                                                       gcggcCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTTATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACaaaa                                                             <  2:1067379/129‑1 (MQ=255)
                                                                                       gcggcCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTTATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACaaaa                                                             >  1:1067436/1‑129 (MQ=255)
                                                                                       gcggcCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTTATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACaaaa                                                             >  1:1067379/1‑129 (MQ=255)
                                                                                         ggcCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGa                                                <  1:537719/140‑1 (MQ=255)
                                                                                         ggcCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGa                                                <  1:507387/140‑1 (MQ=255)
                                                                                         ggcCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGa                                                <  1:1304206/140‑1 (MQ=255)
                                                                                         ggcCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGa                                                <  2:1287506/140‑1 (MQ=255)
                                                                                          gcCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGAt                                               >  1:1761751/1‑140 (MQ=255)
                                                                                          gcCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGAt                                               >  1:1761788/1‑140 (MQ=255)
                                                                                          gcCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGAt                                               >  1:347090/1‑140 (MQ=255)
                                                                                          gcCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGAt                                               >  1:649827/1‑140 (MQ=255)
                                                                                              aTGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGc                                           <  1:428317/140‑1 (MQ=255)
                                                                                                ggTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCaa                                         <  1:1635106/140‑1 (MQ=255)
                                                                                                 gTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAg                                        <  1:1770028/140‑1 (MQ=255)
                                                                                                   cgcgGCAGCCTGTATTTCAGAACCGAAACAGCATTCGGTGATAATGCGTTACAGCACCAGCAGCGGGGCCGACCACAGATTGTTGAGGACGTTCTTTTTGAGCCACACCGAACAAAAATGAGCAACACGCTAGGCAAGaa                                      <  2:156762/140‑1 (MQ=255)
                                                                                                   cgcgGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGaa                                      <  2:1572552/140‑1 (MQ=255)
                                                                                                   cgcgGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGaa                                      <  2:1560520/140‑1 (MQ=255)
                                                                                                     cgGCAGCCTGGATTTCATAACCGATAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACaa                                                               <  1:356318/113‑1 (MQ=255)
                                                                                                     cgGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACaa                                                               >  2:356318/1‑113 (MQ=255)
                                                                                                      ggCAGCCTGGATTTCCGAACCGAAAAAGCATTCGGTGATAATGCGATAGAGCACCAGCAGCGGGGCCGCCCCCAGATTGTTGAGGACGTTCTTGTGGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCa                                   <  2:1215186/140‑1 (MQ=255)
                                                                                                      ggCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCa                                   <  1:199078/140‑1 (MQ=255)
                                                                                                             tGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATgcagca                            <  2:124142/140‑1 (MQ=255)
                                                                                                             tGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATgcagca                            <  1:384668/140‑1 (MQ=255)
                                                                                                             tGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATgcagca                            <  1:1275925/140‑1 (MQ=255)
                                                                                                                     aTAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTa                    >  2:1766161/1‑140 (MQ=255)
                                                                                                                     aTAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTa                    >  2:1766235/1‑140 (MQ=255)
                                                                                                                     aTAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTa                    >  2:210000/1‑140 (MQ=255)
                                                                                                                      tAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGaa                                      >  2:574058/1‑121 (MQ=255)
                                                                                                                      tAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGaa                                      <  1:574058/121‑1 (MQ=255)
                                                                                                                                   aTTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTg      <  1:365339/140‑1 (MQ=255)
                                                                                                                                       ggTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGgcgc  <  1:973852/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                        
AGTGACGCCGTGGCGAATGGTAAAATCATCCCCCGCTACCACGTTTTTATTGATCACGACGGCGTAACCGTGATGGATAGTAAAGCGGCGGCCAATGGTCGCGGCAGCCTGGATTTCATAACCGAAAAAGCATTCGGTGATAATGCGATACAGCACCAGCAGCGGGGCCGCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGC  >  NC_000913/2132265‑2132539

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