Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,050,724 C→T 7.0% intergenic (‑57/‑234) gcvT ← / → yqfI aminomethyltransferase/protein YqfI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,050,7240CT 7.0% 175.9 / 2.6 57intergenic (‑57/‑234)gcvT/yqfIaminomethyltransferase/protein YqfI
Reads supporting (aligned to +/- strand):  ref base C (37/16);  new base T (2/2);  total (39/18)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGCGACAGAATA  >  NC_000913/3050578‑3050871
                                                                                                                                                  |                                                                                                                                                   
gTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGTTAATTTCACCATGaa                                                                                                                                                  <  1:1159411/150‑1 (MQ=255)
gTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGaa                                                                                                                                                  >  1:1121320/1‑150 (MQ=255)
  aaTGCAGCGGCATCATCCAGCCGGGGAACTCCACCAGGCGAGCGCCGCAAAGCGGGTGGTGGTCGTACAAAGGAGTCTGTTGGGCCAGCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAAc                                                                                                                                                <  1:938911/150‑1 (MQ=255)
    tGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCGTGGTGATAATTTCACCATGaaaaag                                                                                                                                              <  1:1112466/150‑4 (MQ=255)
               catcCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGc                                                                                                                                   <  2:99884/150‑1 (MQ=255)
                        gTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCgg                                                                                                                          >  2:214332/1‑150 (MQ=255)
                          ggAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTc                                                                                                                        >  1:601876/1‑150 (MQ=255)
                                ccaccaTGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAAc                                                                                                                  >  1:1139475/1‑150 (MQ=255)
                                        cGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACg                                                                                                                 <  2:310318/143‑1 (MQ=255)
                                        cGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACg                                                                                                                 >  1:310318/1‑143 (MQ=255)
                                          aGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTtc                                                                                                        <  2:1206432/150‑1 (MQ=255)
                                             gccgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTctct                                                                                                     >  1:131719/1‑150 (MQ=255)
                                              ccgcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCtt                                                                                                    >  2:239350/1‑150 (MQ=255)
                                                gcAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCtttt                                                                                                  >  1:537830/1‑150 (MQ=255)
                                                  aaaGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGc                                                                                                <  1:1207103/150‑1 (MQ=255)
                                                          tgttgtTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACAACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACtt                                                                                        >  2:482076/1‑150 (MQ=255)
                                                              gttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTacca                                                                                    >  2:1228664/1‑150 (MQ=255)
                                                              gttCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTacca                                                                                    >  1:747785/1‑150 (MQ=255)
                                                                       aaGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAg                                                                           >  2:268602/1‑150 (MQ=255)
                                                                            gTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGt                                                                      <  1:214332/150‑1 (MQ=255)
                                                                            gTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGt                                                                      >  2:301718/1‑150 (MQ=255)
                                                                              cTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTa                                                                    <  2:799070/150‑1 (MQ=255)
                                                                                gttgtGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAAc                                                                  <  1:423726/150‑1 (MQ=255)
                                                                                     gCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATaa                                                             >  1:888125/1‑150 (MQ=255)
                                                                                       cATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAgg                                                           <  2:1021664/150‑1 (MQ=255)
                                                                                             gTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTGTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAAt                                                     >  2:617650/1‑150 (MQ=255)
                                                                                               ccTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTg                                                   >  2:852954/1‑150 (MQ=255)
                                                                                                cTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGa                                                  >  1:60210/1‑150 (MQ=255)
                                                                                                 tCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGat                                                 >  2:972544/1‑150 (MQ=255)
                                                                                                  cATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatc                                                >  1:1313177/1‑150 (MQ=255)
                                                                                                   aTTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatca                                               <  2:888125/150‑1 (MQ=255)
                                                                                                     tGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGatgatc                                             >  1:619801/1‑150 (MQ=255)
                                                                                                       aaTAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATcac                                           <  1:268602/150‑1 (MQ=255)
                                                                                                       aaTAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATcac                                           >  2:114793/1‑150 (MQ=255)
                                                                                                         tAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTtctc                                                                                                      >  1:1187523/1‑89 (MQ=255)
                                                                                                         tAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTtctc                                                                                                      <  2:1187523/89‑1 (MQ=255)
                                                                                                            ggggggCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCGACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTAGCACCGAAACTGACTGTTAACGATAAGGTAATTTTGATCATCACATTAg                                      <  2:1315906/148‑1 (MQ=255)
                                                                                                              ggggCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCt                                    <  1:239350/150‑1 (MQ=255)
                                                                                                                  cTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATg                                >  1:36349/1‑150 (MQ=255)
                                                                                                                  cTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATg                                >  1:1224524/1‑150 (MQ=255)
                                                                                                                    gACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGt                              >  1:599861/1‑150 (MQ=255)
                                                                                                                      cAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTa                            >  2:571991/1‑150 (MQ=255)
                                                                                                                      cAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCATATGGTTa                            >  1:1257805/1‑150 (MQ=255)
                                                                                                                       aaCTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaa                           >  1:367351/1‑150 (MQ=255)
                                                                                                                        aCTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGt                                                          <  2:358126/118‑1 (MQ=255)
                                                                                                                        aCTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGt                                                          >  1:358126/1‑118 (MQ=255)
                                                                                                                        aCTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaa                          >  2:1054698/1‑150 (MQ=255)
                                                                                                                         cTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTtaaa                         >  2:1063839/1‑150 (MQ=255)
                                                                                                                          tttttCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaaa                        >  1:478034/1‑150 (MQ=255)
                                                                                                                           gtttCGTGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTaaaaaa                       >  1:694745/2‑150 (MQ=255)
                                                                                                                             ttCGTGGTGATAATTTCACCATGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATg                     >  1:699321/1‑150 (MQ=255)
                                                                                                                             ttCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATg                     >  1:25418/1‑150 (MQ=255)
                                                                                                                                  ggTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaaaa                >  1:863666/1‑150 (MQ=255)
                                                                                                                                   gTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCaaaaa               >  1:1070195/1‑150 (MQ=255)
                                                                                                                                       aaTTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATcgc           <  2:848297/150‑1 (MQ=255)
                                                                                                                                         ttATCACCACGAAACCCAGCATCGGAGCCACTCCCGGCCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGCGa         >  1:1329791/1‑150 (MQ=255)
                                                                                                                                           aTCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGCGACa       >  1:1403481/1‑150 (MQ=255)
                                                                                                                                               ccacGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGCGACAGAAt   >  1:915577/1‑150 (MQ=255)
                                                                                                                                                cacGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGCGACAGAATa  >  2:431126/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTAATGCAGCGGCATCATCCAGCCGTGGAAATCCACCATGCGAGCGCCGCAAAGCGTGTGTTGTTCGTACAAAGGAGTCTGTTGTGCCATCTTGTCCTCATTGAATAAGCGGGGCTGACAACTTTTTCATGGTGAAATTATCACCACGAAACCCAGCATCGGAGCCACTCCCGGTCCCCAACGCAATCGTTCTCTTTTGCCTGAACTTACCACCGAAACAGACTGTTAACCATAAGGTAAAATTGATCATCACATTAGCTTATGGTTAAAAAATGCAAAAATCGCGACAGAATA  >  NC_000913/3050578‑3050871

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