Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947196,3060TC31.8% 56.5 / 30.3 44N604N (AAT→AAC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (19/11);  new base C (6/8);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACG  >  NC_002947/196157‑196454
                                                                                                                                                     |                                                                                                                                                    
aCGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAAt                                                                                                                                                      >  2:408911/1‑150 (MQ=255)
      ggTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGcc                                                                                                                                                <  1:408911/150‑1 (MQ=255)
           tttCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGAc                                                                                                                                           >  1:688106/1‑150 (MQ=255)
               gAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCg                                                                                                                                       >  2:132445/1‑150 (MQ=40)
                     cTGACCCCGAACCCGACGCCGGCTCAGACCACGAGCAACGACCCGACGTGTACTACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACGCCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACa                                                                                                                                 <  1:44155/150‑1 (MQ=11)
                      tGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAg                                                                                                                                >  2:226525/1‑150 (MQ=38)
                                cccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTc                                                                                                                      >  1:24109/1‑150 (MQ=255)
                                           cTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCaaa                                                                                                           <  2:688106/150‑1 (MQ=255)
                                           cTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCGGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCaaa                                                                                                           <  2:228450/150‑1 (MQ=255)
                                                aCCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACGATGACCAACGGCCAAGCCAtt                                                                                                      >  1:346848/1‑150 (MQ=255)
                                                  cACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGc                                                                                                                                                 <  1:248345/105‑1 (MQ=35)
                                                  cACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGc                                                                                                                                                 >  2:248345/1‑105 (MQ=38)
                                                  cACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTAc                                                                                                    >  1:155389/1‑150 (MQ=40)
                                                    cGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCg                                                                                                  >  1:621976/1‑150 (MQ=255)
                                                     gATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCGGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAAGGCCAAACCATTACAGt                                                                                                 >  2:557104/1‑150 (MQ=255)
                                                     gATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCGCTTCCACCTTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGt                                                                                                 >  2:587305/1‑150 (MQ=255)
                                                        cAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                              >  1:670759/1‑150 (MQ=40)
                                                          aCGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGGCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAg                                                                                            >  1:674625/1‑150 (MQ=255)
                                                              cTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCgg                                                                                        <  2:514404/150‑1 (MQ=40)
                                                                   tCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCaaga                                                                                   >  1:698502/1‑148 (MQ=38)
                                                                       caccacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCa                                                                               >  2:340222/1‑150 (MQ=21)
                                                                        accacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAg                                                                              >  1:223301/1‑150 (MQ=39)
                                                                            cTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCa                                                                          >  2:588452/1‑150 (MQ=14)
                                                                             caCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAg                                                                         <  2:527974/149‑1 (MQ=255)
                                                                             acacgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAg                                                                         >  1:316123/4‑150 (MQ=11)
                                                                                        tGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTcc                                                              >  1:522579/1‑150 (MQ=35)
                                                                                            ggCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCa                                                                               <  1:410073/129‑1 (MQ=255)
                                                                                            ggCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCa                                                                               >  2:410073/1‑129 (MQ=255)
                                                                                               gAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGAcccc                                                       <  1:26360/150‑1 (MQ=255)
                                                                                                     gTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGaa                                                 <  2:24109/150‑1 (MQ=35)
                                                                                                         gTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGTCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGAc                                             <  1:135633/150‑1 (MQ=11)
                                                                                                          tCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCAAGACCCCGACGAacgacg                                            >  2:395394/1‑150 (MQ=11)
                                                                                                              cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTa                                        >  2:356344/1‑150 (MQ=11)
                                                                                                                 aGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTacaa                                     <  1:85128/150‑1 (MQ=255)
                                                                                                                   gcggcgTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaaca                                   >  2:611263/1‑150 (MQ=35)
                                                                                                                         tactCACTTACACCGTGCCCCTAAGCAATCCTGCCCATACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCAGTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                             <  1:557104/147‑1 (MQ=11)
                                                                                                                         tGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCaa                                                                                                                   >  2:538492/1‑64 (MQ=37)
                                                                                                                         tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCAacaacacccgtt                             <  2:190880/150‑7 (MQ=11)
                                                                                                                         tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGCCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                             <  2:316123/150‑1 (MQ=11)
                                                                                                                          gATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc                            <  1:611263/150‑1 (MQ=35)
                                                                                                                              aCTTACACCGTGAGTCTCATCAATCCTGCCCAGACGCCGGTGACAGTGCCGCCGTGACTTCACAATTCTATAACCGGTAAAGCCGTCAACACCTAGGGCAGCGTCGATTTCCATACCCCGGCGAATGACGTCTACAACAACGGTTCGAcc                        <  1:340222/150‑2 (MQ=2)
                                                                                                                                ttACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCAt                                                                                                       <  2:371722/69‑1 (MQ=14)
                                                                                                                                ttACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCAt                                                                                                       >  1:371722/1‑69 (MQ=9)
                                                                                                                                      cGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAgcgt                >  1:378032/1‑147 (MQ=11)
                                                                                                                                                    aTCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACg  <  2:522579/150‑1 (MQ=39)
                                                                                                                                                     |                                                                                                                                                    
ACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACG  >  NC_002947/196157‑196454

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑