Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,878,855 C→T A28V (GCT→GTT)  yoaK → putative membrane protein YoaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,878,8550CT100.0% 101.5 / NA 35A28V (GCT→GTT) yoaKputative membrane protein YoaK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/19);  total (16/19)

GGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCACGAAGCTGCCAAAATTGAAGCCATCGACT  >  NC_000913/1878609‑1879094
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
ggCATAGTAGAAGTCGTTCCTGAAACAACCGATGCTGACAGTGCCCTGCATGCTGCCGATACTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACATCGGTCGTATTGTTAGCATGGTTTTTTATTGTCGGCTATGTTGcc                                                                                                                                                                                                                                             <  2:426696/251‑1 (MQ=255)
  cATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCGGAATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACTTCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTCTTTTATTGGCGGCTAtgttgcgcc                                                                                                                                                                                                                                           <  2:324955/251‑4 (MQ=255)
        agAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAg                                                                                                                                                                                                                                      <  1:470481/250‑1 (MQ=255)
         gAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCa                                                                                                                                                                                                                                    <  2:364077/251‑1 (MQ=255)
                     gAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGCCGTATCCGGAGAAAAAACACAAACAGACAAAACCTTTTGTAGGGCATCCAGCGCTACATTCCTGAGGCGGATTCACATCCTTTTGATTGGTGATAACATGCGAATCGATATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGGCTAAAGAAGaaa                                                                                                                                                                                                                         >  1:267564/1‑250 (MQ=255)
                         caacCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATGAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCTCATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTGTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAaac                                                                                                                                                                                                                    <  1:395887/251‑1 (MQ=255)
                                       aGTGCCCTGCATGCTGCCGATATTGCGATGGATCAGGAGAAAAAACACAAACAGAAAACCCCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGAGAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGTCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTatga                                                                                                                                                                                                      <  2:485656/251‑1 (MQ=255)
                                       aGTGCCCTGCAGGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTatga                                                                                                                                                                                                      <  1:516213/251‑1 (MQ=255)
                                            ccTGCATGCTGCCGATATTGCGCTGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCCTAAAAATGGAAAAAACAACGATTATTATGATgggg                                                                                                                                                                                                 >  2:522241/1‑250 (MQ=255)
                                                         gATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTg                                                                                                                                                                                    <  2:458414/251‑1 (MQ=255)
                                                            atTGCGATGTATCAGGAGACAAAACACAACCAGAAAACAACTTTTGTCGCGCATCCAGCGCTACATTACCGAGGCGTATTCACATCCTTTTGATTGGCGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGATCATAAAGGTGGAAAAAACAACGATTTTTATGATggggtgggcggatattgtcg                                                                                                                                                                                 >  2:539934/1‑234 (MQ=255)
                                                                       tCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCg                                                                                                                                                                                 >  1:51691/1‑240 (MQ=255)
                                                                       tCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCg                                                                                                                                                                                 <  2:51691/240‑1 (MQ=255)
                                                                                 aaaCACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCAACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCAATGTTGCCCCGGGAGCATAAAGATGGAAAAAACAACGGTTATTATGGTGGGGGAGGCGATTATTGTCGTACTCGGCACTGAGCTggg                                                                                                                                                              >  2:353465/1‑249 (MQ=255)
                                                                                   acccaAACAGAAAACACCTTTTGTCGCGCAGCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTGTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGTCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATCATGAGGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGAtgg                                                                                                                                                          <  1:320363/248‑1 (MQ=255)
                                                                                    cacaAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGAtggg                                                                                                                                                         >  1:203509/1‑250 (MQ=255)
                                                                                                ccacCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACCACGGTTATTTTGATGGGTGTGGCGGTTATTGTCGTACTCGGCACTGAGCTGGGATGGGGGTAACTTCAcc                                                                                                                                             >  2:192896/2‑251 (MQ=255)
                                                                                                            gcgcATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCTGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAg                                                                                                                                 <  1:539866/251‑1 (MQ=255)
                                                                                                                                   aGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCt                                                                                                          <  1:130532/251‑1 (MQ=255)
                                                                                                                                    ggCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACATCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTAGTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGTTGCCTGCGTGCAGcgtt                                                                                                         <  2:182994/251‑4 (MQ=255)
                                                                                                                                     gCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGGACTCGGCACTGAGCTGGGATGGTGGTAACGTCCCCTCTAAAAAATAGCAAAAGCTGCCTGTGTGCAGCCttt                                                                                                        >  2:138903/1‑251 (MQ=255)
                                                                                                                                         aTTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGACCCGAGCACATCAAGATGAAAAAACCCACGCTTATTATAATGGGGGTGACGCTTATTGTCGTACTCGGCACTGAACTGGGCCGTTGGTAACGTCACCTCTAAAAAACAGCAAAGGCTGCCTGTGTGCAGGCTTTGTGc                                                                                                    >  1:571754/1‑251 (MQ=255)
                                                                                                                                                 ccggttGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATGTTTAGCATGGTTTTTTAATGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATTATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCGGGGATGGTGGTAACGTCGCCTCTAAAAAATAGCAAAGTCTTCCTGTGGTCAGCCTTTGTGCAAGTTAAg                                                                                            <  2:373833/247‑1 (MQ=255)
                                                                                                                                                   ttttGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGAGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGGGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCCCTGAGATGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCg                                                                                          >  1:562256/1‑251 (MQ=255)
                                                                                                                                                                  aCATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGTATGGTTTTTTATTGGCGGCTATATTGCCCCGGGAGCATAAAAATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGCACTTGGAGGGTGGTAACGTCACCTCTAAAAAATAGCAAAAGCTGCCTGTGGGCAGCCTTTGTGGACTTTAAACGCTAACTTTTacattt                                                                           >  2:447745/1‑246 (MQ=255)
                                                                                                                                                                                  tattTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCTTTAACTTTTAGTCTTTCTGTAGATAAATAg                                                            <  2:90722/250‑1 (MQ=255)
                                                                                                                                                                                     ttttCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTGTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCGTCCTGTAGATAAATAGCAc                                                         <  2:296891/250‑1 (MQ=255)
                                                                                                                                                                                                 aaTCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTGTATTGGCGGCTATGTTGCCACGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTGTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCa                                            <  2:14452/250‑1 (MQ=255)
                                                                                                                                                                                                                ggTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGGTTATTATGAGGGGTGGGGCGATTATTGTAGTACTCGGTACTGAGCTGGGCTGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGAGTGCAGCCCTTGGGCAAATTAAACATTAACATTTAAACTTCCTGTAGATAAATAGAACGACAATCGCAACAATAACGGCAACCACg                             >  1:102911/1‑251 (MQ=255)
                                                                                                                                                                                                                ggTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGGACTCGTCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCGTCCTGGAGATAAATAGCACGACAATCGCACCAATAACGGCAAcctcg                             <  2:257408/251‑4 (MQ=255)
                                                                                                                                                                                                                ggTCGTATTTTTAGCATGGTTTTGTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGCTTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGGTGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTTTGCAGCCTTTTTGCAATTTAAACGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCACg                             >  1:132757/1‑251 (MQ=255)
                                                                                                                                                                                                                     tATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTTCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGGACTCGGCACTGAGGTGGGATGGTGGTAAAGTCACCTCTTAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCCATTTAAGCGTTAAATTTTTATCTTCCTGTAGATAAATAGCCCGACAATCGCACCAATAAAGGCAACCCCGTAGCt                        >  1:132299/1‑251 (MQ=255)
                                                                                                                                                                                                                      aTTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTATAAAATAGCAAAGGCTTCCTGTGTGCATCCTTTGTGCAATTGAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCGACCACGAAGCt                        <  2:297150/250‑1 (MQ=255)
                                                                                                                                                                                                                             gCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGGGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCCCGAAGCTGCCATAAt                >  1:192750/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ttATTGGCGGCTATGTTGCCCCGGGGGCATAAAGCTGAAAAAAAAAACGATTATTATGATGGGTGTGGCGATTATTGTCGGACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAACAAAAGCTGCCTGTGGGACGCCTTTGGGCAATTTAAGCGTTAACTTTTAATCTTTCTGGAGATAAATAACACCACAATCGCACCAATAACGGCAAACACAAAGCTTCCAAAAATGAAGCCATc      >  2:15657/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ggcggcTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCACGAAGCTGCCAAAATTGAAGCCATCGACt  <  1:30179/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
GGCATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCACGAAGCTGCCAAAATTGAAGCCATCGACT  >  NC_000913/1878609‑1879094

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