Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947197,5720GA22.7% 74.9 / 19.7 44K1026K (AAG→AAA) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (18/16);  new base A (4/6);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGAT  >  NC_002947/197426‑197712
                                                                                                                                                  |                                                                                                                                            
gATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGa                                                                                                                                             >  7:33807/1‑148 (MQ=37)
   gCCACCACCGCGACCCTGACGGCGAGCACGTCGGTCACCGAAGTCGCCGTGATCACCTACACCGTGATCATGAGCAATACTGCCAAGACGACGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAc                                                                                                                                            >  7:77931/1‑146 (MQ=11)
      accaccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCCGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAgg                                                                                                                                       <  5:176837/148‑1 (MQ=37)
          ccGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAg                                                                                                                                   >  4:48655/1‑148 (MQ=37)
              gACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTc                                                                                                                               <  4:178160/148‑1 (MQ=37)
                   tGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTTATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAttt                                                                                                                          <  5:166918/148‑1 (MQ=37)
                      cGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCa                                                                                                                       <  3:1809/148‑1 (MQ=17)
                             cccGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCg                                                                                                                <  8:2134/148‑1 (MQ=37)
                                gTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGc                                                                                                              >  4:69194/1‑147 (MQ=37)
                                    gTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                          <  7:163539/147‑1 (MQ=37)
                                         cGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGt                                                                                                     >  3:64294/1‑147 (MQ=37)
                                         cGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGt                                                                                                     >  5:61167/1‑147 (MQ=37)
                                          gAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCt                                                                                                   <  3:48655/148‑1 (MQ=37)
                                          gAAGGCGGCGTGATCAACTAAACCGTGATCCTGAGCAATCCTGCCAAGACGCAGGTAACCGTGACCCTGTCCAACGGCCAAACCAATAACGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAGTTCCAGACCCCGGCGAATGACGTCt                                                                                                   >  2:155006/1‑148 (MQ=11)
                                                   gTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACg                                                                                          <  8:112225/148‑1 (MQ=37)
                                                           cTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGACTTc                                                                                                                         >  3:164291/1‑109 (MQ=255)
                                                           cTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCACGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAAAAACGGTTCGAcc                                                                                  >  6:113777/1‑147 (MQ=255)
                                                           cTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                  <  8:33807/148‑1 (MQ=38)
                                                            tACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTTACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTcc                                                                                                                        <  4:164291/109‑1 (MQ=11)
                                                                             cAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                                                                                       >  4:44964/1‑125 (MQ=35)
                                                                              aaTCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc                                                                                      <  3:44964/125‑1 (MQ=35)
                                                                                tCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAAc                                                             >  7:158678/1‑148 (MQ=17)
                                                                                  cTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGACGCCGGCAAGACCCAGGGCAGCGTCGCTTTCCAGACCCCGGCGAATGACGTCTACAACACCGGTTCGACCGTCAGCGTCCCCATCGAGAACGc                                                           >  8:67666/1‑148 (MQ=11)
                                                                                   tGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACAGTCAGCGTCACCATCGAGAACGcc                                                          >  8:89089/1‑148 (MQ=21)
                                                                                     cccTGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACGGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCcac                                                        <  1:101750/148‑1 (MQ=255)
                                                                                       cAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGCTTCGACTGTCAGCGTCCCCATCGAGAACGCCACAg                                                      >  7:47014/1‑148 (MQ=255)
                                                                                        aGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAgg                                                     >  8:142436/1‑148 (MQ=17)
                                                                                         gACGCCTGTATCCGTGACCCTGTCCAAAGGACAAAACATTACGGTTGAAGCTTGCAAGACCAAGTTCAGCGTCGATTTGCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTAACCATCGAGAACGCCACAGGc                                                    <  8:77931/148‑1 (MQ=11)
                                                                                          aCGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                   <  6:3631/148‑1 (MQ=255)
                                                                                            gCCGGTGACCGTGACCCTGTCCAACGGCCAACCGATTACCGTTGAAGCCGGCCAACCCCAGGGCAGCGTCGACTTCCAGACACCGGCGAATGACGTCGACAACCACGGGTCGACTGTCAGCGTCACCCGCGAGAACGCCCCAGGCgg                                                  >  7:173307/1‑147 (MQ=11)
                                                                                              cGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTaa                                               <  6:61167/148‑1 (MQ=37)
                                                                                                    ccGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                         <  4:49947/148‑1 (MQ=16)
                                                                                                      gTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAg                                        >  1:107331/1‑147 (MQ=32)
                                                                                                          cccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACt                                   >  4:101246/1‑148 (MQ=32)
                                                                                                                ccAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAccc                              <  2:56317/147‑1 (MQ=32)
                                                                                                                   aCGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGaa                          <  1:165293/148‑1 (MQ=12)
                                                                                                                        cAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGa                     >  2:73802/1‑148 (MQ=34)
                                                                                                                         aaaCGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGAc                    >  6:66542/1‑148 (MQ=255)
                                                                                                                           aCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGc                  <  3:47271/148‑1 (MQ=35)
                                                                                                                             cATGACCGTGGCCGCCGGCAAGACCCACGGCAGCGGCGATTTCCAGCCCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGGCACCATCGAGCCCGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCg                <  5:21214/148‑1 (MQ=11)
                                                                                                                                    gtggAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGa         <  4:149514/145‑1 (MQ=255)
                                                                                                                                     ttGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGa         >  6:13750/1‑147 (MQ=12)
                                                                                                                                        aaGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAccc                       >  8:21346/1‑130 (MQ=25)
                                                                                                                                           ccGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGAt  <  3:117134/148‑1 (MQ=12)
                                                                                                                                                  |                                                                                                                                            
GATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGAT  >  NC_002947/197426‑197712

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