Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 188,613 T→C 12.5% L241P (CTG→CCG)  PP_0165 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947188,6130TC12.5% 134.5 / 10.7 56L241P (CTG→CCG) PP_0165GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (27/22);  new base C (3/4);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

CAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGAC  >  NC_002947/188478‑188761
                                                                                                                                       |                                                                                                                                                    
cAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcc                                                                                                                                         <  6:216040/149‑2 (MQ=255)
cAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcc                                                                                                                                         <  2:44192/149‑2 (MQ=255)
 agagTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcct                                                                                                                                        <  3:284691/149‑3 (MQ=255)
  gagTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGaa                                                                                                                                       >  5:175572/1‑149 (MQ=255)
         cacgcacgCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcctggt                                                                                                                                     >  7:160523/1‑139 (MQ=255)
         cacgcacgCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcctggt                                                                                                                                     <  8:160523/144‑6 (MQ=255)
          acgcacgcCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACGCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACcgccgc                                                                                                                                >  7:268121/1‑148 (MQ=255)
           cgcacgcCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTa                                                                                                                              >  8:309278/1‑149 (MQ=255)
                gcCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGAAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTAtttt                                                                                                                          >  8:247180/1‑148 (MQ=255)
                            ccgcGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTg                                                                                                              >  3:147626/1‑148 (MQ=255)
                                      cAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGTCCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCggg                                                                                                    <  2:135875/148‑1 (MQ=255)
                                      cAAGCGAGGACCCAGATGGTCGCGAACCTCAACGCCCGGTTCCCCGAGCAGGCCGAGCGCAGGGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCggg                                                                                                    <  8:268121/148‑1 (MQ=255)
                                                aCCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCACCCTg                                                                                         >  6:332946/1‑149 (MQ=255)
                                                 ccAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTg                                                                                         >  2:278961/1‑148 (MQ=255)
                                                 ccAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTg                                                                                         >  2:55902/1‑148 (MQ=255)
                                                 ccAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGAAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTgg                                                                                        <  7:247180/149‑1 (MQ=255)
                                                   aGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGaa                                                                                      <  7:282436/149‑1 (MQ=255)
                                                    gATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTgg                                                                                        >  2:344002/1‑146 (MQ=255)
                                                       ggTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAAccgcc                                                                                                                                  >  2:299478/1‑101 (MQ=255)
                                                       ggTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAAccgcc                                                                                                                                  <  1:299478/101‑1 (MQ=255)
                                                          cGAGAAACTCAAAGCCCTGTTCACTGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGt                                                                                                   >  2:236223/1‑129 (MQ=255)
                                                          cGAGAAACTCAAAGCCCTGTTCACTGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGt                                                                                                   <  1:236223/129‑1 (MQ=255)
                                                           gagaAACTCAAAGCCCTGTTCACCGAGGAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAcgcgcg                                                                              >  1:10672/1‑149 (MQ=255)
                                                           gagaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAcgcgcg                                                                              >  7:65888/1‑149 (MQ=255)
                                                            agaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAcgcgcgc                                                                             >  5:292013/1‑149 (MQ=255)
                                                              aaaCTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACgcgcgcgc                                                                           <  7:149526/149‑1 (MQ=255)
                                                               aaCTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACgcgcgcgc                                                                           <  1:344002/148‑1 (MQ=255)
                                                                aCTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCg                                                                         <  7:254201/149‑1 (MQ=255)
                                                                 cTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCgg                                                                        >  1:175966/1‑149 (MQ=255)
                                                                    aaaGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGt                                                                       <  1:278961/147‑1 (MQ=255)
                                                                    aaaGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGt                                                                       <  6:175572/147‑1 (MQ=255)
                                                                    aaaGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTa                                                                     >  5:204011/1‑149 (MQ=255)
                                                                       gCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACct                                                                  <  4:147626/149‑1 (MQ=255)
                                                                        cccTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACctg                                                                 <  1:79256/149‑1 (MQ=255)
                                                                          cTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACctgct                                                               >  6:322875/1‑149 (MQ=255)
                                                                           tGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACctgctg                                                              >  6:250457/1‑149 (MQ=255)
                                                                            gTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGt                                                             >  7:292880/1‑149 (MQ=255)
                                                                                aCCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTtgct                                                         >  3:259992/1‑149 (MQ=255)
                                                                                   gAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCg                                                      >  8:81105/1‑149 (MQ=255)
                                                                                         gCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCa                                                <  1:55902/149‑1 (MQ=255)
                                                                                                     gAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCctggctgg                                    >  8:112198/1‑149 (MQ=255)
                                                                                                      aGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCctggctgg                                    <  7:81105/148‑1 (MQ=255)
                                                                                                      aGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCtggctggc                                   <  7:309278/149‑1 (MQ=255)
                                                                                                        cGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTTCTGCGGGTTCAggg                                             <  5:291965/137‑1 (MQ=255)
                                                                                                        cGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCCCCGGTTACCTGCTGTTGCTGCGGGTTCAggg                                             >  6:291965/1‑137 (MQ=255)
                                                                                                              cgtgcgGAGTCCTACCATGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgc                           <  2:260562/149‑1 (MQ=255)
                                                                                                               gtgcgGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgcc                          >  3:35960/1‑149 (MQ=255)
                                                                                                                tgcgGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGGAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgccc                         <  2:175966/149‑1 (MQ=255)
                                                                                                                tgcgGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgccc                         >  3:287374/1‑149 (MQ=255)
                                                                                                                      gTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTg                   <  8:65888/149‑1 (MQ=255)
                                                                                                                           aCCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcggc              >  3:706/1‑149 (MQ=255)
                                                                                                                            ccAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcggc              >  5:337408/1‑148 (MQ=255)
                                                                                                                             ttcgACAGCCCGGTCAGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcggc              >  2:210423/4‑147 (MQ=255)
                                                                                                                                gACAGCCCGGTCAGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCg         >  4:264317/1‑149 (MQ=255)
                                                                                                                                 aCAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGt        <  2:10672/149‑1 (MQ=255)
                                                                                                                                     ccTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTAcc     <  8:292880/148‑1 (MQ=255)
                                                                                                                                       tGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGAc  >  4:125809/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                                    
CAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTCGAACCGCCGCTATTTTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGAC  >  NC_002947/188478‑188761

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: ‑