Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,064,612 C→A 12.7% G129C (GGT→TGT) ‡ PP_4461 ← MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,064,6120CA12.7% 175.6 / 15.9 63G129C (GGT→TGT) ‡PP_4461MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base A (4/4);  total (30/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.69e-02

GGAATGAGCAAGAGCGCCGAGGCTACTAGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGCGGGCTTGCGCCCCC  >  NC_002947/5064467‑5064758
                                                                                                                                                 |                                                                                                                                                  
ggAATGAGCAAGAGCGCCGAGGCTACTAGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTgc                                                                                                                                                 >  2:280409/1‑149 (MQ=255)
  aaTGAGCAAGAGCGCCGAGGCTACTAGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTgcgg                                                                                                                                               <  3:102918/149‑1 (MQ=255)
      aGCAAGAGCGCCGAGGCTACTAGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGaa                                                                                                                                           >  7:239707/1‑149 (MQ=255)
         aatagCGCATAGTCTACTAGCGAGGACAGCAAAGCGCCCCAATAGTAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGATGTGCCCCATTCGCCGCAGGTGGCGAAtg                                                                                                                                         <  8:280618/145‑2 (MQ=255)
           gagCGCCGAGGCTACTAGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCGCAGGTGGCGAAtgggg                                                                                                                                      <  7:167677/149‑5 (MQ=255)
           gagCGCCGAGGCTACTAGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCtt                                                                                                                                      >  4:308098/1‑149 (MQ=255)
              cgcCGAGGCTACTAGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGga                                                                                                                                   <  8:239707/149‑1 (MQ=255)
              cgcCGAGGCTACTAGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGga                                                                                                                                   <  4:181920/149‑1 (MQ=255)
              cgcCGAGGCTACTAGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGga                                                                                                                                   <  2:26683/149‑1 (MQ=255)
                 cGAGGCTACTAGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGAc                                                                                                                                <  1:280409/149‑1 (MQ=255)
                      ctactaGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACAtctc                                                                                                                           <  6:253102/149‑1 (MQ=255)
                          taGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACa                                                                                                                       >  3:207487/1‑149 (MQ=255)
                            gCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACAt                                                                                                                      <  6:138569/148‑1 (MQ=255)
                               aGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTgcgg                                                                                                                                               >  4:90763/1‑120 (MQ=255)
                               aGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTgcgg                                                                                                                                               <  3:90763/120‑1 (MQ=255)
                                 gcaaGCAAAGCGCACCCAGAGTAGCTAACGCTCTTGAAACTTCCATATAATCCCTTACGACCTGCAGGAGAACATTCCACGATATAAGATGCTGAGGTGCCCCATTCGCCAGCTGCTGCGAATCCTTGGAGACATCTCAACATGACTAg                                                                                                                <  2:203778/146‑1 (MQ=255)
                                 gCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGaa                                                                                                                                           >  3:29458/1‑122 (MQ=255)
                                 gCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGaa                                                                                                                                           <  4:29458/122‑1 (MQ=255)
                                           cgcCCCCAGAGGACCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACTATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGa                                                                                                      <  1:215742/149‑1 (MQ=255)
                                            gcCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAAc                                                                                                                        >  5:270264/1‑130 (MQ=255)
                                            gcCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAAc                                                                                                                        <  6:270264/130‑1 (MQ=255)
                                              ccccAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGAAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCa                                                                                                   >  3:172039/1‑149 (MQ=255)
                                               cccAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGAc                                                                                                                                >  2:70665/1‑119 (MQ=255)
                                               cccAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCAc                                                                                                  <  7:30583/149‑1 (MQ=255)
                                                ccAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCGCAGGTGGCGaa                                                                                                                                           >  2:254555/1‑107 (MQ=255)
                                                ccAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCGCAGGTGGCGaa                                                                                                                                           >  6:66544/1‑107 (MQ=255)
                                                ccAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCGCAGGTGGCGaa                                                                                                                                           <  5:66544/107‑1 (MQ=255)
                                                ccAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCGCAGGTGGCGaa                                                                                                                                           <  1:254555/107‑1 (MQ=255)
                                                   gaggagCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAAc                                                                                              <  3:308098/149‑1 (MQ=255)
                                                         cTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAAt                                                                                        <  3:254101/149‑1 (MQ=255)
                                                         cTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAAt                                                                                        <  7:176140/149‑1 (MQ=255)
                                                            aCGCTCTGGAAACTTCCAAAGAATTCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGt                                                                                     >  1:231203/1‑149 (MQ=255)
                                                                   ggAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTAt                                                                              >  5:145010/1‑149 (MQ=255)
                                                                      aaCTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGTAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGtt                                                                           <  4:172039/149‑1 (MQ=255)
                                                                      aaCTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGtt                                                                           >  6:41430/1‑149 (MQ=255)
                                                                       aCTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACTATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTg                                                                          <  1:269969/149‑1 (MQ=255)
                                                                            cAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCACCATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAAt                                                                      >  7:297282/1‑148 (MQ=255)
                                                                            cAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAAt                                                                      >  5:321431/1‑148 (MQ=255)
                                                                                gAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGac                                                                 <  4:207487/149‑1 (MQ=255)
                                                                                    cccTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGaccaa                                                              >  6:85889/1‑148 (MQ=255)
                                                                                     ccTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAGTTGGAGCACCAACTCCAATAGTCTCGTATCTTTGAATCAGaccaata                                                            >  7:80781/1‑149 (MQ=255)
                                                                                     ccTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGaccaata                                                            <  5:41430/149‑1 (MQ=255)
                                                                                      cTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACt                                                                                                                  >  4:94328/1‑94 (MQ=255)
                                                                                      cTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACt                                                                                                                  <  3:94328/94‑1 (MQ=255)
                                                                                      cTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGaccaatac                                                           <  6:145010/149‑1 (MQ=255)
                                                                                         aCGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGaccaatacca                                                         <  3:119441/148‑1 (MQ=255)
                                                                                          cGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGaccaataccaat                                                       <  3:152397/149‑1 (MQ=255)
                                                                                                               aCGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAAc                                                                                              <  5:117051/89‑1 (MQ=255)
                                                                                                               aCGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAAc                                                                                              >  6:117051/1‑89 (MQ=255)
                                                                                                               aCGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATc                                                                     >  8:237622/1‑114 (MQ=255)
                                                                                                               aCGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATc                                                                     <  7:237622/114‑1 (MQ=255)
                                                                                                               aCGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCaaa                                   >  4:55705/1‑148 (MQ=255)
                                                                                                                cGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAaata                                 <  3:55705/149‑1 (MQ=255)
                                                                                                                   tataAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAaataat                               <  2:231203/148‑1 (MQ=255)
                                                                                                                        gATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTAGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTcag                         >  2:102538/1‑149 (MQ=255)
                                                                                                                                 tccaaggaTTCGCCGCAGGTGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACTCCAATAGTCTCGTATGTTGGAATAAGACCAATACCAATACTGCCGATTGCCATGAAAATAATCGTCAGCAGCATcgc                >  4:254183/8‑149 (MQ=255)
                                                                                                                                    ccccATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGCggg             >  3:156321/1‑149 (MQ=255)
                                                                                                                                      ggaTTCGCCGCAGGTGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACTAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGCggggt           >  2:316157/3‑147 (MQ=255)
                                                                                                                                         ttCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGaccaataccaat                                                       >  4:4964/1‑102 (MQ=255)
                                                                                                                                         ttCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGaccaataccaat                                                       <  3:4964/102‑1 (MQ=255)
                                                                                                                                             ccACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGCGGGCTTGCGccc    >  2:253993/1‑149 (MQ=255)
                                                                                                                                               aCCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGCGGGCTTGCGccccc  >  8:160087/1‑149 (MQ=255)
                                                                                                                                               aCCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACAAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGCGGGCTTGCGccccc  >  5:298839/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGAATGAGCAAGAGCGCCGAGGCTACTAGCGAGGCCAGCAAAGCGCCCCCAGAGGAGCTCACGCTCTGGAAACTTCCAAAGAATCCCTTACGACCTGCAGGAGACCATTCCACGATATAAGATGCAGAGGTGCCCCATTCGCCACCTGCGGCGAATCCTTGGAGACATCTCAACATGACTAGCAGAATTGGAGCACCAACGCCAATAGTCTCGTATGTTGGAATCAGACCAATACCAATACTGCCGATTGCCATGCAAATAATCGTCAGCAGCATCGCGGGCTTGCGCCCCC  >  NC_002947/5064467‑5064758

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