Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,679,7250GA75.0% 56.6 / 19.5 36W70* (TGG→TAG) cbpDpossible secretory antigen
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base A (15/12);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGAACAAACGCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTT  >  CP000730/2679595‑2679861
                                                                                                                                  |                                                                                                                                        
gTCGAACAAACGCATCGATCGCATGCCGCAGTTAATTATTATAGAAAAAACAAATGTACATGGTAGGCATTTAAACGTCGCGCACAAGTCGCTAAACCTGTTCCTAATAGATGGAGCAATGCTAAAAATTGGGATTACAAt                                                                                                                                >  2:357771/1‑141 (MQ=255)
 tCGAACAAACGCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATg                                                                                                                               >  2:378429/1‑141 (MQ=255)
 tCGAACAAACGCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATg                                                                                                                               >  1:351768/1‑141 (MQ=255)
  cGAACAAACGCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGc                                                                                                                              >  2:150673/1‑141 (MQ=255)
     aCAAACGCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACg                                                                                                                           >  1:12272/1‑141 (MQ=255)
           gCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATc                                                                                                                     <  1:34331/141‑1 (MQ=255)
             aTCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCaa                                                                                                                   <  1:228042/141‑1 (MQ=255)
                  tCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAatat                                                                                                              >  1:46834/1‑141 (MQ=255)
                  tCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAatat                                                                                                              <  1:65769/141‑1 (MQ=255)
                       tGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGAc                                                                                                         >  1:313831/1‑141 (MQ=255)
                         ccGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTg                                                                                                       >  2:327467/1‑141 (MQ=255)
                           gCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGt                                                                                                     >  1:233197/1‑141 (MQ=255)
                                tAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTac                                                                                                <  2:49339/141‑1 (MQ=255)
                                     attatAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAACCCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCaa                                                                                           <  2:78239/141‑1 (MQ=255)
                                     attatAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCaa                                                                                           <  2:149362/141‑1 (MQ=255)
                                                 aCCAATGTACATGGTTGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTg                                                                               >  2:296079/1‑141 (MQ=255)
                                                   cAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTc                                                                             >  2:176711/1‑141 (MQ=255)
                                                   cAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGATAAACCTGTTAATAATAGATAGGACAATGCTAAAAATAAGTATTACAATGCACGTAAATAAAAATATGCAACTGGTCGTACACAAAGaaaaatagatata                                                                             >  2:274987/1‑132 (MQ=255)
                                                          cATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAAt                                                                      <  2:289662/141‑1 (MQ=255)
                                                            tggtggGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCa                                                                    <  1:326312/141‑1 (MQ=255)
                                                            tggtggGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCa                                                                    >  1:91544/1‑141 (MQ=255)
                                                            tggtggGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACAAAGAAAATTTGCTGTCATGCAAtta                                                                    >  2:336071/1‑139 (MQ=255)
                                                                       tAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGata                                                         >  2:173627/1‑141 (MQ=255)
                                                                                cgcACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACa                                                <  2:408080/141‑1 (MQ=255)
                                                                                  cacaAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATg                                              <  2:46834/141‑1 (MQ=255)
                                                                                      aGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGc                                          >  1:18054/1‑141 (MQ=255)
                                                                                                     ttCTAATAGATGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGt                           <  2:388037/141‑1 (MQ=255)
                                                                                                       cTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGtat                         <  2:45858/141‑1 (MQ=255)
                                                                                                          aTAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAaaaaaaaaaaa                      <  2:313831/141‑8 (MQ=255)
                                                                                                               tGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAAc                 <  2:219350/141‑1 (MQ=255)
                                                                                                               tGGGGCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAAc                 >  2:320939/1‑141 (MQ=255)
                                                                                                                  ggCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGgt              >  2:258543/1‑141 (MQ=255)
                                                                                                                  ggCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCCTGAAATCAACTGCAGGATATTATGGCCATGTCGCCATTGTTGAAAAAGTATATAAAAACGgt              >  1:334772/1‑141 (MQ=255)
                                                                                                                  ggCAATGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGgt              <  2:273375/141‑1 (MQ=255)
                                                                                                                       tGCTAAAAATTAGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTAt         <  1:41046/141‑1 (MQ=255)
                                                                                                                              aaTTAGTATTACAATGCACGTAAATCAAAAGATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGtt  >  2:86124/1‑141 (MQ=255)
                                                                                                                                  |                                                                                                                                        
GTCGAACAAACGCATCAATCGCATGCCGCAGTTAATTATTATAGTAAAAACCAATGTACATGGTGGGCATTTAAACGTCGCGCACAAGTCGGTAAACCTGTTTCTAATAGATGGGGCAATGCTAAAAATTGGTATTACAATGCACGTAAATCAAAATATGCGACTGGTCGTACACCAAGAAAATTTGCTGTCATGCAATCAACTGCAGGATATTATGGACATGTCGCAGTTGTTGAACAAGTATATAAAAACGGTAGTATTAAAGTT  >  CP000730/2679595‑2679861

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