Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 69,672 +T coding (377/1701 nt) araB ← ribulokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091369,6721.T100.0% 108.4 / NA 33E126? (GAA→GNA) araBribulokinase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (14/19);  total (14/19)

ATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCAC  >  NC_000913/69431‑69921
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         
gcgCCCGGCGCGGCAACGCCCGCAACGAATACCCGGCGGGGGGGGGTTACGGGAAAGCAGAGCGGCCCCCCCGTCGCCCACCGCAATCCACGATGCGGCAGATGCCGCCAGGGCGGTGTCCTGGCGGGTCACAGGCGGTATTTTTGCCCAGACCCATTCGCGGGAATAAATCCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAg                                                                                                                                                                                                                                                  <  2:257950/249‑1 (MQ=255)
         cgtgaaaCGGCCGCGACGAATATCCTGCGGGCGGGTGGCACCGGAAAGCAGAGCTGCCACCCAGTCGAACAGCTCAACCCCCGAGGCGGCAGATTGCGCCACGGCGCTGTCCTGGGGAGTCACATGCAGGATTTTTGCCCAGACCCATCCGCGGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGTTCtt                                                                                                                                                                                                                                         <  2:65206/246‑1 (MQ=255)
             cAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTcc                                                                                                                                                                                                                                      >  1:211181/1‑250 (MQ=255)
               aCGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCcac                                                                                                                                                                                                                                    >  2:212759/1‑250 (MQ=255)
                 gTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAAt                                                                                                                                                                                                                                 >  2:65023/1‑251 (MQ=255)
                      ggacgaATATCCTGCGGGCGGTGGGTACCGGAAAGCAGAGCGGGCACCAATTCGCACAGACAAATCAACGTGGCGGAAGATTGCCCCACGGCGCTGTCCTGGCGGGTAACATGAAGGATTTTTGCCCAGAACCATTCGCTGGAATAAACACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGaa                                                                                                                                                                                                                            <  1:233659/250‑1 (MQ=255)
                      cgacgaATAGCCGGGGGGGGGGTGGGACCGAAAAGCAGAGCTGGCCCCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGaa                                                                                                                                                                                                                            <  2:95784/251‑1 (MQ=255)
                                ccTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTc                                                                                                                                                                                                                  >  1:84797/1‑251 (MQ=255)
                                                  cGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCCCGAGGCGGCAGATTGCGCCACGGCGCTGCCCTGGGGATTCACAGGCAGGATTTTTGCCCAGACCCATGCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTGCGGCAAACTcc                                                                                                                                                                                                <  1:202000/251‑1 (MQ=255)
                                                       aGCAGAGCTGGCCCCCAGTCGCACAGCTCAATCCACGATGCGGCAGATGCCGCCCCGGCGCTTTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGACCCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGgcg                                                                                                                                                                                           <  1:153864/251‑1 (MQ=255)
                                                                                         ccgaggCGGCAGATTGCCCCACGGCGCTGTCCGGGCGATTCACAGGCAGGATTTTTGCCCAAACCCATTCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGt                                                                                                                                                         <  1:10732/246‑1 (MQ=255)
                                                                                             tGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCgg                                                                                                                                                     >  2:142599/1‑251 (MQ=255)
                                                                                                               ggCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGt                                                                                                                                   >  1:38755/1‑251 (MQ=255)
                                                                                                                    tGTCCTGGCGAGCCAAATGATGGTTTTTTGCCAAGAACCTTACCCTGAAAAAAACTCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAAc                                                                                                                               <  2:18064/250‑1 (MQ=255)
                                                                                                                            cGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATccc                                                                                                                      >  1:92513/1‑251 (MQ=255)
                                                                                                                                    aTGCAGGATTTTTGCCCAGAACCATTCGCTGGAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAg                                                                                                              <  1:65023/251‑1 (MQ=255)
                                                                                                                                              tttGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGAGGGAGTAGTCAACATTGCAAGGCGCGTGGCACAACCGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCCCAATACGAACATCGCGTTCGGGCTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACACCAGCTGCGCGCTg                                                                                                    >  1:125842/1‑251 (MQ=255)
                                                                                                                                                 gccaagaacatttCGCTGAAAGAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCTAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGtt                                                                                                  <  1:179603/239‑1 (MQ=255)
                                                                                                                                                            ttCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTAAGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGct                                                                                      >  1:136312/1‑251 (MQ=255)
                                                                                                                                                                 gggAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGGCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCa                                                                                 <  1:10892/250‑1 (MQ=255)
                                                                                                                                                                  gaaataaataCCGCCAATTTAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCGCGTTACCGTCGGCATCAATCGGTGCGGGCGTCGGGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTTTTCGACGCTAAGCTCTGCa                                                                                 <  1:70613/248‑1 (MQ=255)
                                                                                                                                                                  gaaaaaaatcccgccaATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCa                                                                                 <  2:84797/237‑1 (MQ=255)
                                                                                                                                                                    aataaataACGCCAATATAGCGGGAGCGGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGACTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGTGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGc                                                                              >  2:102407/1‑251 (MQ=255)
                                                                                                                                                                                   tatAGCGGGAGCAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGATTCTTTCAACCGCAGGGTGGTCTTTCAACAATACGAACATCGCGTTCGGGTATTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGAGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCCGCGCGCTGTTCGACGCTAAGATCTGCAAGCACGGTTTTCAGTGcg                                                               >  2:237768/3‑250 (MQ=255)
                                                                                                                                                                                                   aacattGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCa                                                >  2:234338/5‑250 (MQ=255)
                                                                                                                                                                                                    acattGCCCGGCGCGTGGCAAAAACGGGTAATCTCTTCGCTTCTTTCAACCGACGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGTTTTTCGGCAAACTCAGGGCGCAGCGCCAACACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAAt                                              >  2:214938/4‑251 (MQ=255)
                                                                                                                                                                                                           ccGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCTTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAATCACGGTTTTCAGTGCCGCTTCCATTGCCTCAATGTAGTCa                                       <  2:20760/251‑1 (MQ=255)
                                                                                                                                                                                                           caggcgcgcagCACAAACGGGTAATATCTTAGCTTCTTTCACCGGCAGTGTGGTCTTTCCACAATACGAACATCGCGATCGGCTTTCCGGCAAACTCAGGGCGCAGCTCCAGCACGTTACCGTCGGCATCAATCGGTGCGGCCGTCGAGCCGGTAATGTCAACGCCAATCCCGACCACAGCTGCGAGCTGTTCGACGCTAACCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCATTGTAGTc                                        >  2:6210/11‑250 (MQ=255)
                                                                                                                                                                                                                    ggCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGgatg                               <  1:191692/250‑1 (MQ=255)
                                                                                                                                                                                                                      cacaAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGatgatg                            <  2:97976/251‑1 (MQ=255)
                                                                                                                                                                                                                        caAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCCCGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGAc                          <  1:138188/251‑1 (MQ=255)
                                                                                                                                                                                                                                         tcGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTAtt             <  1:102096/247‑1 (MQ=255)
                                                                                                                                                                                                                                           gCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCgggg        <  1:183445/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cttctttcAACCGCAGTGTGGTCTTTCCACAATACGACCATCGCGTTCGGGTTCTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGCGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGc       >  1:40099/6‑251 (MQ=255)
                                                                                                                                                                                                                                                 tcAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATcac  <  1:121948/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         
ATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCAC  >  NC_000913/69431‑69921

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