Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,461,375 T→G 7.0% V125G (GTG→GGG)  lolC → lipoprotein releasing system transmembrane protein LolC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,461,3750TG 7.0% 171.0 / 4.8 57V125G (GTG→GGG) lolClipoprotein releasing system transmembrane protein LolC
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base G (2/2);  total (30/27)
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.08e-01

GCAGGCCCTTGCCCAACAAGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCTCACCTTCGTCGCCCCCGAGGTCAGCGT  >  NC_002947/2461232‑2461522
                                                                                                                                               |                                                                                                                                                   
gCAGGCCCTTGCCCAACAAGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGa                                                                                                                                                >  6:278540/1‑149 (MQ=255)
gCAGGCCCTTGCCCAACAAGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGa                                                                                                                                                <  1:60880/149‑1 (MQ=255)
 cAGGCCCTTGCCCAACAAGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGa                                                                                                                                                >  5:60140/1‑148 (MQ=255)
  aGGCCCTTGCCCAACAAGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGat                                                                                                                                               >  7:142767/1‑148 (MQ=255)
    gCCCTTGCCCAACAAGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCTTTCACCCAGATGCAGGGCCTGCTGACCCATGATGTCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGatcat                                                                                                                                            <  8:252166/149‑1 (MQ=255)
    gCCCTTGCCCAACAAGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGatcat                                                                                                                                            >  6:67900/1‑149 (MQ=255)
          gCCCAACAAGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGAc                                                                                                                                        <  3:281290/147‑1 (MQ=255)
          gCCCAACAAGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGAc                                                                                                                                        <  3:153598/147‑1 (MQ=255)
              aacaaGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTc                                                                                                                                  >  4:176871/1‑149 (MQ=255)
              aacaaGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTc                                                                                                                                  <  8:142767/149‑1 (MQ=255)
                 aaGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGTTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTc                                                                                                                               <  5:226132/149‑1 (MQ=255)
                 aaGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTc                                                                                                                               <  5:359163/149‑1 (MQ=255)
                 aaGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTc                                                                                                                               <  3:272508/149‑1 (MQ=255)
                   gTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGatcat                                                                                                                                            >  7:187118/1‑134 (MQ=255)
                   gTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTccc                                                                                                                             >  2:136552/1‑148 (MQ=255)
                           gAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGgcc                                                                                                                     >  3:297018/1‑149 (MQ=255)
                                          aGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTgg                                                                                                      >  4:180136/1‑149 (MQ=255)
                                              gTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGatc                                                                                                                                              <  7:351474/105‑1 (MQ=255)
                                              gTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGatc                                                                                                                                              >  8:351474/1‑105 (MQ=255)
                                                  cgcCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGc                                                                                              <  3:180136/149‑1 (MQ=255)
                                                   gcCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGc                                                                                              >  1:235194/1‑148 (MQ=255)
                                                   gcCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGc                                                                                              >  5:254130/1‑148 (MQ=255)
                                                      gTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGt                                                                                          >  7:346570/1‑149 (MQ=255)
                                                      gTTCACCCAGATGCAGGGCCTGCTGACCCATGAGGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGt                                                                                          <  7:343054/149‑1 (MQ=255)
                                                       ttCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTg                                                                                         >  7:20634/1‑149 (MQ=255)
                                                          aCCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTggg                                                                                       <  5:278540/148‑1 (MQ=255)
                                                            ccAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGca                                                                                    >  8:226347/1‑149 (MQ=255)
                                                            ccAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGca                                                                                    >  7:180293/1‑149 (MQ=255)
                                                                            ctgACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTggg                                                                                       <  4:207059/130‑1 (MQ=255)
                                                                            ctgACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTggg                                                                                       >  3:207059/1‑130 (MQ=255)
                                                                                 ccATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAggcgg                                                               >  8:7414/1‑149 (MQ=255)
                                                                                  cATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAAGTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAggcggc                                                              >  1:118640/1‑149 (MQ=255)
                                                                                               tGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGc                                                 <  8:346570/149‑1 (MQ=255)
                                                                                                gCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCg                                                >  2:272588/1‑149 (MQ=255)
                                                                                                  aGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCCGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGtt                                              <  1:136552/149‑1 (MQ=255)
                                                                                                    aaGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGc                                            <  4:297018/149‑1 (MQ=255)
                                                                                                       gTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCAt                                         <  5:67900/149‑1 (MQ=255)
                                                                                                           tgctCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGc                                     >  2:356442/1‑149 (MQ=255)
                                                                                                              tCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCg                                                                                                 >  3:312363/1‑86 (MQ=255)
                                                                                                              tCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCg                                                                                                 <  4:312363/86‑1 (MQ=255)
                                                                                                                aaTGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGAc                                  >  1:28667/1‑147 (MQ=255)
                                                                                                                aaTGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGAc                                  >  5:123381/1‑147 (MQ=255)
                                                                                                                 aTGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAg                               <  8:20634/149‑1 (MQ=255)
                                                                                                                   atgaTCGACACCTCGCGCTCGCGAGCGGGGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCAt                                         <  4:303110/134‑1 (MQ=38)
                                                                                                                   atgaTCGACACCTCGCGCTCGCGAGCGGGGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCAt                                         >  3:303110/4‑137 (MQ=38)
                                                                                                                    gCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCt                             <  6:254130/148‑1 (MQ=255)
                                                                                                                       tCGACACCTCGCGCTCGCGAGCGGGGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAggcggc                                                              <  4:146854/112‑1 (MQ=37)
                                                                                                                       tCGACACCTCGCGCTCGCGAGCGGGGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAggcggc                                                              >  3:146854/1‑112 (MQ=37)
                                                                                                                            cgcGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCTTCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCTCACCTtc                      <  6:226745/145‑1 (MQ=255)
                                                                                                                            cccGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCTCACCTtcgt                    <  1:356442/149‑1 (MQ=255)
                                                                                                                            cccGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCACTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCTCACCTtcga                    >  5:226745/1‑148 (MQ=255)
                                                                                                                             ccGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCTCACCTtcgtc                   <  2:308039/149‑1 (MQ=255)
                                                                                                                              cGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCTCACCTtcgtcg                  <  3:361640/149‑1 (MQ=255)
                                                                                                                              cGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCTCACCTtcgtcg                  >  3:175893/1‑149 (MQ=255)
                                                                                                                               gCTCGCGAGCGCGAGGTGTCGATGATCGACAACTTCGTCCTGCAAGGACGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCTCACCTTCGTCGc                 >  7:55773/1‑149 (MQ=255)
                                                                                                                                   gcgagtgcgagGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCTCACCTTCGTCGccccc             >  7:159036/1‑149 (MQ=255)
                                                                                                                                            agGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCTCACCTTCGTCGCCCCCGAGGTCAGc    <  8:159036/149‑1 (MQ=255)
                                                                                                                                              gtgtCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCTCACCTTCGTCGCCCCCGAGGTCAGCGt  <  1:272588/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GCAGGCCCTTGCCCAACAAGTAAAGCAGAATCCGCAGGTGGTAGCGGTTGCGCCGTTCACCCAGATGCAGGGCCTGCTGACCCATGATGGCAAGGTGCAGAAGGTGCTGCTCAATGGCATCGACCCCGCTCGCGAGCGCGAGGTGTCGATCATCGACAACTTCGTCCTGCAAGGCCGCCTGGACCAGTTGGCGCCGGGCGAGTGGGGCATCATGATCGGCGACAAGGCGGCGGCCAAGCTCGGCGTTGCCATCGGCGACAAGCTCACCTTCGTCGCCCCCGAGGTCAGCGT  >  NC_002947/2461232‑2461522

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: