Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 76,334 G→A G322G (GGC→GGT sgrR ← DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091376,3340GA97.1% 100.6 / ‑5.6 35G322G (GGC→GGTsgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/10);  minor base C (0/1);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAACAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCGCCCTGTGGACCTTTTAGCATCAGCCCTCCGGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGAATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAGTAAC  >  NC_000913/76100‑76580
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
gCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTa                                                                                                                                                                                                                                        <  1:248266/251‑1 (MQ=255)
       ggTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATc                                                                                                                                                                                                                                   >  2:328478/1‑249 (MQ=255)
              cGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAGTAATAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCCCCCTGTGGACCTTTTAGCATCAGACCTcca                                                                                                                                                                                                                          <  1:29876/251‑2 (MQ=255)
              cGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTcca                                                                                                                                                                                                                          <  1:285608/251‑2 (MQ=255)
               gCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTccag                                                                                                                                                                                                                         <  1:284359/251‑3 (MQ=255)
                           gCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACTTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCtcgtcg                                                                                                                                                                                                              <  2:190445/250‑1 (MQ=255)
                                         aaCGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGc                                                                                                                                                                                                  >  2:168278/1‑248 (MQ=255)
                                                 aGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCTCCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGAc                                                                                                                                                                                       <  2:161583/251‑1 (MQ=255)
                                                     gTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGAAGtt                                                                                                                                                                                    >  1:161066/1‑250 (MQ=255)
                                                          tAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGCTTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCGCACCCTGTGGACCTTTTAGCATCAGACCTCCCGCTGGCTCGTCGGCAATTTGCGGCAGAACCCAGACGTTAACTtc                                                                                                                                                                              >  2:91931/1‑251 (MQ=255)
                                                               tgggAACCACAGTTGCTGGTACTGTTCCTCAGCGAACTAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTGGGCAATTTCCGGCAGAACCCAGCCGTTAACTtcgtcg                                                                                                                                                                          >  2:93273/2‑250 (MQ=255)
                                                                                ggTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCg                                                                                                                                                        >  2:266105/1‑251 (MQ=255)
                                                                                                gAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGTTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCCAATg                                                                                                                                        >  2:86124/1‑251 (MQ=255)
                                                                                                    tAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTg                                                                                                                                    >  1:48965/1‑251 (MQ=255)
                                                                                                                   ggAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGAtt                                                                                                                     >  1:253515/1‑251 (MQ=255)
                                                                                                                            acaTAGCTTACCCAGACCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGtt                                                                                                            >  1:268919/1‑251 (MQ=255)
                                                                                                                                                          ttCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCgat                                                                              >  2:293626/1‑251 (MQ=255)
                                                                                                                                                            cGCCCCGCGATGGGTGAGGCTGTCGAATCGTAAATAGTAGCAACCTTCCTCCAGGCGGATTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCgatgg                                                                            >  2:39901/1‑251 (MQ=255)
                                                                                                                                                                       gggTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGAtggctggc                                                                  <  1:205588/250‑1 (MQ=255)
                                                                                                                                                                             ggCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGtt                                                            >  1:10392/1‑250 (MQ=255)
                                                                                                                                                                                       taGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGttt                                                 >  2:22555/2‑251 (MQ=255)
                                                                                                                                                                                        aGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATg                                                                                                                                        >  2:304695/1‑163 (MQ=255)
                                                                                                                                                                                        aGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATg                                                                                                                                        <  1:304695/163‑1 (MQ=255)
                                                                                                                                                                                           aaaTAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTGATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGGTTCCCa                                             <  2:252920/251‑1 (MQ=255)
                                                                                                                                                                                             aTAGGAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCAtt                                           >  2:90978/1‑251 (MQ=255)
                                                                                                                                                                                               agtagCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTc                                          <  2:221560/250‑1 (MQ=255)
                                                                                                                                                                                                               cAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCg                         >  1:245332/1‑251 (MQ=255)
                                                                                                                                                                                                                aGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGc                        >  2:316577/1‑251 (MQ=255)
                                                                                                                                                                                                                          cAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGc              >  1:43401/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGCCCTCCAGCTGGTTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGccc            >  2:39463/1‑251 (MQ=255)
                                                                                                                                                                                                                             tctcTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGcccc           >  1:27453/1‑251 (MQ=255)
                                                                                                                                                                                                                                ctTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCcagc        >  1:204672/1‑251 (MQ=255)
                                                                                                                                                                                                                                ctTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCcagc        >  1:147930/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ttttCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCcagca       <  1:352205/250‑1 (MQ=255)
                                                                                                                                                                                                                                      ccTCACCCTGTGGTACTCTTAGCATCAGGCCTCCAGCTGGCTCGTCGGCAATTTCCGGCCGAACCCAGATGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACAGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTCCTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCACCAGCAATAAc  <  2:196393/251‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
GCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAACAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCGCCCTGTGGACCTTTTAGCATCAGCCCTCCGGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGAATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAGTAAC  >  NC_000913/76100‑76580

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