Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,740,474 C→A 21.5% intergenic (+39/+32) hutF → / ← PP_5037 formiminoglutamate deiminase/lipocalin family lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,740,4740CA21.5% 111.2 / 24.4 55intergenic (+39/+32)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (18/25);  new base A (6/6);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

ATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGT  >  NC_002947/5740332‑5740617
                                                                                                                                                |                                                                                                                                               
aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGaa                                                                                                                                           <  2:257603/149‑1 (MQ=255)
 tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                          <  2:100575/149‑1 (MQ=255)
  gtgAGGGTGACTGGGCACTGGGTGGTTCGCCAGGGGCGGCATGCCCAGGAAGAGGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTTAGCGTACCGGGG‑TGGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                          <  2:247755/148‑1 (MQ=37)
   tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAggg                                                                                                                                        <  3:54414/149‑1 (MQ=255)
   tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAgg                                                                                                                                         <  2:244165/148‑1 (MQ=255)
     aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGGTTTGGGTTG‑AGCCCCTAGCCTGGAAGCCCGCCAAAGGtt                                                                                                                                      >  2:356442/1‑149 (MQ=255)
     aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggc                                                                                                                                       >  1:191283/1‑147 (MQ=38)
     aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggc                                                                                                                                       <  2:191283/148‑2 (MQ=38)
      ttttGAATGGGGAATGGGTGTTGCGTCAGGGGCGGCATGCGCAGGTAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGTCTGAGCGTACCGGGG‑TTGCGTTGCAACGCTTCGCGGGGTTGCACGCGAAGGGGTg                                                                                                                                     <  8:373868/146‑1 (MQ=255)
      tGGTGAATGGGCAATGGGGGGTGCGCCAGGGGCGGCAAGCCCATGAAGTTGTAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTTCTGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                      <  1:304224/149‑1 (MQ=255)
       ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTTCGTTTGAGCCCATACCCTGGAAGCCCCGGAAAGGTTg                                                                                                                                     >  8:261094/1‑148 (MQ=255)
       ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAggg                                                                                                                                        <  1:349803/145‑1 (MQ=255)
       ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAggg                                                                                                                                        >  2:349803/1‑145 (MQ=255)
         tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGGGGTTGAGCCCCTGGCGGGGCAGCCCGCGCAGGGGTggca                                                                                                                                  >  5:345102/1‑146 (MQ=255)
         tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                           <  8:356195/140‑1 (MQ=255)
         tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                           >  7:356195/1‑140 (MQ=255)
           aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTTGAGCCCATAGCGTGGAAGCCCGCGAAGGGTTtgaacg                                                                                                                                >  2:21489/1‑149 (MQ=255)
           aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCCTTGCAGCCCATAGCGTGCCACCCCGCGCAGGGTTgccacg                                                                                                                                >  6:173765/1‑149 (MQ=255)
                       gtggggCGCCAGGGGCGGCATGCCCAGGAAGAGGAACGCGGGCGGGCGTTTGCCAAGGTGGTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTGGCAGCCCTGCGCGTGCAACCCCGCGAAGGGTTGCAACGCAACCCCaaaaa                                                                                                                    <  7:93242/149‑1 (MQ=255)
                        tcgtccgcCAGGGGCGGCATGCCCAGGAAGATGAACGCGGCCGGGCGTTTGCCAAGGTGTTGCCCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCaaaaaa                                                                                                                   <  5:198105/144‑1 (MQ=255)
                              gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTTCAGCCCATAGCGTGGAAGGCCGCGCAGGGGTGCAACGCCACCCCACAAAACTCAGt                                                                                                             >  6:360659/1‑149 (MQ=255)
                              gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                          >  1:135291/1‑120 (MQ=38)
                              gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                          <  2:135291/120‑1 (MQ=38)
                               ccAGGGGCCGCATGCCCACGAAGATGAAAGCGGGCGGGCGTTTGCCAAGCTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTAACCGGGCAACGCCCCCCAGGGTTGCAAAGCCACCCCAAAAAACtagctg                                                                                                            >  8:283920/1‑144 (MQ=255)
                                        gcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCCATTCAGCCCATACCGTGGCACCCCCCCAAGGGTTGCAACCCAACCCCAAAAAACAGAGAGCAtttttt                                                                                                    >  6:183937/1‑147 (MQ=255)
                                             cccAGGAAGATGAAAGGGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAtt                                                                                               <  7:261094/148‑1 (MQ=255)
                                             cccAGGAAGATGAAAGCGGGCGGTCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCTCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAttt                                                                                              <  5:183937/149‑1 (MQ=255)
                                             cccAGGAAGATGAAAGCGGGCGGGCTTTTGCCAAGGTGTTGCGCCTGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAttt                                                                                              <  5:173765/149‑1 (MQ=255)
                                             cccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAttt                                                                                              <  1:21489/149‑1 (MQ=255)
                                              ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTg                                                                                             <  6:377057/149‑1 (MQ=255)
                                               cAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTTCCGCCCCTCGCCGGGCAGCCCCCGCCCGGGTGCCCCGCCACCCCAAAAAACACAGTGCCATTTTAATAttttc                                                                                            >  7:4128/1‑147 (MQ=255)
                                               cAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc                                                                                            <  5:360659/149‑1 (MQ=255)
                                                 ggAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTg                                                                                          <  6:374808/149‑1 (MQ=255)
                                                  taaGATCAAAGCGGGGGGGCGTTTGCCATGGTGTTGCGCCTGG‑TTTTGGGCTGAGCGTACGGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGa                                                                                         <  6:131083/148‑1 (MQ=255)
                                                         aaaGCGGGCGGGCGTTTGCAAAGGTGTTGCGCCAGC‑TTTTGGGCTTAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAc                                                                                   <  1:8980/148‑1 (MQ=255)
                                                               ggcgggcgTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCAAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgcc                                                                            <  1:356442/149‑1 (MQ=255)
                                                               ggcgggcgTTTGCCAAGGTGTGGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TGGCGTTGCAGCCCTTCGCGGCCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgcc                                                                            <  3:360534/149‑1 (MQ=255)
                                                                                gtgtTGCGCCAGC‑TTGTGGGCTGCGCGTACCGGGG‑TTGGGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                            <  8:4128/148‑1 (MQ=255)
                                                                                     gcgcCAGC‑TTTTGGGCTGCGCGTACCGGGG‑TGGCGTGGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCAt                                                       <  7:283920/148‑1 (MQ=255)
                                                                                      cgcCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc                                                     <  3:40293/149‑1 (MQ=255)
                                                                                                                   g‑ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCACTGGACAGCAGGt                         >  2:34171/1‑149 (MQ=255)
                                                                                                                       gcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtct                      >  3:362277/1‑149 (MQ=255)
                                                                                                                         gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcg                                          <  8:360020/127‑1 (MQ=255)
                                                                                                                         gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcg                                          >  7:360020/1‑127 (MQ=255)
                                                                                                                             cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt                >  7:283667/1‑149 (MQ=255)
                                                                                                                             cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt                >  1:224756/1‑149 (MQ=255)
                                                                                                                              aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc               >  4:211815/3‑149 (MQ=38)
                                                                                                                              aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc               >  3:63730/1‑149 (MQ=255)
                                                                                                                               acccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGACGCGCCCTGGACAGCAGGTTCTCACGCTcc              >  7:367206/2‑149 (MQ=39)
                                                                                                                                  cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg           >  4:18890/1‑149 (MQ=39)
                                                                                                                                   ttCGCGGGCTTGCACGCGAAGGGTCGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg          >  2:365499/1‑149 (MQ=255)
                                                                                                                                     cgcgTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc         <  7:39653/148‑1 (MQ=255)
                                                                                                                                           cAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgt  <  4:362277/149‑1 (MQ=255)
                                                                                                                                              gCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAca                               <  8:317463/117‑1 (MQ=255)
                                                                                                                                              gCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAca                               >  7:317463/1‑117 (MQ=255)
                                                                                                                                                |                                                                                                                                               
ATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGT  >  NC_002947/5740332‑5740617

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