Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,422,207 A→G 7.2% K122R (AAG→AGG)  PP_3038 → TraC domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,422,2070AG 7.2% 158.2 / 2.4 55K122R (AAG→AGG) PP_3038TraC domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (29/22);  new base G (2/2);  total (31/24)
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.02e-01

CCGACGAGCGTGAAGTGATGCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTCGGCCTGCAGGTAATTTTCCCGGCTGTCCAAGAGG  >  NC_002947/3422081‑3422353
                                                                                                                              |                                                                                                                                                   
ccGACGAGCGTGAAGTGATGCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTTCCTGGAt                                                                                                                              <  4:303503/149‑1 (MQ=255)
ccGACGAGCGTGAAGTGATGCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCATAGCGCCTACCTGGAt                                                                                                                              >  1:118327/1‑149 (MQ=255)
     gaGCGTGAAGTGATGCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGa                                                                                                                           >  5:27350/1‑147 (MQ=255)
        cGTGAAGTGATGCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCa                                                                                                                      <  7:101368/149‑1 (MQ=255)
          tGAAGTGATGCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAgg                                                                                                                    >  8:44524/1‑149 (MQ=255)
            aaGTGATGCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGtt                                                                                                                  >  2:66716/1‑149 (MQ=255)
                 aTGCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCg                                                                                                              <  4:319268/148‑1 (MQ=255)
                 aTGCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCg                                                                                                              >  3:319268/1‑148 (MQ=255)
                  tGCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCACGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGt                                                                                                            >  1:77578/1‑149 (MQ=255)
                     gcgcCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTc                                                                                                         >  6:123310/1‑149 (MQ=255)
                      cgcCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCg                                                                                                        >  3:30545/1‑149 (MQ=255)
                            gCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGGCCGAAAAGGGGCGCTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTc                                                                                                   <  1:66716/149‑1 (MQ=255)
                                gCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTa                                                                                              >  3:49258/1‑149 (MQ=255)
                                    aCGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATgc                                                                                           <  8:315928/148‑1 (MQ=255)
                                      gCCAAACGGCGTGCTGCAGGGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGcc                                                                                        >  1:117575/1‑149 (MQ=255)
                                      gCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGcc                                                                                        >  2:32286/1‑149 (MQ=255)
                                        cAAACGGCGGGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGCAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTCCCTGGATCGAAAGCAGGTTGGCGGTTTCGGCGTCCGCTATGCGCCg                                                                                       >  4:220360/1‑148 (MQ=255)
                                                    tgcAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCgg                                                                          >  2:278933/1‑149 (MQ=255)
                                                        gggaTCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCt                                                                      <  2:117575/147‑1 (MQ=255)
                                                           aTCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAGGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGt                                                                   >  6:54417/1‑149 (MQ=255)
                                                           aTCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAGGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGt                                                                   >  6:49352/1‑149 (MQ=255)
                                                             cgccgcGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCgg                                                                                <  5:53261/134‑1 (MQ=255)
                                                             cgccgcGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCgg                                                                                >  6:53261/1‑134 (MQ=255)
                                                             cgccgcGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTc                                                                 >  6:234470/1‑149 (MQ=255)
                                                              gccgcGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTc                                                                 <  4:30545/148‑1 (MQ=255)
                                                               ccgcGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGcgcg                                                                           >  7:271737/1‑137 (MQ=255)
                                                               ccgcGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGcgcg                                                                           <  8:271737/137‑1 (MQ=255)
                                                               ccgcGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCg                                                               <  1:278933/149‑1 (MQ=255)
                                                                cgcgAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTAAGCGCCCAAGTAGGGCGCGGTGCTGGTTCCGa                                                              >  8:279714/1‑149 (MQ=255)
                                                                 gcgAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGGGCCGAAGTCGGGCGCGGTGCTGGTTCCGAt                                                             <  1:32286/149‑1 (MQ=255)
                                                                    aaCGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGaa                                                                                     <  4:135580/122‑1 (MQ=255)
                                                                    aaCGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGaa                                                                                     >  3:135580/1‑122 (MQ=255)
                                                                    aaCGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCa                                                          <  2:77578/149‑1 (MQ=255)
                                                                         tGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGtt                                                                  <  6:335500/136‑1 (MQ=255)
                                                                         tGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCCAAGACGGGCGCGGTGCTGGtt                                                                  >  5:335500/1‑136 (MQ=255)
                                                                          gCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAg                                                         <  8:158063/144‑1 (MQ=255)
                                                                          gCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAg                                                         >  7:158063/1‑144 (MQ=255)
                                                                              gCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAg                                                <  7:44524/149‑1 (MQ=255)
                                                                                      gCCAAACGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTCGCGCATTCGGCGTCCGCTATGCGCCGAAGACGGGCGCGGTGCAGGTACCGATGCAGAATGCTCAGGCCGca                                          >  2:178638/1‑146 (MQ=255)
                                                                                           gCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTCg                                   <  4:49258/149‑1 (MQ=255)
                                                                                             gTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTc                                    <  5:9602/146‑1 (MQ=255)
                                                                                             gTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTc                                    >  6:9602/1‑146 (MQ=255)
                                                                                             gTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTCGGc                                 <  5:234470/149‑1 (MQ=255)
                                                                                                 ggAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGCCGACAATGGCGTCCTga                             <  1:178638/149‑2 (MQ=255)
                                                                                                 ggAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTCGGCCTg                              >  7:322403/1‑148 (MQ=255)
                                                                                                 ggAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTCGGCCTg                              >  2:183345/1‑148 (MQ=255)
                                                                                                  gAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTCGGCCTg                              >  1:117553/1‑147 (MQ=255)
                                                                                                           ttCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGCCGACAATGGCGTCCTga                             <  6:104033/139‑2 (MQ=255)
                                                                                                           ttCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGCCGACAATGGCGTCCTga                             >  5:104033/1‑138 (MQ=255)
                                                                                                              aaGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGCATGCTCAGGACGCCATTCTCGTCCGGCAGGTAATTTAccc                >  3:190277/1‑149 (MQ=255)
                                                                                                              aaGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTCGGCCTGCAGGTAATTTTcc                 <  2:142725/148‑1 (MQ=255)
                                                                                                              aaGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTCGGCCTGCAGGTAATTTTcc                 >  1:142725/1‑148 (MQ=255)
                                                                                                                  gTATGCCCGAAAGGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTCGGCCTGCAGGTAATTTTCCCGGCt            <  5:54417/149‑1 (MQ=255)
                                                                                                                  gTATGCCCGAAAGGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTCGGCCTGCAGGTAATTTTCCCGGCt            <  5:49352/149‑1 (MQ=255)
                                                                                                                   tATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTCGGCCTGCAGGTAATTTTCCCGGCTg           >  7:59898/1‑149 (MQ=255)
                                                                                                                            aaaGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGATGTCGGGCGCGGTGCTGGTTCCGTTGCAGAATGCTCAGGACGCCATTGTCGGCCTGCAGGTAATTTTCCCGGCTGTCCAAGAgg  <  5:79291/149‑1 (MQ=255)
                                                                                                                              |                                                                                                                                                   
CCGACGAGCGTGAAGTGATGCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCTCAGGACGCCATTGTCGGCCTGCAGGTAATTTTCCCGGCTGTCCAAGAGG  >  NC_002947/3422081‑3422353

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