Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,440,087 Δ2 bp coding (942‑943/990 nt) rpoA ← RNA polymerase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,440,0860C.100.0% 242.4 / NA 52coding (944/990 nt)rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (28/24);  total (28/24)
*NC_000913_3_pae_pgi3,440,0870C.100.0% 242.4 / NA 52coding (943/990 nt)rpoARNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (28/24);  total (28/24)

ACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGT  >  NC_000913_3_pae_pgi/3439939‑3440227
                                                                                                                                                    ||                                                                                                                                             
aCATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGGCAGCGATGCTTCCCGGTGGCCAGTTTTCCAGGCGCATGca                                                                                                                                               <  1:314442/149‑2 (MQ=255)
  aTAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTACTAGCAGGTGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGcgga                                                                                                                                             >  1:475140/1‑146 (MQ=255)
     gCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGcagagac                                                                                                                                          >  1:214492/1‑143 (MQ=255)
      ccTGGGGATGGCTGCTGTTGCGGTTCAGTTGACGACCACGCTGACGATGGCGCATGACCTTATCCCTCTCAGTAAGACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGcagagaca                                                                                                                                         <  1:482071/149‑8 (MQ=255)
        tGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGt                                                                                                                                     >  1:411582/1‑149 (MQ=255)
             aTGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAACTCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACg                                                                                                                                <  1:654850/149‑1 (MQ=255)
              tGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACgg                                                                                                                               <  2:214492/149‑1 (MQ=255)
                gctgctGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTTGCCCGTTTTCCAGGCGCCTG‑‑CAGAGCCATTCCACGGGa                                                                                                                             >  2:171242/1‑149 (MQ=255)
                   gctgTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGc                                                                                                                          >  1:931310/1‑149 (MQ=255)
                      gTTGCGGTTCAGTTGACGACCACTCTTACGACGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTGCTCGTCCGCGATGCTTGCCGGTGGCCAGTGTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAg                                                                                                                       <  2:666183/149‑1 (MQ=255)
                         gCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCGGCAc                                                                                                                    >  1:445026/1‑149 (MQ=255)
                         gCGGTTCAGTGGACGACCACTCGTACGATGCCGCATGACCTCATCCTTCTCAGTAAAACATTAACCTGTGATCCGGTTACTCGTCAGCGTTGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCAc                                                                                                                    <  1:419590/149‑1 (MQ=255)
                           ggTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGt                                                                                                                  >  1:305586/1‑149 (MQ=255)
                             ttCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCt                                                                                                                <  2:445026/149‑1 (MQ=255)
                                  ttGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAt                                                                                                           >  1:241291/1‑149 (MQ=255)
                                    gacgacCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtct                                                                                                         <  2:1070621/149‑1 (MQ=255)
                                       gacCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAg                                                                                                      >  1:890713/1‑149 (MQ=255)
                                            ctctTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAaga                                                                                                 <  1:888251/149‑1 (MQ=255)
                                               ttACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATACGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCCGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAt                                                                                              >  1:484738/1‑149 (MQ=255)
                                                     tGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTa                                                                                        <  2:411582/149‑1 (MQ=255)
                                                            tGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAt                                                                                                           <  2:449591/123‑1 (MQ=255)
                                                            tGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAt                                                                                                           >  1:449591/1‑123 (MQ=255)
                                                                 ttATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGc                                                                            >  1:107020/1‑149 (MQ=255)
                                                                    tCCTTCTCAGTAAATCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGtt                                                                         <  2:308441/149‑1 (MQ=255)
                                                                     ccTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGttt                                                                        <  2:931310/149‑1 (MQ=255)
                                                                        tctcAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGTCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTaa                                                                     <  1:266550/149‑1 (MQ=255)
                                                                        tctcAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTaa                                                                     <  2:241291/149‑1 (MQ=255)
                                                                         ctcAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAg                                                                    <  2:305586/149‑1 (MQ=255)
                                                                         cgcAGGAAAACCTGAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCTTTTTAAg                                                                    <  2:460835/147‑1 (MQ=255)
                                                                                aaaCCTTAACCTGTGATCAGGTTACTCGTCAGCGATGCTTGCCGGTGACCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCCGCACGCCTTTAATCGCAGTAAGAGATTTTTTACCGAGGTTAGGCGTTTTAAGGTGCTCa                                                             >  2:302421/1‑149 (MQ=255)
                                                                                 aaacTTAACCTGTGATCCGGTTACTCGTCACCGATTCTTGCCGGTTGCGAGATTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCaa                                                            >  1:155666/4‑149 (MQ=255)
                                                                                   ccTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAAcc                                                          >  1:1054735/1‑149 (MQ=255)
                                                                                         ccTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTa                                                   >  2:157668/1‑149 (MQ=255)
                                                                                           tgtgATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCGGGCGCGTG‑‑CAGAGACTGTCCACGGGAAGCCAGCACGTCTTTAGTCTCAGTAAGAGGTTTTTTGTCAAGGTTAGGCGTTTTTATGAGCTCAACCTCGGtggg                                                  >  1:1016655/1‑146 (MQ=255)
                                                                                              gATCCGGTTGCTGGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGTAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTg                                              <  2:792436/149‑1 (MQ=255)
                                                                                                tcggtttACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCGGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTa                                            >  2:190540/6‑149 (MQ=255)
                                                                                                tCCGGTTACTCGTCAGCGATTCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTa                                            <  2:475140/149‑1 (MQ=255)
                                                                                                     ttACTCGTCAGCGATGCTTGCCGGTGGTCAGTTTTCCAGGCGCATG‑‑TAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTGAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGa                                       <  1:743024/149‑1 (MQ=255)
                                                                                                         tcggaaGCGATGCTTGCCGGTGGCAAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCATATCAc                                   <  2:405380/144‑1 (MQ=255)
                                                                                                           gTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCg                                 >  1:285123/1‑149 (MQ=255)
                                                                                                            tCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGa                                >  1:867652/1‑149 (MQ=255)
                                                                                                                       ttGCCGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGAt                      >  1:50097/1‑149 (MQ=255)
                                                                                                                           cGGTGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCt                  >  1:452223/1‑149 (MQ=255)
                                                                                                                              tGGCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCt               <  1:464085/149‑1 (MQ=255)
                                                                                                                               ggCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAg                                                                    >  1:964311/1‑96 (MQ=255)
                                                                                                                               ggCCAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAg                                                                    <  2:964311/96‑1 (MQ=255)
                                                                                                                                 ccAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAGTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGGTATAGTGGATAGCTTCTGCt            >  2:799017/1‑149 (MQ=255)
                                                                                                                                  cAGTTTTCCAGGCGCATG‑‑CAGAGACCAACCGCGTTTATCGAGCACGTCTTTAATCTCAGTAAGACATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCtt           >  2:1040181/1‑149 (MQ=255)
                                                                                                                                  cAGTTTTCCAGGCGCATG‑‑CAGAGACAGTCGACGGGAAGCCACGCCGCCCTTAATCTCAGACAGAGAATTTTTACCCAGGTCACGCGTTTTAGGGAGGTCAACCTCGGTACGCTGTACCAGCTCCCCGCTGCAGTGGATCGCTTCTGCtt           >  2:979675/1‑149 (MQ=255)
                                                                                                                                     ttttCCAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTaa        <  1:920473/149‑1 (MQ=255)
                                                                                                                                          cAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCttt          >  2:478187/1‑142 (MQ=255)
                                                                                                                                          cAGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCttt          <  1:478187/142‑1 (MQ=255)
                                                                                                                                          cAGGCGCATG‑‑CAAAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAGTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATAGAGTGGATAGCTTCTGATTTAAGGCAg   >  2:18484/1‑149 (MQ=255)
                                                                                                                                           aGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAATGTTAGGCGTTTTAAGCAGCTCAACCTCGGTACGCTGTACCAGATAACCGATATAGTGGATAGCTTCTGATTTAAGGCAGt  >  2:639958/1‑149 (MQ=255)
                                                                                                                                           aGGCGCATG‑‑CAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGt  <  1:633820/149‑1 (MQ=255)
                                                                                                                                           aGGCGCATG‑‑CAGAGACAGGCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGATCACGACAAGAATATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGt  >  2:1031575/1‑149 (MQ=255)
                                                                                                                                                    ||                                                                                                                                             
ACATAGCCTGGCGATGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGT  >  NC_000913_3_pae_pgi/3439939‑3440227

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