Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1140GT25.0% 128.3 / 11.8 56intergenic (‑319/+71)ybgD/gltAputative fimbrial protein YbgD/citrate synthase
Reads supporting (aligned to +/- strand):  ref base G (26/16);  new base T (0/14);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 2.67e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGC  >  NC_000913/752976‑753250
                                                                                                                                          |                                                                                                                                        
cATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTgtgt                                                                                                                                         >  2:193453/1‑140 (MQ=255)
 aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGtt                                                                                                                                        <  1:1015907/140‑1 (MQ=255)
     tCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAgg                                                                                                                                    <  1:193453/140‑1 (MQ=255)
               ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                                          >  2:754153/1‑140 (MQ=255)
               ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                                          >  1:612337/1‑140 (MQ=255)
                taAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAAt                                                                                                                         >  2:859717/1‑140 (MQ=255)
                     gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATATAAAAATCAAcc                                                                                                                    <  2:89287/140‑1 (MQ=255)
                     gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                    >  1:18570/1‑140 (MQ=255)
                      cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                                   >  1:290157/1‑140 (MQ=255)
                        aCTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTTTTCTATTTTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGc                                                                                                                 <  1:976722/140‑1 (MQ=255)
                                       ttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCggg                                                                                                  >  2:87115/1‑140 (MQ=255)
                                            taaaaGTATTGGGTTGCTGATAATTTGAGTTTTTTTTTTTTTTTAAAATATCTATATAGGTCTTTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTaaa                                                                                             <  2:820524/138‑1 (MQ=255)
                                            tGAAAGTATTGGTTTGTTGATAATTTGAGCTGTTCTATTTTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTaaa                                                                                             <  1:338472/140‑1 (MQ=255)
                                                   tttGGGTTGCTGATATTTTGAGTTTTTTTATTTTTTTTAAATATCTATATAGGTTTGTTAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGAATAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAc                                                                                     <  2:776323/139‑1 (MQ=255)
                                                         ttGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTa                                                                               >  1:776323/1‑140 (MQ=255)
                                                            cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                             >  1:479439/1‑140 (MQ=255)
                                                            cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATGTTTACAATTTTTTGTGTTTAGGTATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                             >  2:910518/1‑140 (MQ=255)
                                                             tGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCa                                                                            <  1:87115/140‑1 (MQ=255)
                                                                 aaTTTGAGTTTTTCTTTTTTTTTTAAATATCTATATAGGTCTTTTAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGGAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCa                                                                        <  2:246025/140‑1 (MQ=255)
                                                                 aaTTTGAGTTGTTTTTTTCTTTTAAAATATCTATATAGGTTTTTTAATGGTTTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCa                                                                        <  1:507823/140‑1 (MQ=255)
                                                                 aaTTTGAGCTTTTCTTTTCTTTTTAAATATCTATATAGGTTTGTTAATGGTTTTTATTTTTACATTTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCa                                                                        <  2:972396/140‑1 (MQ=255)
                                                                      gAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                   <  1:859717/140‑1 (MQ=255)
                                                                       ggCTGTTTTATTTTTTTAAAATTTCTTTAGAGGTCTGTAAATGGATTTTTTTTTTACAATTTTTTTTGTTTAGGTATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAtt                                                                  <  1:818434/139‑1 (MQ=255)
                                                                           tttatATTCTTTTAAAATATTTATATAGGTCGGTTAATGGATTTTATTTTTATATTTTTTTGTGTTTAGCCAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                              <  2:1155292/136‑1 (MQ=255)
                                                                               tATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCAGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                   <  2:774658/131‑1 (MQ=255)
                                                                               tATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCAGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                   >  1:774658/1‑131 (MQ=255)
                                                                                tttttttatAAATATCTATATAGGTTTGTTAATGTATTTTTTTTTTACATTTTTTTTTTTTTAGGCATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGa                                                         <  2:1032616/136‑1 (MQ=255)
                                                                                      tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAAACAACCATTAAAGCTTGATATCGc                                                   >  2:13026/1‑140 (MQ=255)
                                                                                      tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                                     <  1:888483/138‑1 (MQ=255)
                                                                                      tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                                     >  2:888483/1‑138 (MQ=255)
                                                                                      tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                   >  2:1024308/1‑140 (MQ=255)
                                                                                      tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAACCATCAACCATTAAAGCTTCataaacc                                                   >  1:39174/1‑136 (MQ=255)
                                                                                      tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGCTAAAATATTTACAACCTAGCACTCAACCATTAACGCTTGATATAGc                                                   >  1:1155292/1‑140 (MQ=255)
                                                                                                 atataGTTTTTTTAATGTTTTTTTTTTTTAAATTTTTTTGTTTTTAGAAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcg                                        <  2:1169386/140‑1 (MQ=255)
                                                                                                   ataGTTTTGTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGAGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgcg                                      <  2:612337/140‑1 (MQ=255)
                                                                                                     aGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtt                                    >  2:514069/1‑140 (MQ=255)
                                                                                                      ggTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTAAAGTCGCGttt                                   >  1:145488/1‑140 (MQ=255)
                                                                                                       gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                  >  1:848702/1‑140 (MQ=255)
                                                                                                       gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                  >  2:634094/1‑140 (MQ=255)
                                                                                                       gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                  >  1:966479/1‑140 (MQ=255)
                                                                                                        tCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttttt                                 >  1:97401/1‑140 (MQ=255)
                                                                                                          tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttttta                               >  1:1169386/1‑138 (MQ=255)
                                                                                                          tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                               >  2:1132054/1‑140 (MQ=255)
                                                                                                          tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAATCGCGttttttt                               >  1:314992/1‑138 (MQ=255)
                                                                                                           tttAATGTTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                              <  2:1078633/139‑1 (MQ=255)
                                                                                                           tttAATGTATTTTATTTTTACAATTTTTTTTTTTTAGACATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                              <  2:931801/139‑1 (MQ=255)
                                                                                                           tttAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGAAAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                              <  2:82515/139‑1 (MQ=255)
                                                                                                           tttAATGGTTTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                              <  1:456910/139‑1 (MQ=255)
                                                                                                           tttAATGGTTTTTATTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                              <  1:857672/139‑1 (MQ=255)
                                                                                                            ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                             >  2:486850/1‑140 (MQ=255)
                                                                                                             tAATGGTTTTTATTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCatat                            <  1:484804/140‑1 (MQ=255)
                                                                                                              aatgttttttttttttACATTTTTTTGTGTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                           <  2:120845/134‑1 (MQ=255)
                                                                                                              aatgttttttttttttACAATTTTTTTTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                           <  2:654164/134‑1 (MQ=255)
                                                                                                              aaTGGTTTTTATTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                           <  2:97401/140‑1 (MQ=255)
                                                                                                              aaTGGATTTTATTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                           <  1:207043/140‑1 (MQ=255)
                                                                                                                   tttttttttttAAATTTTTTTGTGTTTAGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta                      <  2:18570/139‑1 (MQ=255)
                                                                                                                   tttttattttTCAAATTTTTGGTGTTAAGCCATAAAAAAATAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta                      <  2:39174/139‑1 (MQ=255)
                                                                                                                     tttattttTAAAATTTTTTGTTTTTAGGAAAAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtata                    <  2:314992/140‑1 (MQ=255)
                                                                                                                      ttattttTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAc                   >  2:721170/1‑140 (MQ=255)
                                                                                                                                aaTTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgcg         <  1:514069/140‑1 (MQ=255)
                                                                                                                                       tgttTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc  <  2:679543/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                        
CATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGC  >  NC_000913/752976‑753250

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