Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,032,797 T→C 7.3% intergenic (+61/+75) dppA‑IV → / ← PP_5284 dipeptide transport protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,032,7970TC 7.3% 149.4 / 2.2 57intergenic (+61/+75)dppA‑IV/PP_5284dipeptide transport protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/30);  new base C (2/2);  total (24/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 = 5.84e-01

GCGACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATC  >  NC_002947/6032650‑6032937
                                                                                                                                                    |                                                                                                                                            
gcgACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGtt                                                                                                                                             <  3:225344/149‑1 (MQ=255)
gcgACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGtt                                                                                                                                             >  2:240042/1‑149 (MQ=255)
 cgACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCCGCCTGGCAGTTg                                                                                                                                            >  1:22463/1‑149 (MQ=255)
  gACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGa                                                                                                                                           <  6:185073/149‑1 (MQ=255)
    ccTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACg                                                                                                                                         <  5:177805/149‑1 (MQ=255)
     cTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGc                                                                                                                                        >  3:261272/1‑149 (MQ=255)
        gctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTa                                                                                                                                      <  1:193231/148‑1 (MQ=255)
        gctgcGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTa                                                                                                                                      <  8:134889/148‑1 (MQ=255)
              ccAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGt                                                                                                                               <  8:151173/149‑1 (MQ=25)
              ccAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGt                                                                                                                               <  3:81415/149‑1 (MQ=255)
               cAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTa                                                                                                                              <  1:240042/149‑1 (MQ=255)
                aaGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAg                                                                                                                             <  6:250270/149‑1 (MQ=255)
                aaGGCGGCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAg                                                                                                                             <  7:67045/149‑1 (MQ=255)
                     gTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCAAGCCCGGTAGGAGCg                                                                                                                        >  1:110626/1‑149 (MQ=255)
                      tCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCgg                                                                                                                       <  1:263321/149‑1 (MQ=255)
                          ggggTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGtt                                                                                                                    >  7:225807/1‑148 (MQ=255)
                             gTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCGGTAGGAGCGGGGTTAcc                                                                                                                >  8:34786/1‑149 (MQ=255)
                             gTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCCCTCCTAAAGGGGTAGCGCCAGCCTGGGAGTTGACGCTACCCCTGTAGGAGCGGGTTTaac                                                                                                                >  8:63018/1‑147 (MQ=255)
                                ccAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCcgc                                                                                                             <  7:34786/149‑1 (MQ=255)
                                   ggTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGCCGCTACCCCTGTAGGAGCGGGTTTACCCGCGa                                                                                                           >  1:52560/1‑148 (MQ=255)
                                    gcgcgCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCACGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGaaa                                                                                                         <  8:290320/145‑1 (MQ=255)
                                    gTGACCCCGTTCGGGCGGCTGTATTTCAGCACGGTGACGGTGACCAAGTCAGCCTGTCCCGGCGTCTTCGCGGGTAACCCCGCTCCTACCGGGGTCGCGCCAGCCTGGCTGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGaaa                                                                                                         <  2:256161/149‑1 (MQ=37)
                                                 ggcggcTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGTCCCGTACTGa                                                                                            >  5:172717/1‑149 (MQ=255)
                                                   ggggtgtaTTTCAGAATGTTTACGTTTACCATTTAAGCCTGTCCGGCCCTCTTCGCGGGTAAACCGGCTCCTACCGGGTTAGCGCCAGCTTGGATGTTGACGCTACCCGTGTAGTAGGGGGTTTACCCGCGAAAGGGCCCGTACTGAcc                                                                                          <  6:172717/142‑1 (MQ=19)
                                                      cTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAACGGCCCGTACTGACCCt                                                                                        >  7:290834/1‑148 (MQ=255)
                                                       tGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCCGTACTGACCCtat                                                                                      <  4:261272/149‑1 (MQ=255)
                                                         gATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAAAGGGGTAGCGCCAGCCGGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGGAAggg                                                                                                      >  7:172465/1‑131 (MQ=255)
                                                           tttCAGCAAGGTGACGGTGACCAAGAAACCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCAACCGGGGAAGCGCCAGCCTGGCAGTTGCCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAgg                                                                                  <  7:63018/149‑1 (MQ=38)
                                                               aGCAAGGTGACGGTTACCAAGTACGCGTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTTGCGCCAGCCTGGCTGTTGACGCTACCCCTGTAGTAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGctct                                                                              <  6:126406/149‑1 (MQ=38)
                                                                  atggtgacggtcacCAAGTACGCCTGTACGGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAg                                                                           <  2:52560/147‑1 (MQ=255)
                                                                       gacggtgacCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTc                                                                      >  1:50513/1‑149 (MQ=255)
                                                                              acCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAg                                                                <  2:237789/148‑1 (MQ=255)
                                                                                cAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCGTGGCAGTTAACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGAcc                                                             >  1:99802/1‑149 (MQ=255)
                                                                                cAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCGACCCCGGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGGACGGACCCTATAAGGCTCTCAGTACTACCACACAcc                                                             >  4:214318/1‑149 (MQ=37)
                                                                                 aaGTATGCCTGTACCGGCCTCTTCGCGGTTAAACCCGCTCCTACCGGGGTAGCGCCAGCTTGGCTGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACca                                                            <  2:90771/149‑1 (MQ=38)
                                                                                      aGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATc                                                       <  2:50513/149‑1 (MQ=255)
                                                                                       cccTGTACCGCCCGCTTCGCTGTTAAACCCTATCCTACAGGGGTAGCCCCAGCCTTGCCGATGACGCTACCCCTTTCGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAACGCTCTCAGTACTCCCTCAGACCACAATcc                                                      <  4:214228/148‑1 (MQ=255)
                                                                                         cTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGTGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCAc                                                    <  8:225807/149‑1 (MQ=255)
                                                                                            tcccGGCCTCTTCGCGGCTAAACCCGCTCCTACCGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAg                                                                                                        <  5:192907/92‑1 (MQ=255)
                                                                                                     tttCGCGGGTAAACCCGCTCCTACAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGCAGGGGCCCGTACTGACCCTATAAGGCTCTCAGGACTCCCTCAGACCACCATCCACCCATGCTCggc                                         >  5:178820/2‑146 (MQ=14)
                                                                                                          cgGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGAc                                                                                           <  8:238297/93‑1 (MQ=255)
                                                                                                          cgGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGAc                                                                                           >  7:238297/1‑93 (MQ=255)
                                                                                                             gTAAACCCGCTCCTACAGGGGTAGCGCCAGCTTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACaa                                <  2:99802/149‑1 (MQ=255)
                                                                                                                   ccGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTc                          >  5:329935/1‑149 (MQ=255)
                                                                                                                        ccTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTTCCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGAcc                                        <  5:270135/130‑1 (MQ=255)
                                                                                                                        ccTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTTCCCCTGTAGGAGCGGCTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCAATGCTCGAcc                                        >  6:270135/1‑130 (MQ=255)
                                                                                                                           aCAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTc                                                                      <  7:177603/97‑1 (MQ=14)
                                                                                                                           aCAGGGGTAGCGTCAACTGCCAGGCTGGCGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTc                                                                      >  8:177603/1‑97 (MQ=14)
                                                                                                                               gggTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCCGACCACACTCCACCCATGCTCGACCATCGACAACGGCACCCCATCACCa                >  5:235370/1‑147 (MQ=255)
                                                                                                                                     cgcCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGTCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAatgatgat        >  3:234149/1‑149 (MQ=255)
                                                                                                                                       ccAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTcc                         >  2:190881/1‑130 (MQ=255)
                                                                                                                                       ccAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTcc                         <  1:190881/130‑1 (MQ=255)
                                                                                                                                       ccAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                  >  5:241247/1‑137 (MQ=255)
                                                                                                                                       ccAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCAc                  <  6:241247/137‑1 (MQ=255)
                                                                                                                                           ccTGGCCGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATc  <  3:214318/149‑1 (MQ=255)
                                                                                                                                           ccTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAATGGCCCGTACTGACCCTATAAGGctc                                                                               <  7:177863/72‑1 (MQ=255)
                                                                                                                                           ccTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAATGGCCCGTACTGACCCTATAAGGctc                                                                               >  8:177863/1‑72 (MQ=255)
                                                                                                                                             tGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCa                   >  5:129333/1‑130 (MQ=255)
                                                                                                                                             tGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCa                   <  6:129333/130‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GCGACCTTGCTGCGCCAAGGCGTCGAGGGGTACCAGGTGAGCCCGTTCGGGCGGCTGGATTTCAGCAAGGTGACGGTGACCAAGTAAGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGTAGCGCCAGCCTGGCAGTTGACGCTACCCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCCGTACTGACCCTATAAGGCTCTCAGTACTCCCTCAGACCACAATCCACCCATGCTCGACCATCGACAACGGCTCCCCATCACCAATGATGATGTGATC  >  NC_002947/6032650‑6032937

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