Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,350,540 C→T 7.4% intergenic (+29/+53) PP_4704 → / ← PP_4705 hypothetical protein/glyoxalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,350,5400CT 7.4% 157.1 / 2.8 54intergenic (+29/+53)PP_4704/PP_4705hypothetical protein/glyoxalase
Reads supporting (aligned to +/- strand):  ref base C (31/19);  new base T (2/2);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCACCA  >  NC_002947/5350394‑5350689
                                                                                                                                                  |                                                                                                                                                       
cTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTg                                                                                                                                                    <  1:621545/150‑1 (MQ=255)
       gTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTa                                                                                                                                             >  2:1323889/1‑150 (MQ=255)
        ggtgacggcgcacgcagcttccacctgctggacccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCTCCCACAGGCCTGCGCACTAc                                                                                                                                            >  2:731697‑M1/124‑150 (MQ=255)
          gACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGc                                                                                                                                                         >  2:587981/1‑137 (MQ=255)
          gACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGc                                                                                                                                                         <  1:587981/137‑1 (MQ=255)
             aaGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCtgtg                                                                                                                                       >  2:861421/1‑150 (MQ=255)
              aGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTgg                                                                                                                                      <  1:125165/150‑1 (MQ=255)
                cgcacgcagcttccacctgctggacccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGa                                                                                                                                    >  2:1249711‑M1/116‑150 (MQ=255)
                  cAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGc                                                                                                                                  <  1:659021/150‑1 (MQ=255)
                   aaGCTTGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCg                                                                                                                                 <  1:340934/150‑1 (MQ=255)
                        ggAAGTCGACTAGGTGTACTTCATTGATACTAGCCTGGGGCAGCCTCAAAGAGATGACCGGGTAGTGTCACCGATGTGCCTGGATACAGAGTCATCCCCTCTGGCGGGACAGCCCCGCGAGGCCTGCGCGCTATCTGTGGGAGGGGGCGt                                                                                                                            >  2:681911/1‑150 (MQ=255)
                        gcttccacctgctggacccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGt                                                                                                                            >  2:782840‑M1/108‑150 (MQ=255)
                                      tgtACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGc                                                                                                              <  1:1323889/150‑1 (MQ=255)
                                        tACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCa                                                                                                            >  1:762406/1‑150 (MQ=255)
                                         ccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCtgtg                                                                                                                                       >  1:363118‑M1/91‑122 (MQ=255)
                                         ccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCtgtg                                                                                                                                       <  2:363118‑M1/32‑1 (MQ=255)
                                            gccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTgg                                                                                                                                      <  1:451920‑M1/33‑1 (MQ=255)
                                            gccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTgg                                                                                                                                      >  2:451920‑M1/88‑120 (MQ=255)
                                              aTCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACgg                                                                                                      >  2:1396755/1‑150 (MQ=255)
                                               tCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGc                                                                                                                                                         >  1:180035/1‑100 (MQ=255)
                                               tCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGc                                                                                                                                                         <  2:180035/100‑1 (MQ=255)
                                               tCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACggg                                                                                                     >  1:1344005/1‑150 (MQ=255)
                                               tCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACggg                                                                                                     <  2:606593/150‑1 (MQ=255)
                                                                  gCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCCCTGcc                                                                                  >  1:141852/1‑150 (MQ=255)
                                                                    cTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGcccc                                                                                >  2:1254468/1‑150 (MQ=255)
                                                                     tCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTTGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCg                                                                               >  1:336708/1‑150 (MQ=255)
                                                                      cAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGa                                                                              <  1:1254468/150‑1 (MQ=255)
                                                                        acccggcggtgtcggggcagcggttgttttgaagcccgtgctggcctcttcgcgggcacGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCa                                                                                                            >  2:903558‑M1/60‑118 (MQ=255)
                                                                             aTGGCCGGGTTGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACAccgcc                                                                       >  1:1265637/1‑150 (MQ=255)
                                                                                 ccGGGTGGTGACACCGATGCACCTGGATGCCGAGTCAGCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGt                                                                   <  2:336708/150‑1 (MQ=255)
                                                                                            cacCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCAAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGt                                                        <  2:346893/150‑1 (MQ=255)
                                                                                             acCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGtt                                                       >  1:33114/1‑150 (MQ=255)
                                                                                                 aTGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAAc                                                   <  1:8298/150‑1 (MQ=255)
                                                                                                 aTGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAAc                                                   <  2:1344005/150‑1 (MQ=255)
                                                                                                       cTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGtt                                             >  2:618929/1‑150 (MQ=255)
                                                                                                       cTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGtt                                             >  1:992257/1‑150 (MQ=255)
                                                                                                       cTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGtt                                             <  2:762406/150‑1 (MQ=255)
                                                                                                               cGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCccggccg                                     >  2:124273/1‑150 (MQ=255)
                                                                                                               cGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCccggccg                                     <  2:1210988/150‑1 (MQ=255)
                                                                                                                 tggcctcttcgcgggcacGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGt                                                                   >  1:903558‑M1/19‑118 (MQ=255)
                                                                                                                      tCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCt                                                                                     >  1:1149929/1‑95 (MQ=255)
                                                                                                                      tCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCt                                                                                     <  2:1149929/95‑1 (MQ=255)
                                                                                                                        cccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcagca                            >  2:1215650/1‑150 (MQ=255)
                                                                                                                        cccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcagca                            >  1:217633/1‑150 (MQ=255)
                                                                                                                         ccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcagcag                           >  1:125291/1‑150 (MQ=255)
                                                                                                                          ctctGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAgg                          >  1:58632/1‑150 (MQ=255)
                                                                                                                            cgggcacGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCCTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTg                        <  2:731720‑M1/143‑1 (MQ=255)
                                                                                                                             gggcacGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGt                                  <  1:705306‑M1/133‑1 (MQ=255)
                                                                                                                             gggcacGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGt                                  >  2:705306‑M1/7‑139 (MQ=255)
                                                                                                                                cAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCg                                                          <  1:1031191/112‑1 (MQ=255)
                                                                                                                                cAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCg                                                          >  2:1031191/1‑112 (MQ=255)
                                                                                                                                 aGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGc                   >  1:940848/1‑150 (MQ=255)
                                                                                                                                 aGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGc                   >  2:1211884/1‑150 (MQ=255)
                                                                                                                                   GCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCtg                 <  1:851214‑M1/150‑1 (MQ=255)
                                                                                                                                        CCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCtgcgtgc            >  1:1396755‑M1/1‑150 (MQ=255)
                                                                                                                                           ccACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGcc         >  2:1081903/1‑150 (MQ=255)
                                                                                                                                             acaGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGt       >  1:711041/1‑150 (MQ=255)
                                                                                                                                               aGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCa     <  2:58632/150‑1 (MQ=255)
                                                                                                                                                  ccTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCACCa  <  1:1215650/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                       
CTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAGCTGCGTGCGCCGTCACCA  >  NC_002947/5350394‑5350689

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