Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,782,5090AC23.2% 135.5 / 10.4 56V73G (GTA→GGA) ypjAadhesin‑like autotransporter YpjA
Reads supporting (aligned to +/- strand):  ref base A (30/13);  new base C (0/13);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 5.50e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAG  >  NC_000913/2782372‑2782641
                                                                                                                                         |                                                                                                                                    
tcgacccTGTGAAAATCCCCCCGTCATGAATTGACTGTGTGCCACCGTTAATGGTTGTATTATAAGCCTCCCCCACAAAATTGTTATGCCTTCCTATATCTTGATATCCACCAGATTCAAGAAGACTCCCATTAGATGcc                                                                                                                                    <  2:476633/135‑1 (MQ=255)
tCGTACCTGTGGAAATCCCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATAcc                                                                                                                                    <  1:684211/140‑1 (MQ=255)
   tACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCg                                                                                                                                 >  2:407025/1‑140 (MQ=255)
       tgtgGAAATCCCCCCGTCATGAATTGACTGTCTCCCCCCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATACCCACCAGATCCAATAAGCCTTCCATTAGATACCCCGCCAt                                                                                                                             <  2:264819/140‑1 (MQ=255)
       tgtgAAAATCCCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTAAATAAGACTTCCATTAGATCCCCCGCCAt                                                                                                                             <  2:2333/140‑1 (MQ=255)
          ggAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGaa                                                                                                                          >  2:729916/1‑140 (MQ=255)
            aaaTCCCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACa                                                                                                                        <  1:593101/140‑1 (MQ=255)
              aTACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACAtt                                                                                                                      >  2:649855/1‑140 (MQ=255)
                 cccccGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTctg                                                                                                                   <  2:80899/140‑1 (MQ=255)
                      gTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGc                                                                                                              >  1:387446/1‑140 (MQ=255)
                       tCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGcc                                                                                                             <  1:41756/140‑1 (MQ=255)
                       ataggAATTGACTGTCTGCCACGGTTAATGGTTATATTTTTAGCCTCCCCCACAAAATTGTTATGACTTCCTATATCTGGATACCCACCAAATTCAATAAGACTTCCTTTAATTCCCCCGCCATGAACATTCTGCTGGcc                                                                                                             <  2:692884/136‑1 (MQ=255)
                         aTGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCAt                                                                                                           <  2:219510/140‑1 (MQ=255)
                         aTGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTAAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCAt                                                                                                           <  1:609105/140‑1 (MQ=255)
                           gAATTGACTGTCTGCCACCGTTAATGTTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCCCCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATgg                                                                                                         <  2:542910/140‑1 (MQ=255)
                              ttGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTg                                                                                                      >  1:663337/1‑140 (MQ=255)
                               tGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCCCCAGATTCAATAAGTCTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGa                                                                                                     <  2:34766/140‑1 (MQ=255)
                                 aCTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGata                                                                                                   >  1:74427/1‑140 (MQ=255)
                                 aCTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGata                                                                                                   >  1:42684/1‑140 (MQ=255)
                                     tttGCCGCGGTTAATGGTTGTATTATTACCCTCCCACACAAAATTTTTATGCCTTCTTATATCTTGATATCCACCAGTTTCAATAAGCCTCCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAatat                                                                                               <  2:350080/138‑1 (MQ=255)
                                     tCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAatat                                                                                               <  2:663337/140‑1 (MQ=255)
                                       tGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGa                                                                                             <  2:387446/140‑1 (MQ=255)
                                            ccGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATTTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGtt                                                                                        <  1:614817/140‑1 (MQ=255)
                                             cGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTa                                                                                       <  1:221674/140‑1 (MQ=255)
                                               ttAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTAtt                                                                                     >  1:390230/1‑140 (MQ=255)
                                                       tGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTTGCCATGGTTGATAATATGAGTGTTATTTGTTGTAc                                                                             >  1:59517/1‑140 (MQ=255)
                                                       tGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTAc                                                                             >  2:412729/1‑140 (MQ=255)
                                                       tGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTCCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTAc                                                                             >  1:508439/1‑140 (MQ=255)
                                                         tattatTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCt                                                                           >  1:61815/1‑140 (MQ=255)
                                                                   cTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCa                                                                 >  1:72828/1‑140 (MQ=255)
                                                                    tGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCAc                                                                >  1:534295/1‑140 (MQ=255)
                                                                      ccccACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACtt                                                              >  1:604236/1‑140 (MQ=255)
                                                                      ccccACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACtt                                                              >  2:228273/1‑140 (MQ=255)
                                                                      ccccACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACtt                                                              >  2:118488/1‑140 (MQ=255)
                                                                      ccccACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACtt                                                              >  1:90915/1‑140 (MQ=255)
                                                                      ccccACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACtt                                                              >  1:402102/1‑140 (MQ=255)
                                                                           caAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTgg                                                         <  1:119019/140‑1 (MQ=255)
                                                                           aaaaaTTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTgg                                                         <  1:189913/139‑1 (MQ=255)
                                                                                ttGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGcc                                                    >  2:358628/1‑140 (MQ=255)
                                                                                ttGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGcc                                                    <  2:42684/140‑1 (MQ=255)
                                                                                ttGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGcc                                                    >  2:95862/1‑140 (MQ=255)
                                                                                ttGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGcc                                                    >  1:219822/1‑140 (MQ=255)
                                                                                   ttATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATc                                                 <  1:649855/140‑1 (MQ=255)
                                                                                     aTGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTa                                               <  2:267192/140‑1 (MQ=255)
                                                                                     aTGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTa                                               <  2:74427/140‑1 (MQ=255)
                                                                                        aCTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAg                                            >  1:292616/1‑140 (MQ=255)
                                                                                        aCTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAg                                            >  1:471238/1‑140 (MQ=255)
                                                                                           tCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGtct                                         <  2:215555/140‑1 (MQ=255)
                                                                                             cTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCa                                       >  1:406653/1‑140 (MQ=255)
                                                                                             cTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCa                                       >  1:527704/1‑140 (MQ=255)
                                                                                               atatCTTGATATCCACCAGATTCAATAAGATTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATc                                     <  2:224209/140‑1 (MQ=255)
                                                                                                  tCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGtt                                  <  1:412729/140‑1 (MQ=255)
                                                                                                     taatacCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTAc                               <  2:480586/138‑1 (MQ=255)
                                                                                                      gATATCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTAcc                              <  1:228273/140‑1 (MQ=255)
                                                                                                       atatCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACca                             <  1:614154/140‑1 (MQ=255)
                                                                                                       atatCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACca                             >  1:67150/1‑140 (MQ=255)
                                                                                                       atatCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACca                             <  2:90915/140‑1 (MQ=255)
                                                                                                        tatCCACCAGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACcac                            <  1:158411/140‑1 (MQ=255)
                                                                                                            caccaGATTCAATAAGACTTCCATTAGATCCCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCa                        <  2:61815/140‑1 (MQ=255)
                                                                                                             accaGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCaa                       >  2:415291/1‑140 (MQ=255)
                                                                                                                        tAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATtagta            >  1:681026/1‑140 (MQ=255)
                                                                                                                              ttCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAg      >  2:151872/1‑140 (MQ=255)
                                                                                                                                  aTTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAGCATTAGTAGTGAAGGCAg  <  2:536467/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
TCGTACCTGTGGAAATACCCCCGTCATGAATTGACTGTCTGCCACCGTTAATGGTTGTATTATTAGCCTGCCCCACAAAATTGTTATGACTTCCTATATCTTGATATCCACCAGATTCAATAAGACTTCCATTAGATACCCCGCCATGAACATTCTGCTGGCCATGGTTGATAATATGAGTGTTATTTGTTGTACCTCGTTCATCCACTTTTTGGTTGCCATCTACAGTCTCATCGTTTACCACACCAATAACATTAGTAGTGAAGGCAG  >  NC_000913/2782372‑2782641

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: