Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,601,452 C→A 7.1% G84G (GGC→GGA hyfA → hydrogenase 4 component A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,601,4520CA 7.1% 185.7 / 3.3 56G84G (GGC→GGAhyfAhydrogenase 4 component A
Reads supporting (aligned to +/- strand):  ref base C (23/29);  new base A (2/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCT  >  NC_000913/2601304‑2601598
                                                                                                                                                    |                                                                                                                                                  
cgcgCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCt                                                                                                                                                   <  2:1000574/150‑1 (MQ=255)
 gcgcCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTg                                                                                                                                                  <  1:1127480/150‑1 (MQ=255)
 gcgcCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTAAACGAAAGCCTCTGTATTGGCTg                                                                                                                                                  >  2:1307281/1‑150 (MQ=255)
   gcCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAAGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCa                                                                                                                                                <  1:454477/150‑1 (MQ=255)
     ctggccctggcGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAg                                                                                                                                              <  1:1019691/150‑1 (MQ=255)
     ctggccctggcGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAg                                                                                                                                              <  1:779791/150‑1 (MQ=255)
     ctggccctggcGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAg                                                                                                                                              <  1:770482/150‑1 (MQ=255)
        gccctggcGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGc                                                                                                                                                 <  2:942595/144‑1 (MQ=255)
        gccctggcGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGc                                                                                                                                                 >  1:942595/1‑144 (MQ=255)
          cctggcGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTg                                                                                                                                         <  2:990151/150‑1 (MQ=255)
           ctggcGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTgc                                                                                                                                        >  1:1028329/1‑150 (MQ=255)
             ggcGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTgcgc                                                                                                                                      >  2:626463/1‑150 (MQ=255)
                    aCGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTc                                                                                                                               <  1:3041/150‑1 (MQ=255)
                    aCGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTc                                                                                                                               <  2:1116740/150‑1 (MQ=255)
                    aCGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTc                                                                                                                               <  1:282483/150‑1 (MQ=255)
                     cGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCt                                                                                                                              >  2:205622/1‑150 (MQ=255)
                     cGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCt                                                                                                                              <  2:78542/150‑1 (MQ=255)
                          acaaTCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCa                                                                                                                         >  2:1281013/1‑150 (MQ=255)
                              tCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTg                                                                                                                     <  1:557626/150‑1 (MQ=255)
                               cACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTgg                                                                                                                    <  2:822065/150‑1 (MQ=255)
                                  tGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGgcgc                                                                                                                 <  2:387135/150‑1 (MQ=255)
                                  tGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTACAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGgcgc                                                                                                                 <  1:97577/150‑1 (MQ=255)
                                          tCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAgcg                                                                                                         <  2:172418/150‑1 (MQ=255)
                                              gtgtCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTc                                                                                                     <  1:1072957/150‑1 (MQ=255)
                                                  catcaCTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGa                                                                                                 <  1:1281013/150‑1 (MQ=255)
                                                  catcaCTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGa                                                                                                 >  1:988343/1‑150 (MQ=255)
                                                    tcaCTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAg                                                                                               >  1:874532/1‑150 (MQ=255)
                                                           gAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCg                                                                                        <  1:626463/150‑1 (MQ=255)
                                                              gAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGTCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCGGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTg                                                                                     <  1:1043829/150‑1 (MQ=255)
                                                                          cTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATgcg                                                                         >  2:1130135/1‑150 (MQ=255)
                                                                          cTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATgcg                                                                         <  1:1350592/150‑1 (MQ=255)
                                                                           tGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATgcgc                                                                        <  2:988343/150‑1 (MQ=255)
                                                                            acAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCa                                                                       >  1:1157049/2‑150 (MQ=255)
                                                                             cAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGAGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCaa                                                                      >  1:976769/1‑150 (MQ=255)
                                                                               ggTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAata                                                                    >  2:140383/1‑150 (MQ=255)
                                                                                   tGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATg                                                                           <  1:346034/139‑1 (MQ=255)
                                                                                   tGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATg                                                                           >  2:346034/1‑139 (MQ=255)
                                                                                        ggTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCa                                                           >  1:534844/1‑150 (MQ=255)
                                                                                               gCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGaa                                                    >  2:775268/1‑150 (MQ=255)
                                                                                                 cATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCGGCAAGCGTTGCGCCGGGGGCTGCCCATTTGGCGCAATCAGCGCTGCAGGAAGCCGTCCGGGGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAgg                                                  <  2:1220290/150‑1 (MQ=255)
                                                                                                     tctcAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGctca                                              >  2:354906/1‑150 (MQ=255)
                                                                                                     tctcAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGctca                                              >  2:359437/1‑150 (MQ=255)
                                                                                                        cAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGctcactc                                           >  2:334053/1‑150 (MQ=255)
                                                                                                           agGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCaaa                                        >  1:813032/1‑150 (MQ=255)
                                                                                                           agGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGATGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCaaa                                        >  2:4391/1‑150 (MQ=255)
                                                                                                           agGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGATGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCaaa                                        >  2:9670/1‑150 (MQ=255)
                                                                                                                       aTCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCgaagaa                            <  1:775268/150‑1 (MQ=255)
                                                                                                                          cAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAAc                         <  2:1028329/150‑1 (MQ=255)
                                                                                                                          cAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAAc                         >  1:685207/1‑150 (MQ=255)
                                                                                                                             cTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAAcgcg                      <  1:1130135/150‑1 (MQ=255)
                                                                                                                                  cGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCcaac                 >  2:6466/1‑150 (MQ=255)
                                                                                                                                      aGCCTCTGTATTGGATGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCcaacacaa             <  1:9670/150‑1 (MQ=255)
                                                                                                                                      aGCCTCTGTATTGGATGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCcaacacaa             <  1:4391/150‑1 (MQ=255)
                                                                                                                                        ccTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCcaacacaaca           >  2:1389784/1‑150 (MQ=255)
                                                                                                                                                ttGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTtgc   <  2:874532/150‑1 (MQ=255)
                                                                                                                                                 tGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACCACATGCTTtgct  >  1:522859/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGCGCCTGGCCCTGGCGAAGACGTCAACAATCACTGCCCCTGTCGTGTGTCATCACTGTGAGGAAGCCCCTTGCCTGCAGGTCTGCCCGGTCAATGCCATCTCTCAGAGGGATGATGCGATCCAACTCAACGAAAGCCTCTGTATTGGCTGCAAGCTTTGCGCCGTGGTCTGCCCATTTGGCGCAATCAGCGCTTCAGGAAGCCGTCCGGTGAATGCCCATGCGCAATATGTTTTTCAGGCTGAAGGCTCACTCAAAGACGGCGAAGAAAACGCGCCAACACAACATGCTTTGCT  >  NC_000913/2601304‑2601598

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