Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 621,600 G→A 7.0% A244T (GCC→ACC)  PP_0535 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947621,6000GA 7.0% 163.0 / 2.4 58A244T (GCC→ACC) PP_0535ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base A (2/2);  total (29/29)
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.44e-01

TTCAGGGCCACATCGAAGGCAGCTACGGCACCTGGGACACCAACCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTCCGGCGAGGTGCCCGACTTCAG  >  NC_002947/621454‑621748
                                                                                                                                                   |                                                                                                                                                    
ttCAGGGCCACATCGAAGGCAGCTACGGCACCTGGGACACCAACCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGc                                                                                                                                                     >  3:300399/1‑148 (MQ=255)
     ggCCACATCGAAGGCAGCTACGGCACCTGGGACACCAACCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAacac                                                                                                                                               <  4:300399/149‑1 (MQ=255)
       ccACATCGAAGGCAGCTACGGCACCTGGGACACCAACCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCa                                                                                                                                             >  5:67598/1‑149 (MQ=255)
         acaTCGAAGGCAGCTACGGCACCTGGGACACCAACCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGc                                                                                                                                           <  4:134214/149‑1 (MQ=255)
           aTCGAAGGCAGCTACGGCACCTGCGACACCAACCAAGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGAGGACTACAACGCCAACAccacca                                                                                                                                          >  2:298643/1‑148 (MQ=255)
                     gCTACGGCACCTGGGACACCAACCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTAc                                                                                                                               <  5:167900/149‑1 (MQ=255)
                        aCGGCACCTGGGACACCAACCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGt                                                                                                                            >  3:127869/1‑149 (MQ=255)
                             aCCTGGGACACCAACCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCt                                                                                                                       >  3:169037/1‑149 (MQ=255)
                                  ggACACCAATCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGTGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGa                                                                                                                    <  1:105248/147‑1 (MQ=255)
                                  ggACACCAACCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGat                                                                                                                   <  5:89749/148‑1 (MQ=255)
                                  ggACACCAACCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTAACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGat                                                                                                                   <  6:228436/148‑1 (MQ=255)
                                   gACACCAACCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGatat                                                                                                                 <  4:127869/149‑1 (MQ=255)
                                      accaaccaTGAATAAATCGTACTTCTCGGCCCGCTGGAAGCCGGCGGTCGCGTGCGGGGGCGGTTCATTGCATCCCGTGACGACACCAATTGCCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCg                                                                                                               <  1:238480/149‑1 (MQ=255)
                                       ccaaccaTGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGTATATCGAc                                                                                                             <  6:179869/149‑1 (MQ=255)
                                        caaccaTGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGaac                                                                                                            >  4:186880/1‑147 (MQ=255)
                                          accaTGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCt                                                                                                                       <  4:77642/136‑1 (MQ=255)
                                          accaTGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCt                                                                                                                       >  3:77642/1‑136 (MQ=255)
                                            caTGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTgg                                                                                                                     <  6:107027/136‑1 (MQ=255)
                                            caTGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTgg                                                                                                                     >  5:107027/1‑136 (MQ=255)
                                            caTGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGga                                                                                                        >  7:199097/1‑149 (MQ=255)
                                             aTGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACACCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGgac                                                                                                       >  8:1827/1‑149 (MQ=255)
                                             aTGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACACCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGgac                                                                                                       >  8:6344/1‑149 (MQ=255)
                                              tGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGgacg                                                                                                      >  7:200173/1‑149 (MQ=255)
                                                 agaaATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGacgacg                                                                                                   >  3:191944/1‑149 (MQ=255)
                                                  gaaATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGc                                                                                                  >  1:289104/1‑149 (MQ=255)
                                                        gtccTGAGCGTCCCGCTGGACGCCGCCGGTCGCTTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTa                                                                                             <  3:229191/146‑1 (MQ=255)
                                                                cGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGt                                                                                    >  4:130080/1‑149 (MQ=255)
                                                                     cGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGc                                                                               <  6:67598/149‑1 (MQ=255)
                                                                      gCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGcc                                                                              <  4:169037/149‑1 (MQ=255)
                                                                               cgccgGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACg                                                                     <  4:278989/149‑1 (MQ=255)
                                                                                     tCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTga                                                               <  2:289104/149‑1 (MQ=255)
                                                                                       gcgTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTgaag                                                             <  8:199097/149‑1 (MQ=255)
                                                                                        cgTGCGCGGGCGCTTCATTGCATCCCGTGACGACACGAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTgaaga                                                            <  1:298643/149‑1 (MQ=255)
                                                                                        cgTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTgaag                                                             >  7:118/1‑148 (MQ=255)
                                                                                         gTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTgaagaa                                                           <  8:200173/149‑1 (MQ=255)
                                                                                         gTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTgaagaa                                                           <  3:130080/149‑1 (MQ=255)
                                                                                                ggCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTg                                                                <  1:129819/137‑1 (MQ=255)
                                                                                                ggCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTg                                                                >  2:129819/1‑137 (MQ=255)
                                                                                                  cgcTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACACCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCAcccc                                                  <  7:6344/149‑1 (MQ=255)
                                                                                                  cgcTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACACCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCAcccc                                                  <  7:1827/149‑1 (MQ=255)
                                                                                                     ttCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTCAAGAACATCACCCCGca                                               >  1:305096/1‑149 (MQ=255)
                                                                                                     ttCATTGCATCCCGTGAAGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCACCAGCTATTAAGGTGCGCTGGATATCGACCTGGACGACGCGACTACGGTGACGTTAGGCCTGATGCACGAAGTGAAGAAAATCACCCCTca                                               >  6:22448/1‑149 (MQ=255)
                                                                                                            cATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGc                                                                               >  3:239826/1‑110 (MQ=255)
                                                                                                            cATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGc                                                                               <  4:239826/110‑1 (MQ=255)
                                                                                                                 cGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATACACGAAGTGAAGAACATCACCCCGCACAACGGCTAccc                                   >  1:15447/1‑149 (MQ=255)
                                                                                                                       gacACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCAc                             <  4:191944/149‑1 (MQ=255)
                                                                                                                          acCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGtt                          >  5:127189/1‑149 (MQ=255)
                                                                                                                            cAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGCTGCACGAAGTGAAGAACCTCACCCCGCACAACGCATACCCGCCCACGttc                         >  8:222336/1‑147 (MQ=255)
                                                                                                                             aa‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTc                       <  6:127189/149‑1 (MQ=255)
                                                                                                                              a‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTcc                      >  2:289615/1‑149 (MQ=255)
                                                                                                                              a‑TGGAGAAGTGGACTACAACGCCAACAACAGCTACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTcc                      <  5:14910/149‑1 (MQ=255)
                                                                                                                                     aaGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTCCGGCGAg                >  8:262080/1‑149 (MQ=255)
                                                                                                                                      aGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCTCTTTATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTCCGGCGAg                <  3:186880/148‑1 (MQ=255)
                                                                                                                                        tGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTCCGGCGAGGTg             >  2:275542/1‑149 (MQ=255)
                                                                                                                                           aCTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTCCgg                    <  5:169119/139‑1 (MQ=255)
                                                                                                                                           aCTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTCCgg                    >  6:169119/1‑139 (MQ=255)
                                                                                                                                             tACAACCCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACACGGCCACGTTGTCCTGCAAGGTGACCGa        >  6:127932/1‑149 (MQ=255)
                                                                                                                                                aaCGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTCCGGCGAGGTGCCCGACtt     >  3:256247/1‑149 (MQ=255)
                                                                                                                                                   gCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTCCGGCGAGGTGCCCGACTTCAg  <  8:118/149‑1 (MQ=255)
                                                                                                                                                   gCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTCCGGCGAGGTGCCCGACTTCAg  <  7:262080/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTCAGGGCCACATCGAAGGCAGCTACGGCACCTGGGACACCAACCATGAAGAAATCGACCTGAGCGGCCCGCTGGACGCCGCCGGTCGCGTGCGCGGGCGCTTCATTGCATCCCGTGACGACACCAA‑TGGCGAAGTGGACTACAACGCCAACACCAGCAACAGCTATTACGGTGCGCTGGATATCGACCTGGACGACGCCACTACGCTGTCGTTCGGCCTGATGCACGAAGTGAAGAACATCACCCCGCACAACGGCTACCCGGCCACGTTGTCCGGCGAGGTGCCCGACTTCAG  >  NC_002947/621454‑621748

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: ‑