Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 196,390 T→C 19.6% D632D (GAT→GAC)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947196,3900TC19.6% 104.1 / 22.3 56D632D (GAT→GAC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (19/26);  new base C (4/7);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

CCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCAC  >  NC_002947/196242‑196533
                                                                                                                                                     |                                                                                                                                                 
cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGAt                                                                                                                                                   >  7:264378/1‑149 (MQ=34)
cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGAt                                                                                                                                                   >  3:316120/1‑149 (MQ=34)
   tGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTc                                                                                                                                                >  6:93315/1‑149 (MQ=34)
     aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTc                                                                                                                                                <  8:264378/147‑1 (MQ=34)
     aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTc                                                                                                                                                <  2:48619/147‑1 (MQ=34)
     aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGACTTcc                                                                                                                                               <  1:204511/148‑1 (MQ=255)
       ggCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGa                                                                                                                                            <  5:93315/149‑1 (MQ=34)
          gAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGAAAAACCCAGGGCAGCGT‑CGATTTCCAGAccc                                                                                                                                         >  4:232463/1‑149 (MQ=21)
                gTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCg                                                                                                                                       >  8:298773/1‑145 (MQ=34)
                gTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCg                                                                                                                                       <  7:298773/145‑1 (MQ=34)
                gTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCg                                                                                                                                    >  4:308163/1‑148 (MQ=34)
                gTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGACTTCCAGACCCCGGCg                                                                                                                                    >  3:207648/1‑148 (MQ=255)
                    gTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTAGCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGa                                                                                                                               <  1:88460/149‑1 (MQ=21)
                    gTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGACTTCCAGACCCCGGCGAATGa                                                                                                                               >  3:127656/1‑149 (MQ=255)
                     tCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGACTTCCAGACCCCGGCGAATGAc                                                                                                                              <  3:169766/149‑1 (MQ=255)
                       aCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGACTTCCAGACCCCGGCGAATGACGt                                                                                                                            <  3:243995/149‑1 (MQ=255)
                        ccGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTc                                                                                                                           <  3:308163/149‑1 (MQ=34)
                         cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTc                                                                                                                           >  8:122867/1‑148 (MQ=34)
                           aaGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTac                                                                                                                        <  4:291652/149‑1 (MQ=34)
                           aaGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGACTTCCAGACCCCGGCGAATGACGTCTac                                                                                                                        <  8:191570/149‑1 (MQ=255)
                             ggcggcGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                                                      <  3:232463/149‑1 (MQ=34)
                             ggcggcGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                                                      <  2:103012/149‑1 (MQ=34)
                                ggcgTGATCACATACACCGTCACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCTGCGAATGACGTCTAcaacaa                                                                                                                   <  7:132569/149‑1 (MQ=12)
                                  cgTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGACTTCCAGACCCCGGCGAATGACGTCTACAACAACg                                                                                                                 <  3:60587/149‑1 (MQ=255)
                                         aCTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGa                                                                                                           <  1:166209/148‑1 (MQ=34)
                                         aCTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGa                                                                                                           <  8:143490/148‑1 (MQ=255)
                                          cTTTCACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGATAGTGACCCTGTCCAACGGCCATACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGAGCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                         <  8:28154/149‑1 (MQ=255)
                                          cTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                         <  4:316120/149‑1 (MQ=34)
                                             acacCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGt                                                                                                       >  3:143729/1‑148 (MQ=255)
                                             acacCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGACTTCCAGAACCCGGCGAATGACGTCTACAACAACGGTTCGACTGt                                                                                                       >  6:132150/1‑148 (MQ=11)
                                              cacCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTc                                                                                                      <  2:322582/148‑1 (MQ=34)
                                               acCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCa                                                                                                     <  1:79558/149‑1 (MQ=21)
                                                ccGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑gc                                                                                                  <  6:197631/149‑2 (MQ=255)
                                                 cGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑gtg                                                                                                 >  8:239901/1‑149 (MQ=35)
                                                     aCCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCa                                                                                              <  7:239901/148‑1 (MQ=37)
                                                               aTCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAAc                                                                                   <  1:210760/149‑1 (MQ=37)
                                                                tCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACg                                                                                  >  5:61564/1‑149 (MQ=37)
                                                                   tGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCa                                                                               >  2:13296/1‑149 (MQ=37)
                                                                       cAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACcgg                                                                           >  2:75215/1‑149 (MQ=39)
                                                                           cGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGt                                                                       <  6:61564/149‑1 (MQ=39)
                                                                            gCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGt                                                                       >  8:268914/1‑148 (MQ=39)
                                                                              cGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACcggcgg                                                                        <  6:158860/145‑1 (MQ=39)
                                                                              cGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACcggcgg                                                                        >  5:158860/1‑145 (MQ=39)
                                                                                           ccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGa                                                      >  5:238403/1‑149 (MQ=37)
                                                                                                ccAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCACGGCAGCGG‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGGTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGAcccc                                                  <  1:75215/149‑1 (MQ=255)
                                                                                                 cAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGa                                                              <  1:204285/135‑1 (MQ=37)
                                                                                                 cAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGa                                                              >  2:204285/1‑135 (MQ=37)
                                                                                                  aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCa                                                                                                     <  2:220911/97‑1 (MQ=33)
                                                                                                  aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCa                                                                                                     >  1:220911/1‑97 (MQ=33)
                                                                                                  aaCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGACCCCGaa                                               >  5:26593/1‑149 (MQ=37)
                                                                                                             cATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAACTTTCGAGCAACTGACCCCGAATCCGACGCCg                                     <  1:36314/149‑1 (MQ=255)
                                                                                                                              gCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGATGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAAc                   >  5:56398/1‑149 (MQ=255)
                                                                                                                              gCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAAc                   <  7:122867/149‑1 (MQ=39)
                                                                                                                              gCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAAc                   <  1:13296/149‑1 (MQ=39)
                                                                                                                                          gCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAt       >  3:88902/1‑149 (MQ=39)
                                                                                                                                           cAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATg      >  1:317737/1‑149 (MQ=39)
                                                                                                                                             gCGT‑CGATTTCCAGACCCCGGCGAATGATGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGcc    <  6:56398/149‑1 (MQ=255)
                                                                                                                                                t‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGccac  >  6:212113/1‑148 (MQ=39)
                                                                                                                                                     |                                                                                                                                                 
CCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGT‑CGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCA‑GTGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCAC  >  NC_002947/196242‑196533

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