Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,888,564 T→C T512T (ACA→ACG cysI ← sulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,888,5640TC100.0% 192.8 / NA 58T512T (ACA→ACGcysIsulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/31);  total (27/31)

CGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTT  >  NC_000913/2888420‑2888709
                                                                                                                                                |                                                                                                                                                 
cGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCa                                                                                                                                             <  1:538787/151‑1 (MQ=255)
  ggCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAAt                                                                                                                                           >  1:937383/1‑151 (MQ=255)
   gCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATg                                                                                                                                          >  2:600177/1‑151 (MQ=255)
    cGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATg                                                                                                                                          <  1:479969/150‑1 (MQ=255)
                   cgcacCGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGAt                                                                                                                          <  2:986927/151‑1 (MQ=255)
                     cacCGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGCTACGCGTCCCAATGCGGTTGCCGACAAGATg                                                                                                                         >  1:836/1‑150 (MQ=255)
                       ccGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAg                                                                                                                      >  2:642568/1‑151 (MQ=255)
                       ccGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAg                                                                                                                      >  1:484904/1‑151 (MQ=255)
                          tAAAGTCGCCGAAGCCTTCACCCGCTTCNNGCNNNTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGt                                                                                                                    <  1:1017588/150‑1 (MQ=255)
                          tAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGt                                                                                                                    <  1:692322/150‑1 (MQ=255)
                          tAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGt                                                                                                                    <  2:1019159/150‑1 (MQ=255)
                                        ccTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGgcg                                                                                                     >  2:811895/1‑151 (MQ=255)
                                        ccTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGgcg                                                                                                     >  2:726852/1‑151 (MQ=255)
                                            cACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCttt                                                                                                 >  1:174886/1‑151 (MQ=255)
                                                   ttCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACccacc                                                                                          <  1:274852/151‑1 (MQ=255)
                                                   ttCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACccacc                                                                                          <  1:68840/151‑1 (MQ=255)
                                                   ttCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACccacc                                                                                          <  2:703875/151‑1 (MQ=255)
                                                   tgcgcgCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACAACCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACccacc                                                                                          <  1:771411/149‑1 (MQ=255)
                                                           ctTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCa                                                                                   >  2:234523/1‑150 (MQ=255)
                                                           ctTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCa                                                                                   >  2:137450/1‑150 (MQ=255)
                                                             ttCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACtt                                                                                <  2:633273/151‑1 (MQ=255)
                                                              tCGCCCAGCGCCCTATCAGTTCATCAAGCGTCGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTc                                                                               <  1:972831/151‑1 (MQ=255)
                                                              tCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTc                                                                               <  1:89533/151‑1 (MQ=255)
                                                                gCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCCGGCCCACTTCCg                                                                             >  2:997367/1‑151 (MQ=255)
                                                                     gcgcCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAg                                                                         >  2:472189/1‑150 (MQ=255)
                                                                     gcgcCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGGGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCCCTTTACCCACCAGGCCCACTTCCGCCAg                                                                         >  2:1054692/1‑150 (MQ=255)
                                                                       ccccTATCAGTTCATCAAGCGACGCCAGTATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCTCTTTACCCACCAGGCCCACTTCCGCCAGCAt                                                                      <  2:485105/150‑1 (MQ=255)
                                                                        cccTATCAGTTCATCAAGCGGCGCCAGGATTTCCGGCTCGGTGATGGTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCAt                                                                      <  2:322310/150‑1 (MQ=255)
                                                                        cccTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCAt                                                                      <  1:984021/150‑1 (MQ=255)
                                                                              cAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGa                                                               >  1:302887/1‑151 (MQ=255)
                                                                                  tcatcaAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACca                                                            >  2:33775/1‑150 (MQ=255)
                                                                                   catcaAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACcaca                                                          >  1:1031080/1‑151 (MQ=255)
                                                                                   catcaAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCAGCGCGCGACcaca                                                          >  2:1030404/1‑151 (MQ=255)
                                                                                    atcaAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACaa                                                         >  2:230756/1‑151 (MQ=255)
                                                                                     tcaAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAAc                                                        <  1:86188/151‑1 (MQ=255)
                                                                                     tcaAGCGACGCCAGGATTTCCGGCGCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCCCCAGGCCCACTTCCGCCAGCATCGCGCGACCACAAc                                                        <  2:353278/151‑1 (MQ=255)
                                                                                            aCGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCgg                                                 >  2:95103/1‑151 (MQ=255)
                                                                                               ccAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGc                                               <  1:175993/150‑1 (MQ=255)
                                                                                                    aTTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGc                                            >  1:714755/1‑148 (MQ=255)
                                                                                                        ccGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTac                                     <  1:635742/151‑1 (MQ=255)
                                                                                                        ccGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTac                                     <  2:291450/151‑1 (MQ=255)
                                                                                                           gCTCGGTGATGTTTTCTTTATACATCGGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTacac                                   <  1:325421/150‑1 (MQ=255)
                                                                                                           gCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCCGCCTGTTACACg                                  >  2:373901/1‑151 (MQ=255)
                                                                                                                gTGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCa                             >  2:392062/1‑151 (MQ=255)
                                                                                                                 tGATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCAc                            <  2:702344/151‑1 (MQ=255)
                                                                                                                  gATGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACg                           >  1:368330/1‑151 (MQ=255)
                                                                                                                   aTGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACg                           <  1:814428/150‑1 (MQ=255)
                                                                                                                    tGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGat                         <  2:612306/151‑1 (MQ=255)
                                                                                                                    tGTTTTCTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGat                         <  2:736418/151‑1 (MQ=255)
                                                                                                                          cTTTATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTc                   <  1:498255/151‑1 (MQ=255)
                                                                                                                            ttATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTc                                                                               <  2:738181/89‑1 (MQ=255)
                                                                                                                            ttATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTc                                                                               >  1:737920/1‑89 (MQ=255)
                                                                                                                            ttATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCAt                 >  1:324636/1‑151 (MQ=255)
                                                                                                                            ttATACATCCGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCAt                 >  1:279722/1‑151 (MQ=255)
                                                                                                                             tataCTACCGTGGGAGACGCGTCCCATTGCGGTTGCCGCCAAGCTGCAGGTTGTAGCGACCCGTCGCTTTACCCACCAGCCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATc                <  1:45909/151‑1 (MQ=255)
                                                                                                                                   tCCGTTGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGCTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCCCATCGCTGa           >  1:421887/1‑150 (MQ=255)
                                                                                                                                    ccGTGGGATACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCCCCAGGCCCCCTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGaca         <  2:261271/151‑1 (MQ=255)
                                                                                                                                           aTACGCGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGtt  <  1:642319/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTTTCGCCCAGCGCCCTATCAGTTCATCAAGCGACGCCAGGATTTCCGGCTCGGTGATGTTTTCTTTATACATCCGTGGGATACGTGTCCCAATGCGGTTGCCGCCAAGATGCAGGTTGTAGCGACCCGGCGCTTTACCCACCAGGCCCACTTCCGCCAGCATCGCGCGACCACAACCGTTCGGGCAGCCTGTTACACGCATCACGATATGCTCATCGCTGACACCATGTT  >  NC_000913/2888420‑2888709

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