Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,936,285 C→G 8.3% A24G (GCC→GGC)  PP_3469 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,936,2850CG 8.3% 176.1 / 6.1 60A24G (GCC→GGC) PP_3469hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/29);  new base G (2/3);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GAGTCAACTGTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAACACCTCCTTACCGCAGGCTTA  >  NC_002947/3936138‑3936432
                                                                                                                                                   |                                                                                                                                                    
gagTCAACTGTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGcc                                                                                                                                                     <  6:144787/149‑1 (MQ=255)
         gTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACg                                                                                                                                            >  4:67413/1‑149 (MQ=255)
         gTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACg                                                                                                                                            >  4:359776/1‑149 (MQ=255)
          tCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGGCAGCGCttgt                                                                                                                                            <  5:67877/148‑4 (MQ=255)
                cTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTcc                                                                                                                                     >  6:128707/1‑149 (MQ=255)
                  aTCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAgc                                                                                                                                                  <  5:192959/134‑1 (MQ=255)
                  aTCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAgc                                                                                                                                                  >  6:192959/1‑134 (MQ=255)
                      gTCTCCCAAGAGATATTTCTCGACCATACTCCAGACAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATTCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACa                                                                                                                              <  4:52672/149‑1 (MQ=255)
                      gTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGGGCTACGTTCTTCCGTGACa                                                                                                                              <  4:57654/149‑1 (MQ=255)
                      gTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACa                                                                                                                              >  7:297262/1‑149 (MQ=255)
                       tCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAg                                                                                                                             >  4:374752/1‑149 (MQ=255)
                           ccAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTAAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTcc                                                                                                                         <  2:124055/149‑1 (MQ=255)
                            cAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAgc                                                                                                                                                  >  8:341213/1‑124 (MQ=255)
                            cAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAgc                                                                                                                                                  <  7:341213/124‑1 (MQ=255)
                            cAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCg                                                                                                                        >  2:321839/1‑149 (MQ=255)
                              agagATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCATCCCGCAGCATCTACGACACTTCTACATGCATAAGAACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCAGTGAAAGCTcag                                                                                                                        >  5:144755/1‑145 (MQ=255)
                              agagATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCgta                                                                                                                      <  3:270397/149‑1 (MQ=255)
                                    atTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTg                                                                                                                >  5:355594/1‑149 (MQ=255)
                                     tttCTCGACCATACTCCATCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTg                                                                                                                <  3:257454/148‑1 (MQ=255)
                                     tttCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTgg                                                                                                               <  6:355594/149‑1 (MQ=255)
                                        ctcGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGcc                                                                                                             <  5:128707/148‑1 (MQ=255)
                                           gACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGaa                                                                                                          <  4:93729/148‑1 (MQ=255)
                                               aTACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAattatt                                                                                                     <  7:23271/149‑1 (MQ=255)
                                                     ccgCCAAAACCTCGGTAGATCTCAATGACCAGGCCGCAGCATGTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATtgtgt                                                                                                 <  3:195471/145‑1 (MQ=255)
                                                     cAGCCAACACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATtgtgt                                                                                                 <  4:293034/147‑1 (MQ=255)
                                                     cAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATtgtgt                                                                                                 <  8:221610/147‑1 (MQ=255)
                                                       gCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAg                                                                                             >  5:84436/1‑149 (MQ=255)
                                                            aaCCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAAt                                                                                        <  2:239982/149‑1 (MQ=255)
                                                                  ggTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACCAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGa                                                                                    <  1:270033/148‑1 (MQ=255)
                                                                  ggTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGa                                                                                    <  3:67413/147‑1 (MQ=255)
                                                                       aTCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCgcc                                                                             >  4:106766/1‑147 (MQ=255)
                                                                            aaTGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCa                                                                           <  5:36085/146‑1 (MQ=255)
                                                                            aaTGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGc                                                                        >  3:187074/1‑149 (MQ=255)
                                                                                  cAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAg                                                                  <  8:297262/149‑1 (MQ=255)
                                                                                  cAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAg                                                                  >  6:73197/1‑149 (MQ=255)
                                                                                   aGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAgg                                                                 >  3:309072/1‑149 (MQ=255)
                                                                                    gCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGc                                                                        >  3:84413/1‑141 (MQ=255)
                                                                                    gCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGc                                                                        <  4:84413/141‑1 (MQ=255)
                                                                                       cGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGctct                                                             <  7:71683/149‑1 (MQ=255)
                                                                                                        ttCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACa                                                                                  <  5:208297/111‑1 (MQ=255)
                                                                                                        ttCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACa                                                                                  >  6:208297/1‑111 (MQ=255)
                                                                                                        ttCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATcccc                                            >  2:270062/1‑149 (MQ=255)
                                                                                                          cTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGACAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCAccc                                           <  3:359776/148‑1 (MQ=255)
                                                                                                                gCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAAccc                                    >  8:143658/1‑149 (MQ=255)
                                                                                                                gCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGAGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAAccc                                    >  8:217185/1‑149 (MQ=255)
                                                                                                                 cATACGCACTTACACTCGGTGCTGACAATCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAAccc                                    <  6:144755/148‑1 (MQ=255)
                                                                                                                   tACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCa                                 <  4:309072/149‑1 (MQ=255)
                                                                                                                          ttACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACgg                           >  5:160122/1‑148 (MQ=255)
                                                                                                                            acacTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACACATAACCCGCAGTACGGaa                         >  8:323132/1‑148 (MQ=255)
                                                                                                                              acTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAac                      <  3:374752/149‑1 (MQ=255)
                                                                                                                                  ggTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAACAcct                  <  7:288545/149‑1 (MQ=255)
                                                                                                                                         aCAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATc                                                                                       <  8:88975/73‑1 (MQ=255)
                                                                                                                                         aCAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATc                                                                                       >  7:88975/1‑73 (MQ=255)
                                                                                                                                          gtaGCGCTGGCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGa                                                           >  3:116188/3‑100 (MQ=255)
                                                                                                                                          gtaGCGCTGGCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGa                                                           <  4:116188/98‑1 (MQ=255)
                                                                                                                                           aaGCGCTGCCAGCGCTACGTTCTTCCGTGACAGATCCGTAGTACTGGCCGAATTATTGTGTGACGCCAATCAGCCAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACACGAAGTACGGAAAACACCTCCTTACCg          >  1:236869/1‑148 (MQ=255)
                                                                                                                                             gcgcTGGCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCAc                                             <  8:142236/111‑1 (MQ=255)
                                                                                                                                             gcgcTGGCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCAc                                             >  7:142236/1‑111 (MQ=255)
                                                                                                                                             gcgcTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAACACCTCCTTACCGCAg       <  5:73197/149‑1 (MQ=255)
                                                                                                                                               gcTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAACACCTCCTTACCGCAGGc     <  1:321839/149‑1 (MQ=255)
                                                                                                                                                cTGCCAGCGCTACGTTCTTCCGTGACCGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGAGAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGCAAAACACCTCCTTACCGCAGGCt    >  1:316391/1‑149 (MQ=255)
                                                                                                                                                   ccAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAACACCTCCTTACCGCAGGCTTa  >  2:25832/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAGTCAACTGTCAAATCTATCAGTCGCCCAAGAGATATTTCTCGACCATACTCCAGCCAAAACCTCGGTAGATCTCAATGACCAGCCCGCAGCATCTACGACACTTCTACATGCATACGCACTTACACTCGGTGCTGACAAGCGCTGCCAGCGCTACGTTCTTCCGTGACAGCTCCGTAGTACTGGCCGAATTATTGTGTGAAGCCAATCAGACAGCGACCACGCTTGTAGGCTCTGACATGGCCAGATCACCCATAACCCGCAGTACGGAAAACACCTCCTTACCGCAGGCTTA  >  NC_002947/3936138‑3936432

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: