Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947195,5590CA41.2% 20.1 / 29.3 34T355T (ACC→ACA) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (7/13);  new base A (7/7);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGA  >  NC_002947/195424‑195695
                                                                                                                                       |                                                                                                                                        
gaccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGGCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTc                                                                                                                              <  6:88422/148‑1 (MQ=25)
  ccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                                                                                                            <  3:52873/148‑1 (MQ=34)
      gaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGCCCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAACTTCGAGCat                                                                                                                        >  7:69893/1‑147 (MQ=12)
         ccTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTg                                                                                                                     >  4:101488/1‑148 (MQ=39)
          cTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGa                                                                                                                    <  6:10709/148‑1 (MQ=255)
                                 tGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCg                                                                                                                             >  5:60445/1‑116 (MQ=18)
                                     gAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACACCGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCa                                                                                         >  7:78711/1‑148 (MQ=37)
                                           ggtaaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGAc                                                                                           >  4:4524/4‑140 (MQ=37)
                                            gtaaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGAcc                                                                                          <  3:4524/138‑1 (MQ=17)
                                                aaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCTGCGGTAATTTCGAGCAACTGACCCCGAATCAGACCCCGGCTCAGACCACGATCAACGAc                                                                              <  6:33834/148‑1 (MQ=25)
                                                      gggCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTc                                                                        <  6:143556/148‑1 (MQ=34)
                                                      gggCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTc                                                                        <  1:8398/148‑1 (MQ=34)
                                                                           cccGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTg                                                   <  2:24317/148‑1 (MQ=39)
                                                                                cgaaTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGc                                              <  1:75016/146‑1 (MQ=2)
                                                                                          cTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCg                                    <  3:101488/148‑1 (MQ=37)
                                                                                                caacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCAc                               >  6:114746/1‑147 (MQ=37)
                                                                                                           ccGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                   <  5:114746/148‑1 (MQ=37)
                                                                                                           ccGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                   <  3:10828/148‑1 (MQ=17)
                                                                                                            cGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTg                  <  8:78711/148‑1 (MQ=38)
                                                                                                               catcGTCACCATCTACAACGCCACCGGCGGTAATTTCGAGCAACTGTCACAGAATCCCACCCCCGCTCCGACCTCGCTGATCGACTCTTTCGACACCACAACCGCGACCCTGACGGCGAGCCCTTCGGTCACCGAAGGCGGCGTGATc               <  8:24627/145‑1 (MQ=11)
                                                                                                                 gCGTCACCATCGAGAACGCCACCGGCGGTAATTTCg                                                                                                                             >  4:84562/1‑36 (MQ=37)
                                                                                                                  cGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                                                                                                            <  3:84562/36‑1 (MQ=33)
                                                                                                                  cGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGgcggcg                    >  6:102234/1‑140 (MQ=17)
                                                                                                                   gTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                   <  5:102234/140‑1 (MQ=17)
                                                                                                                    tCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACGACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTa          >  1:91300/1‑148 (MQ=12)
                                                                                                                      aCCATCGAGAACCCCACAGGCGGTAATTTCGAGCAACTGCCCCCGACTCCGACCCCGCCACAGACCACGATCAACGACTCGGTCGACCCGACTACCGTGCCCCTCACGGCGAGCCCGTCCGGCCCCGAAGGCGGCGAGAACACTTaca        >  8:110062/1‑148 (MQ=11)
                                                                                                                        cATCGAGAACGCCACCGGCGGTAATTTCGAGCATCTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACAcc      >  7:75350/1‑148 (MQ=21)
                                                                                                                        cATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACAACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGa                            >  3:23280/1‑126 (MQ=17)
                                                                                                                         aTCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGaa                           <  4:23280/126‑1 (MQ=17)
                                                                                                                            gagaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGCTCAACGACTCGGTCGACACCACTACCGCGACCCTGCCGGCGAGCCCGTCGGTCACCGAAGgcggcg                    >  5:89372/1‑130 (MQ=255)
                                                                                                                             agaACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGa  <  2:169410/147‑1 (MQ=39)
                                                                                                                             agaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                   <  6:89372/130‑1 (MQ=17)
                                                                                                                               aaCGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTac         >  7:69370/1‑138 (MQ=21)
                                                                                                                                aCGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTaca        <  8:69370/138‑1 (MQ=21)
                                                                                                                                       cggcggTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCg                                             >  1:60796/1‑94 (MQ=38)
                                                                                                                                       |                                                                                                                                        
GACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGA  >  NC_002947/195424‑195695

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