Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,725,040 G→A 14.2% intergenic (+25/+43) PP_5661 → / ← PP_4183 putative lipoprotein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,725,0400GA14.2% 130.7 / 12.8 57intergenic (+25/+43)PP_5661/PP_4183putative lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (24/24);  new base A (3/5);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

AAGCCGCCGAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACGG  >  NC_002947/4724894‑4725185
                                                                                                                                                  |                                                                                                                                                   
aaGCCGCCGAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGAt                                                                                                                                                   >  6:198688/1‑149 (MQ=255)
      ccgAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCACCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGAAGTC‑‑ggg                                                                                                                                           >  3:219029/1‑149 (MQ=255)
        gAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTC‑‑GGGtt                                                                                                                                         <  6:353731/149‑1 (MQ=255)
           gAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGttttt                                                                                                                                      >  1:291371/1‑149 (MQ=255)
             aaaCGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGttttttt                                                                                                                                    <  1:214238/149‑1 (MQ=255)
             aaaCGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGttttttt                                                                                                                                    >  8:245697/1‑149 (MQ=255)
              aaCGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTtg                                                                                                                                   <  2:291371/149‑1 (MQ=255)
                cGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTtgtg                                                                                                                                 >  8:132704/1‑149 (MQ=255)
                   cGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCt                                                                                                                              >  5:24291/1‑149 (MQ=255)
                   cGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTC‑‑GGGTTTTTTtgtttta                                                                                                                              <  1:209555/149‑5 (MQ=255)
                            cAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACttt                                                                                                                     >  3:154237/1‑149 (MQ=255)
                            cAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACttt                                                                                                                     >  5:268522/1‑149 (MQ=255)
                            cAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTC‑‑GGGttttttt                                                                                                                                    >  7:322460/1‑134 (MQ=255)
                            cAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTC‑‑GGGttttttt                                                                                                                                    <  8:322460/134‑1 (MQ=255)
                             aaaGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTc                                                                                                                    <  7:132704/149‑1 (MQ=255)
                             aaaGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTc                                                                                                                    <  5:198688/149‑1 (MQ=255)
                             aaaGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTC‑‑GGGTTTTTTtgtttta                                                                                                                              >  8:96907/1‑135 (MQ=255)
                             aaaGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTC‑‑GGGTTTTTTtgtttta                                                                                                                              <  7:96907/139‑5 (MQ=255)
                               aGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAACCCGACGATGTC‑‑GGGt                                                                                                                                          <  4:275037/124‑1 (MQ=255)
                               aGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAACCCGACGATGTC‑‑GGGt                                                                                                                                          >  3:275037/1‑124 (MQ=255)
                                   cAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGtttt                                                                                                                                       >  3:11004/1‑124 (MQ=255)
                                   cAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGtttt                                                                                                                                       <  4:11004/124‑1 (MQ=255)
                                    aGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTaa                                                                                                              <  7:245697/148‑1 (MQ=255)
                                      gCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCg                                                                                                           >  3:230109/1‑149 (MQ=255)
                                         gAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGttt                                                                                                        <  8:238867/149‑1 (MQ=255)
                                             ccGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGt                                                                                                    <  6:268522/149‑1 (MQ=255)
                                                 gCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTgg                                                                                                 <  6:24291/148‑1 (MQ=255)
                                                    gAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAaac                                                                                             <  2:331301/149‑1 (MQ=255)
                                                     aaGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGttt                                                                                                        >  7:111590/1‑137 (MQ=255)
                                                     aaGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGttt                                                                                                        <  8:111590/137‑1 (MQ=255)
                                                         gagCCGGAGTAGCCGGAGCGGCTTCTTCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGttttttt                                                                                                                                    <  2:68636/103‑1 (MQ=37)
                                                         gagCCGGAGTAGCCGGAGCGGCTTCTTCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGttttttt                                                                                                                                    >  1:68636/3‑105 (MQ=37)
                                                           gCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGTCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAg                                                                                      >  3:55864/1‑149 (MQ=255)
                                                                   agCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCtt                                                                              <  4:230109/149‑1 (MQ=255)
                                                                    gCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTATCCTTg                                                                             >  7:142816/1‑149 (MQ=255)
                                                                        cTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTTGAAACAACCGAGCTTTCCTTGGCCg                                                                         >  7:77066/1‑149 (MQ=255)
                                                                         agcGGCTTCTTCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGTCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGcc                                                                          <  2:193633/145‑1 (MQ=255)
                                                                           cGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAATCATTAATTGGCTAAAAAAAAAAAACCCGACTATGTCTACGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGcc                                                                          <  1:111661/147‑1 (MQ=255)
                                                                                  tCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAg                                                               >  3:286744/1‑149 (MQ=255)
                                                                                    cGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGTCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAg                                                               <  4:55864/147‑1 (MQ=255)
                                                                                                cacCGGCAGAACCTGCCAAGCAGTAATTGGATAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtcggtcg                                                 <  8:21591/149‑1 (MQ=255)
                                                                                                              gCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCt                                                                      <  6:293419/114‑1 (MQ=255)
                                                                                                              gCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCt                                                                      >  5:293419/1‑114 (MQ=255)
                                                                                                                 aaGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCAt                                    <  7:169797/145‑1 (MQ=255)
                                                                                                                 aaGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCAt                                    >  8:169797/1‑145 (MQ=255)
                                                                                                                     aGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTTGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGa                            <  8:77066/149‑1 (MQ=255)
                                                                                                                     aGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATc                              >  1:356383/1‑147 (MQ=255)
                                                                                                                     aGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATc                              <  2:356383/147‑1 (MQ=255)
                                                                                                                      gTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTc                                                                                                       >  1:170675/1‑73 (MQ=255)
                                                                                                                      gTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTc                                                                                                       <  2:170675/73‑1 (MQ=255)
                                                                                                                               tAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAAACGAGCTTTCCTTTGCCGGCTCCGACAGTTCACTTTCGGTCGGTGTGAACACCATCACATAGAGTACATACg                   >  6:165216/1‑148 (MQ=255)
                                                                                                                                     aaaaaaaCCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGaa                 <  3:130175/144‑1 (MQ=255)
                                                                                                                                     aaaaaaaCCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGaa                 >  4:130175/1‑144 (MQ=255)
                                                                                                                                     aaaaaaaCCCGACATCGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATg               <  2:176482/146‑1 (MQ=255)
                                                                                                                                     aaaaaaaCCCGACATCGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATg               >  1:176482/1‑146 (MQ=255)
                                                                                                                                            cccGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGa     <  4:154237/149‑1 (MQ=255)
                                                                                                                                             ccGACGAAGTC‑‑GGATTTTTTTGTTACTTCTAACTATCGATTAATCGTTTCAGTTGCAAACAACAGAGCAATCGTTGGACGGCTCAGACAGTACACTTTAGGACGGTGTGAACACCATAACATAGATTACAAACCAAAGGAACGCTCTAc    >  2:57070/1‑149 (MQ=255)
                                                                                                                                               gacgaTGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACgg  >  8:279325/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAGCCGCCGAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTC‑‑GGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACGG  >  NC_002947/4724894‑4725185

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