Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,213,225 T→A 8.2% V17V (GTT→GTA)  PP_1956 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,213,2250TA 8.2% 177.2 / 4.9 61V17V (GTT→GTA) PP_1956hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/32);  new base A (2/3);  total (26/35)
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.29e-01

TCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTCCTACCTACTGGACAGAG  >  NC_002947/2213084‑2213372
                                                                                                                                             |                                                                                                                                                   
tctTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATTTGGTTTTTgg                                                                                                                                                <  5:76337/147‑1 (MQ=255)
  tttGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAAt                                                                                                                                            <  4:238057/149‑1 (MQ=255)
    tGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGaaa                                                                                                                                             <  5:119120/146‑1 (MQ=255)
      cAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATAcc                                                                                                                                         <  5:330765/148‑1 (MQ=255)
         gTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACAACAt                                                                                                                                     >  6:279209/1‑149 (MQ=255)
           aaTCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACAtt                                                                                                                                    <  2:55242/148‑1 (MQ=255)
            aTCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTtg                                                                                                                                  <  4:238291/149‑1 (MQ=255)
             tCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTATTTCCAAAAACCACAtgc                                                                                                                                   <  1:227193/147‑3 (MQ=255)
             tCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTATTTCCAAAAACCACAtgc                                                                                                                                   >  2:227193/1‑145 (MQ=255)
                ggTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTtgctgc                                                                                                                              >  6:139252/1‑149 (MQ=255)
                 gTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGcc                                                                                                                             >  3:261929/1‑149 (MQ=255)
                 gTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGcc                                                                                                                             >  4:128468/1‑149 (MQ=255)
                 gTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTATTTCCAAAAACCACAtgcag                                                                                                                                 <  7:272388/145‑5 (MQ=255)
                 gTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTATTTCCAAAAACCACAtgcag                                                                                                                                 <  5:208999/145‑5 (MQ=255)
                     aaaCCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGaaa                                                                                                                         <  5:139252/149‑1 (MQ=255)
                      aaCCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGaaaa                                                                                                                        <  6:128929/149‑1 (MQ=255)
                       aCCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGaaaa                                                                                                                        >  6:328907/1‑148 (MQ=255)
                         cGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAatat                                                                                                                     <  7:43992/149‑1 (MQ=255)
                           gAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCa                                                                                                                   <  1:199518/149‑1 (MQ=255)
                            accAACCGCCTCCAAAAATAAAAAATTCCTGCCTCAAAAGGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCaa                                                                                                                  <  1:111099/147‑1 (MQ=255)
                                aCCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAagcca                                                                                                              >  8:106317/1‑149 (MQ=255)
                                      tCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGc                                                                                                        >  6:207735/1‑149 (MQ=255)
                                        cAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGccc                                                                                                      <  4:261929/149‑1 (MQ=255)
                                         aataatAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGccc                                                                                                      >  1:215907/1‑148 (MQ=255)
                                              taaCAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTgg                                                                                                >  7:93297/1‑149 (MQ=255)
                                                aCAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTggcg                                                                                              <  5:207735/149‑1 (MQ=255)
                                                 cAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAAtt                                                                                                  <  6:113089/144‑1 (MQ=255)
                                                 cAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAAtt                                                                                                  >  5:113089/1‑144 (MQ=255)
                                                  aaTATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTggcgg                                                                                             >  1:245427/1‑148 (MQ=255)
                                                    tatGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCAc                                                                                          <  3:128468/149‑1 (MQ=255)
                                                           ccGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGata                                                                                     <  6:327092/147‑1 (MQ=255)
                                                           ccGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATAtt                                                                                   >  6:154838/1‑149 (MQ=255)
                                                            cGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTc                                                                                  <  6:115353/149‑1 (MQ=255)
                                                               cATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGAc                                                                               <  5:301074/149‑1 (MQ=255)
                                                                aTTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACa                                                                              <  5:154838/149‑1 (MQ=255)
                                                                     cgaagcgaaTGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGa                                                                         <  8:93297/149‑1 (MQ=255)
                                                                      gaagcgaaTGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTgg                                                                                                >  8:262570/1‑125 (MQ=255)
                                                                      gaagcgaaTGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTgg                                                                                                <  7:262570/125‑1 (MQ=255)
                                                                                 tGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCt                                                                     <  5:175207/141‑1 (MQ=255)
                                                                                 tGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCt                                                                     >  6:175207/1‑141 (MQ=255)
                                                                                 tGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTccac                                                             >  4:71176/1‑149 (MQ=255)
                                                                                        tttATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGTGGCACTGATATTCGATAAGCGATTCTACCTCCACCAGgtgt                                                      >  6:189654/1‑149 (MQ=255)
                                                                                               ataataAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATgg                                                <  5:328907/148‑1 (MQ=255)
                                                                                                   taaGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTgg                                           >  4:278829/1‑149 (MQ=255)
                                                                                                         cgcacgcTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTg                                     <  2:215907/149‑1 (MQ=255)
                                                                                                               cTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAg                               <  7:106317/149‑1 (MQ=255)
                                                                                                                tATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTgg                                           >  7:216600/1‑136 (MQ=255)
                                                                                                                tATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTgg                                           <  8:216600/136‑1 (MQ=255)
                                                                                                                    tcGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCtt                          >  6:309300/1‑149 (MQ=255)
                                                                                                                      gcgGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTg                         >  7:201926/1‑148 (MQ=255)
                                                                                                                            aTTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTggc                                                                                               >  7:146140/1‑72 (MQ=255)
                                                                                                                            aTTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTggc                                                                                               <  8:146140/72‑1 (MQ=255)
                                                                                                                             ttGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGTGGCACTGATATTCGATAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTcct                 <  5:189654/149‑1 (MQ=255)
                                                                                                                                cagcaaaTGTGGTATTTCCAAAAACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTcctacc              >  7:180202/7‑149 (MQ=255)
                                                                                                                                 ccTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTcctacc              >  1:81674/1‑148 (MQ=255)
                                                                                                                                    gCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGc                                                                           <  6:114772/84‑1 (MQ=255)
                                                                                                                                    gCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGc                                                                           >  5:327446/1‑84 (MQ=255)
                                                                                                                                    gCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGc                                                                           <  6:327446/84‑1 (MQ=255)
                                                                                                                                    gCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGc                                                                           >  5:114772/1‑84 (MQ=255)
                                                                                                                                     cATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTCctacctac           <  5:279209/147‑1 (MQ=255)
                                                                                                                                      aTGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGTCCCAATTGCCGGCACTGATCTTCGACAAGCGATTCTACCTCCACCAGTTGATTATGGCCTGGGCGTTGTTGTAGGGCTTGGATTTCCTACCTACTg         <  8:97069/148‑1 (MQ=255)
                                                                                                                                            tttttGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTCCTACCTACTGGACagag  <  3:71176/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTCCTACCTACTGGACAGAG  >  NC_002947/2213084‑2213372

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