Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 1,494,281 G→T 7.3% A264D (GCC→GAC)  ydcI ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA1,494,2810GT 7.3% 141.8 / 2.3 56A264D (GCC→GAC) ydcIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base T (2/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

GATCCCTTCACCTTTCAGAACGGCATTGATTTTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGA  >  NC_000913_3_pae_tpiA/1494134‑1494429
                                                                                                                                                   |                                                                                                                                                    
gATCACTACACGTTTAAGAACAGCATTAATTATCGAATAGCGTTAATCATCAACTGTGAAACAGAGAAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGTTTCGCCTGCACAATGTCCCGGAACAAACAGGGACACCAGCGTGGc                                                                                                                                                     >  2:689764/1‑149 (MQ=255)
    ccTTCACCTTTCAGAACGGCATTGATTTTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCCTGa                                                                                                                                                 >  1:815591/1‑149 (MQ=255)
      ttCACCTTTCAGAACGGCATTGATTTTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACg                                                                                                                                               >  1:240706/1‑149 (MQ=255)
      ttCACCTTTCAGAACGGCATTGATTTTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACg                                                                                                                                               >  1:545265/1‑149 (MQ=255)
      ttCACCTTTCAGAACGGCATTGATTTTCGAATAGCGTTAATCATCAACTGCCAACGATGATAGAACGTCGCATCTACGCGGGTCGGTATTCCAATCGGTTAGCCTGCACCATGTCCCGGACCAGGCAGGGCCACCACCGTGGCacgacg                                                                                                                                               >  1:185758/1‑149 (MQ=255)
         aCCTTTCAGAACGGCATTGATTTTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAg                                                                                                                                            <  1:708380/149‑1 (MQ=255)
            tttCAGAACGGCATTGATTTTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtc                                                                                                                                         >  2:921337/1‑149 (MQ=255)
               cAGAACGGCATTGATTTTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtcgtc                                                                                                                                      <  1:753919/149‑1 (MQ=255)
                aGAACGGCATTGATTTTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTACCGGAACAGGAAGCGCCACCAGCGTGGCCTGAGGCAGgtcagct                                                                                                                                     >  2:731867/1‑145 (MQ=255)
                    cGGCATTGATTTTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTa                                                                                                                                 <  2:719582/149‑1 (MQ=255)
                         ttGATTTTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCACATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGTCAGGGCCACCAGCGTGGAATGACGCAGGTCGGCTGTTAAAgag                                                                                                                            >  2:94650/1‑147 (MQ=255)
                           gATTTTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCGTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGGCTTTTACAGCGcc                                                                                                                          >  1:513430/1‑149 (MQ=255)
                               ttCGAATAGCGTTAATCATCAACTGGTAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGaa                                                                                                                      <  1:389572/149‑1 (MQ=255)
                               ttCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGaa                                                                                                                      <  2:246122/149‑1 (MQ=255)
                               gtCGAATAGCTTTAATCACCCACTGGCACCCAGAAGAGAACGTCGCATCTTCGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACGGGGGTGGCATGACGTAGGTCGTCTTTTACAGCGCCAGaa                                                                                                                      <  2:8002/148‑1 (MQ=255)
                                tCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAg                                                                                                                     <  1:825799/149‑1 (MQ=255)
                                 cGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAgg                                                                                                                    <  1:673343/149‑1 (MQ=255)
                                    aTAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGAc                                                                                                                 <  1:921337/149‑1 (MQ=255)
                                     tAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGACCAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACa                                                                                                                >  1:67634/1‑149 (MQ=255)
                                            aaTCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACcaca                                                                                                         <  1:762802/149‑1 (MQ=255)
                                               catcaACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTACAACCAGTCCCCCTGCACCATGACACGGACCAGGCAGGGACACCAGCGTGGCATGACGTAGGTCGTGGTGTACAGCGCCAGAAGGGACAAACCACACga                                                                                                      >  2:82637/1‑148 (MQ=255)
                                                atcaACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGGCAGTATTCCAATCGGTTCGCCTGCCCCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCCGGTCGTCTTTTACAGCGCCAGAAGAGAAAAACCAAACGTa                                                                                                     >  2:356256/1‑149 (MQ=255)
                                                 tcaACTGGCAACCAGAAGAGCACGAGGCATCCATGTAGGGCAGTATTCCAACCGGTTCGGTTGCACCAGGTGGCGGCAGAGGCAGTACAACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGtca                                                                                                    >  2:9034/1‑147 (MQ=255)
                                                          aaCCAGAAGAGAAAGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCa                                                                                           >  1:139176/1‑149 (MQ=255)
                                                           aCCAGAAGAGAACGTCGCATCAACGCGGGTAAGTATAGAAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCaa                                                                                          >  2:631705/1‑149 (MQ=255)
                                                               gaagaGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGt                                                                                      >  2:680439/1‑149 (MQ=255)
                                                                            cATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGAt                                                                         >  1:393884/1‑149 (MQ=255)
                                                                             aTCTACCCGGGTCAGCATTCAGGCCCGTGCGCCGGCACCATGTCCCGGCACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATa                                                                        <  1:258917/149‑1 (MQ=255)
                                                                                tACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCACGTCGTCTTTTACAGCGCCAGAAGGGACAAACGACACGTAATCGTATTCGACCGTAAATTGCGGAGATAACg                                                                     >  2:33460/1‑149 (MQ=255)
                                                                                  cgcgGGTGAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGaa                                                                   <  2:472443/149‑1 (MQ=255)
                                                                                    cgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCAACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGc                                                                 <  2:702562/149‑1 (MQ=255)
                                                                                     gggTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCCGGCTGACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCa                                                                <  1:727238/149‑1 (MQ=255)
                                                                                     gggACAGTATTCCAATCCGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCCTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCa                                                                <  2:183756/149‑1 (MQ=255)
                                                                                        tCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGAc                                                             <  2:760682/149‑1 (MQ=255)
                                                                                         cAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACa                                                            <  2:572337/149‑1 (MQ=255)
                                                                                               tCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGttt                                                      <  1:291448/149‑1 (MQ=255)
                                                                                                  aaTCGGTTCTCCTGAAAACTATTCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGa                                                   >  2:670712/1‑149 (MQ=255)
                                                                                                  aaTCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCTGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCATACGTAATCGTATTCAACCGTAAGTTGACAAGTTAGCGAAGCAGACAGCGTTTCGa                                                   <  1:536917/149‑1 (MQ=255)
                                                                                                    tCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTAGATa                                                 >  2:807007/1‑149 (MQ=255)
                                                                                                       gTTCGCGTGCACCGGGTCCCGGAACAGGCAGGAGCAGTAGCGTGCGATGACGCAGGTAGTCTTATACAGCGCCAGAAGGGACAAACCACACGTAATCGTATACAACGGTAAGATGACGAGATAGCGAAGCAGGAAGCGTGTCGGTACaa                                              >  2:292889/1‑149 (MQ=255)
                                                                                                        ttCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAAc                                             >  1:474376/1‑149 (MQ=255)
                                                                                                            ccTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGTCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGa                                         >  2:639899/1‑149 (MQ=255)
                                                                                                            ccTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGTCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGa                                         >  2:639960/1‑149 (MQ=255)
                                                                                                            ccTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGa                                         >  1:161240/1‑149 (MQ=255)
                                                                                                                   cATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAAtt                                  <  1:826939/149‑1 (MQ=255)
                                                                                                                    aTGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGAAAGCGTTTCGATACAACCCGAAGGAAttt                                 >  2:877981/1‑149 (MQ=255)
                                                                                                                      gTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTTACGAGATAGCGAAGCAGACAGCGTTTCGATAGAACCCGAAGGACTTTTa                               >  2:3852/1‑149 (MQ=255)
                                                                                                                        cccGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACa                             >  2:261545/1‑149 (MQ=255)
                                                                                                                           ggAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACAGcc                          <  1:857492/149‑1 (MQ=255)
                                                                                                                            gAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAAATTTACAGccc                         >  1:154867/1‑149 (MQ=255)
                                                                                                                                     gggCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACAGCCCTGGCTCTGc                <  1:244482/149‑1 (MQ=255)
                                                                                                                                       gCCACCAGCGTGTCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACAGCCCTGGCTCTGCAc              <  1:639899/149‑1 (MQ=255)
                                                                                                                                       gCCACCAGCGTGTCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACAGCCCTGGCTCTGCAc              <  1:639960/149‑1 (MQ=255)
                                                                                                                                       gCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACAGCCCTGGCTCTGCAc              <  1:759406/149‑1 (MQ=255)
                                                                                                                                       gCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACAGCCCTGGCTCTGCAc              >  1:43108/1‑149 (MQ=255)
                                                                                                                                             agcttcgCATGACGCATGTCGTCTTTGCCAGCGCCAGAAGTGACAATGCACACGGAGTCGTATTCAGTCGTTAGTTTTCGATATAGCTAAGCATACAGCGTTTCCATCCAACCCGAAGGAATTTTACAGCCCTGGCTCTGCACTAATGc        <  2:287967/143‑1 (MQ=255)
                                                                                                                                             aGCGTGGCATGCCGCAGGTCGTCTTTCACAGCGGCGGAAGTGCCAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACAGCCCTGGCTCTGCACTAATGc        <  2:67634/149‑1 (MQ=255)
                                                                                                                                                   gCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGa  <  2:42860/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GATCCCTTCACCTTTCAGAACGGCATTGATTTTCGAATAGCGTTAATCATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAGAAGGGACAAACCACACGTAATCGTATTCAACCGTAAGTTGACGAGATAGCGAAGCAGACAGCGTTTCGATACAACCCGAAGGAATTTTACAGCCCTGGCTCTGCACTAATGCATCTGA  >  NC_000913_3_pae_tpiA/1494134‑1494429

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