Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,745,446 T→C 6.8% V230A (GTG→GCG)  paoC → aromatic aldehyde dehydrogenase molybdopterin‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,745,4460TC 6.8% 176.0 / 2.6 59V230A (GTG→GCG) paoCaromatic aldehyde dehydrogenase molybdopterin‑binding subunit
Reads supporting (aligned to +/- strand):  ref base T (30/25);  new base C (2/2);  total (32/27)
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.82e-01

CACAACCCGATGGAGCCCCATGCGTCAACCGTACTCTACAAAGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCCAGCAGATGTTCACCTTTGGCTACCGGCCC  >  NC_002947/3745304‑3745594
                                                                                                                                              |                                                                                                                                                    
cacaACCCGATGGAGCCCCATGCGTCAACCGTACTCTACAAAGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGcgcgcg                                                                                                                                               >  2:622979/1‑150 (MQ=255)
cacaACCCGATGGAGCCCCATGCGTCAACCGTACTCTACAAAGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGCGcacgcg                                                                                                                                               >  2:946477/1‑150 (MQ=255)
      ccGATGGAGCCCCATGCGTCAACCGTACTCTACAAAGTCGACGGCACGCTGGAAATCCACGACAAGACCTAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCAGCGTGCGCGCGGCATtc                                                                                                                                         >  1:975418/1‑150 (MQ=255)
        gATGGAGCCCCATGCGTCAACCGTACTCTACAAAGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGCGCACGCggatatttg                                                                                                                                       <  1:946477/150‑8 (MQ=255)
            gAGCCCCATGCGTCAACCGTACTCTACAAAGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTcggc                                                                                                                                   <  1:20268/150‑1 (MQ=255)
                            ccGTACTCTACAAAGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCgg                                                                                                                   <  1:572227/150‑1 (MQ=255)
                             cGTACTCTACAAAGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCggg                                                                                                                  >  1:425073/1‑150 (MQ=255)
                                  tctACAAAGGTGACAGCAATCTGGAAATCCACGACATGACCGAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCAGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTgc                                                                                                             <  2:975418/150‑1 (MQ=255)
                                  tctACAAAGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGACCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTc                                                                                                                          >  2:784433/1‑137 (MQ=255)
                                  tctACAAAGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGACCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTc                                                                                                                          <  1:784433/137‑1 (MQ=255)
                                      cAAAGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGccc                                                                                                         >  1:1009057/1‑150 (MQ=255)
                                         aGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCgg                                                                                                                                              <  2:858688/110‑1 (MQ=255)
                                         aGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCgg                                                                                                                                              >  1:858688/1‑110 (MQ=255)
                                                     gCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGcc                                                                                          >  1:256698/1‑150 (MQ=255)
                                                     gCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGcc                                                                                          >  2:251854/1‑150 (MQ=255)
                                                       tGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGt                                                                                              >  2:939281/1‑144 (MQ=255)
                                                       tGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGt                                                                                              <  1:939281/144‑1 (MQ=255)
                                                              ccACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGt                                                                                 >  2:609655/1‑150 (MQ=255)
                                                              ccACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGt                                                                                 <  2:1009057/150‑1 (MQ=255)
                                                                aCGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGa                                                                               <  1:609655/150‑1 (MQ=255)
                                                                         cccAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGt                                                                      <  2:155453/150‑1 (MQ=255)
                                                                         cccAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGt                                                                      <  2:689950/150‑1 (MQ=255)
                                                                          ccAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTgcg                                                                                                            <  1:793059/111‑1 (MQ=255)
                                                                          ccAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTgcg                                                                                                            >  2:793059/1‑111 (MQ=255)
                                                                          ccAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCa                                                                                                <  2:681163/123‑1 (MQ=255)
                                                                          ccAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCa                                                                                                >  1:681163/1‑123 (MQ=255)
                                                                          ccAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAg                                                                                               >  1:569263/1‑124 (MQ=255)
                                                                          ccAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAg                                                                                               <  2:569263/124‑1 (MQ=255)
                                                                               ggacacaGACCTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCCTTCGTCGCCGGTGCCTTCGGCTCGGGCTTGAGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCAcc                                                                >  2:140065/3‑150 (MQ=255)
                                                                                 acacaGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTg                                                              >  2:520998/1‑150 (MQ=255)
                                                                                   acaGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGTCACTGGCCGTGATGGCGGCGTTGCACCTGCa                                                            >  1:607125/1‑150 (MQ=255)
                                                                                   acaGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCa                                                            <  2:700360/150‑1 (MQ=255)
                                                                                              aGGACTACCTGCACAAGGTGTTTGGCCTGCCCAATGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTTGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTgc                                                 <  2:144863/150‑1 (MQ=255)
                                                                                              aGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTgc                                                 >  2:164078/1‑150 (MQ=255)
                                                                                                     ccTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCg                                                     >  2:505077/1‑139 (MQ=255)
                                                                                                     ccTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCg                                                     <  1:505077/139‑1 (MQ=255)
                                                                                                      cTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACg                                         >  1:188582/1‑150 (MQ=255)
                                                                                                      cTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACg                                         >  1:170002/1‑150 (MQ=255)
                                                                                                        gCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCt                                       >  2:761785/1‑150 (MQ=255)
                                                                                                         cacaAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTg                                      >  1:1082821/1‑150 (MQ=255)
                                                                                                          acaAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGa                                     >  2:220532/1‑150 (MQ=255)
                                                                                                             aGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTc                                  <  2:256698/150‑1 (MQ=255)
                                                                                                              ggTGTTTGGCCTGCCCAAGGCAAATATCCGCATGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGcc                                                                                          <  2:962750/93‑1 (MQ=255)
                                                                                                              ggTGTTTGGCCTGCCCAAGGCAAATATCCGCATGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGcc                                                                                          >  1:962750/1‑93 (MQ=255)
                                                                                                                  tttGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCa                                                            >  1:1002452/1‑119 (MQ=255)
                                                                                                                  tttGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCa                                                            <  2:1002452/119‑1 (MQ=255)
                                                                                                                  tttGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCcag                             >  1:371880/1‑150 (MQ=255)
                                                                                                                    tGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGTCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCcagca                           <  2:607125/150‑1 (MQ=255)
                                                                                                                       ccTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCCAGCAGAt                        >  2:1149310/1‑150 (MQ=255)
                                                                                                                           cccAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCCAGCAGATGTTc                    <  1:164078/150‑1 (MQ=255)
                                                                                                                              aaGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCCAGCAGATGTTCAcc                 >  1:753445/1‑150 (MQ=255)
                                                                                                                                 gCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCCAGCAGATGTTCACCttt              <  2:753445/150‑1 (MQ=255)
                                                                                                                                  cAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCCAGCAGATGTTCACCTTTg             <  1:220532/150‑1 (MQ=255)
                                                                                                                                   aaaTATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCCAGCAGATGTTCACCTTTgg            >  1:105211/1‑150 (MQ=255)
                                                                                                                                      atgccGCGCGCACGCGGCATTCGTCGGCGGTGCCATCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGc                                                                   <  2:815279/89‑1 (MQ=255)
                                                                                                                                      atgccGCGCGCACGCGGCATTCGTCGGCGGTGCCATCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGc                                                                   >  1:815279/4‑92 (MQ=255)
                                                                                                                                       atCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCCAGCAGATGTTCACCTTTGGCTAc        <  2:105211/150‑1 (MQ=255)
                                                                                                                                         ccGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCCAGCAGATGTTCACCTTTGTCTACCg      >  2:117772/1‑150 (MQ=255)
                                                                                                                                             gTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCCAGCAGATGTTCACCTTTGGCTACCGGccc  <  1:251854/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CACAACCCGATGGAGCCCCATGCGTCAACCGTACTCTACAAAGGCGACGGCACGCTGGAAATCCACGACAAGACCCAAGGCACACAGAACTGCCAGGACTACCTGCACAAGGTGTTTGGCCTGCCCAAGGCAAATATCCGCGTGCGCGCGGCATTCGTCGGCGGTGCCTTCGGCTCGGGCTTGCGCCCGCAGTACCAGTTGCCACTGGCCGTGATGGCGGCGTTGCACCTGCAGCGCTCGGTGCGCGTCACGCTGACTCGCCAGCAGATGTTCACCTTTGGCTACCGGCCC  >  NC_002947/3745304‑3745594

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