Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,505,488 A→G V178A (GTG→GCG)  yhfS ← FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,505,4880AG100.0% 216.1 / NA 58V178A (GTG→GCG) yhfSFNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/30);  total (28/30)

TCTGGCTACCGCCGGAGTAAAGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCACCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGCGCCGCATCCGGTTGTTGCTCGTCGACG  >  NC_000913/3505342‑3505634
                                                                                                                                                  |                                                                                                                                                  
tCTGGCTACCGCCGGAGTAAAGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCtt                                                                                                                                                >  2:799540/1‑151 (MQ=255)
 cTGGCTACCGCCGGAGTAAAGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTtc                                                                                                                                               >  1:534969/1‑151 (MQ=255)
   ggCTACCGCCGGAGTAAAGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCtt                                                                                                                                                <  2:801249/148‑1 (MQ=255)
       aCCGCCGGAGTAAAGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTtcatca                                                                                                                                           <  1:139307/149‑1 (MQ=255)
       aCCGCCGGAGTAAAGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCAcc                                                                                                                                         >  1:305218/1‑151 (MQ=255)
        ccgccgGAGTAAAGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCg                                                                                                                                        <  1:241891/151‑1 (MQ=255)
          gccgGAGTAAAGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCa                                                                                                                                      <  1:973625/151‑1 (MQ=255)
           ccgGAGTAAAGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCAt                                                                                                                                     <  2:656192/151‑1 (MQ=255)
            cgGAGTAAAGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATa                                                                                                                                    <  1:44162/151‑1 (MQ=255)
                 tAAAGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGt                                                                                                                               >  1:765568/1‑151 (MQ=255)
                  aaaGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTc                                                                                                                              >  2:362514/1‑151 (MQ=255)
                          gcgcGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCgg                                                                                                                         >  2:1042143/1‑148 (MQ=255)
                           cgcgAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTa                                                                                                                      >  1:308722/1‑150 (MQ=255)
                           cgcgAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTa                                                                                                                      >  2:1039081/1‑150 (MQ=255)
                             cgAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACg                                                                                                                   >  1:521863/1‑151 (MQ=255)
                              gAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGc                                                                                                                  <  1:396953/151‑1 (MQ=255)
                                  aCGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTgg                                                                                                               <  1:92621/150‑1 (MQ=255)
                                  aCGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTgg                                                                                                               >  1:861444/1‑150 (MQ=255)
                                        gATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCt                                                                                                        <  1:276969/151‑1 (MQ=255)
                                         aTAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTg                                                                                                       <  2:854196/151‑1 (MQ=255)
                                           aaCATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGcc                                                                                                     >  1:880260/1‑151 (MQ=255)
                                                cAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCcgcgc                                                                                                 <  1:111389/150‑1 (MQ=255)
                                                   cATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCa                                                                                              >  2:915746/1‑150 (MQ=255)
                                                   cATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCa                                                                                              >  1:948668/1‑150 (MQ=255)
                                                       gCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGtt                                                                                         <  1:970843/151‑1 (MQ=255)
                                                        ccGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGACGCGCGCAACGTTg                                                                                        >  1:356746/1‑151 (MQ=255)
                                                         cGACCACTGCACCAACACCCTCTGACCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGc                                                                                       <  2:115749/151‑1 (MQ=255)
                                                           aCCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCa                                                                                     <  2:629316/151‑1 (MQ=255)
                                                            ccACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAg                                                                                    <  2:848094/151‑1 (MQ=255)
                                                               cTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGcac                                                                                 >  2:8836/1‑151 (MQ=255)
                                                                  cACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACAt                                                                               <  1:899897/150‑1 (MQ=255)
                                                                  cACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATc                                                                              >  2:794641/1‑151 (MQ=255)
                                                                       acacCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCa                                                                         <  1:776145/151‑1 (MQ=255)
                                                                                ggCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTa                                                                  <  1:862344/149‑1 (MQ=255)
                                                                                         aGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCg                                                       <  1:775661/151‑1 (MQ=255)
                                                                                         aGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCg                                                       >  1:616230/1‑151 (MQ=255)
                                                                                         aGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCg                                                       <  1:672699/151‑1 (MQ=255)
                                                                                              gAAGCAGGAAAATGTCGAGACATGCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGgctg                                                   <  2:639534/150‑1 (MQ=255)
                                                                                                     gAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGcgcg                                           <  2:703210/151‑1 (MQ=255)
                                                                                                      aaaaTGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGt                                          >  2:700625/1‑151 (MQ=255)
                                                                                                        aaTGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTAt                                        <  1:717993/151‑1 (MQ=255)
                                                                                                          tGTCGAGACATTCGCGCCGCATTCACAGCCGCTTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTAtgc                                      >  1:635635/1‑151 (MQ=255)
                                                                                                           gTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTAtgct                                     <  1:709791/151‑1 (MQ=255)
                                                                                                             cGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTAtgctg                                    <  1:597785/150‑1 (MQ=255)
                                                                                                             cGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTAtgctg                                    <  2:957378/150‑1 (MQ=255)
                                                                                                              gagaCATTCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTAtgctgc                                   >  1:268975/1‑150 (MQ=255)
                                                                                                                    ttCGCGCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAgcg                            >  1:991593/1‑151 (MQ=255)
                                                                                                                      cgcgCCGCATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGcgcc                          >  1:717504/1‑151 (MQ=255)
                                                                                                                        cgccgcATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGcgccgc                        <  2:280072/151‑1 (MQ=255)
                                                                                                                        cgccgcATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGcgccgc                        <  1:419254/151‑1 (MQ=255)
                                                                                                                        cgccgcATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGcgccgc                        >  1:176335/1‑151 (MQ=255)
                                                                                                                            gcATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGCGCCGCATcc                    >  2:494475/1‑151 (MQ=255)
                                                                                                                            gcATTCACAGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGCGCCGCATcc                    >  1:932210/1‑151 (MQ=255)
                                                                                                                                 cacaGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGCGCCGCATCCGgtt                <  2:409058/150‑1 (MQ=255)
                                                                                                                                 cacaGCCGATTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGACAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGCGCCGCATACGgt                 >  2:185598/1‑149 (MQ=255)
                                                                                                                                         aTTCGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGCGCCGCATCCGGTTGTTGCtcgt       >  2:90134/1‑151 (MQ=255)
                                                                                                                                            cGAGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGCGCCGCATCCGGTTGTTGCtcgtcg     <  1:192222/150‑1 (MQ=255)
                                                                                                                                              aGCCGCCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGCGCCGCATCCGGTTGTTGCTCGTcgacg  >  1:214610/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCTGGCTACCGCCGGAGTAAAGCGTGGCGCGAATACGGTTGATAACATCAGCATCGCCGACCACTGCACCAACACCCTCTGGCCCAAATAGCTTGAAGCAGGAAAATGTCGAGACATTCGCGCCGCATTCACAGCCGATTCGAGCCACCTTCATCACCGCATAGTTGTCATCGGTTAACGCTGGAACACCTGCCGCGCGCAACGTTGCCAGCACATCTGCCAGCACGTAGCTGTCCTGCGGCTGCTGGCGCGTATGCTGCACCAGCGCCGCATCCGGTTGTTGCTCGTCGACG  >  NC_000913/3505342‑3505634

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