Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA46.2% 19.7 / 13.9 27intergenic (‑330/+60)ybgD/gltAputative fimbrial‑like adhesin protein/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (10/4);  new base A (0/12);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 2.01e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATT‑GAAAGTATTGGGTTGCTG‑ATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTC  >  NC_000913/752982‑753265
                                                                                                                                                 |                                                                                                                                            
cATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATT‑GAAAGTATTGGGTTGCTG‑ATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATc                                                                                                                                       <  1:150353/151‑1 (MQ=255)
                        aTGAAGTTTTTAATT‑AAAAGTATTGGTTTGCTT‑ATAATTTGTGCTTTTCTTTTCTTTTTAAATATCTATATAGGTTTTTTAATGTATTTTATTTTTACAATTTTTGGTGTTGAGGCAAAAAAAAATCAACCCGCCATATGAACGGCGGGt                                                                                                                <  1:325475/150‑1 (MQ=255)
                           aattttttATTT‑GAAAGTATTGGG‑TGCGGCATAATTTAGGTTGTTATATCTTTTAAAAATATGGTTATCGATTTGTTAATGATTTTTATTTTTAAATTTTTTTGTGTAAAGCAATTTAAAAACCACCCCGCCATATGAACGGCGGGtt                                                                                                               <  2:542800/145‑1 (MQ=255)
                           aattttttAATT‑AAAATTATTGGGTTGCTA‑ATATTTTGAGTTTTTTTTTTTTTTTAAAATATTTATATGGTTCTGTAAATGGTTTTTTTTTTTACAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTa                                                                                                              <  2:155257/146‑1 (MQ=255)
                               ttttAATTGGAAAGTATTGGGTTGGTG‑ATATTTTGAGTTGTTCTATTTTTTTTAAATATCTATAGAGGTGGGTAAAGGGATTTTATTTTTACATTTTTTTGTGTTTAGGATTATAAAAATAAACCCGCCATATGAACGGCGGGTTaaaa                                                                                                           <  1:171882/149‑1 (MQ=255)
                                               ttGGGTTGCTG‑ATATTTTGAGCTTTTTTATTTTTTTTAAATATTTATATAGGTTTGTTAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCa                                                                                           <  2:544867/149‑1 (MQ=255)
                                                    ttGCTG‑ATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAAcc                                                                                    >  2:301480/1‑151 (MQ=255)
                                                                ttGAGCTTTTTTATTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCtt                                                                          <  2:431149/150‑1 (MQ=255)
                                                                     cTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGC‑TGATATCg                                                                   >  2:58030/1‑151 (MQ=255)
                                                                       gTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                  >  1:346060/1‑151 (MQ=255)
                                                                              ttttttaAAATATCTATATAGTTCTGTAAATGGATTTTATTTTTAAAATTTTTTGTGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTaa                                                            <  2:89453/148‑1 (MQ=255)
                                                                              tCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTaa                                                            >  2:285544/1‑150 (MQ=255)
                                                                                    tAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcgcg                                                     >  2:287676/1‑151 (MQ=255)
                                                                                          tttatataGTTTTGTTAATGGTTTTTATTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                 <  2:410032/147‑1 (MQ=255)
                                                                                                ttggttttttAATGGTTTTTATTTTAACATTTTTTTGTGTTTAGTCATAAAAAAAGCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtt                                                   <  1:167325/133‑1 (MQ=255)
                                                                                                  ggTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCt                                        >  2:42292/1‑150 (MQ=255)
                                                                                                        ttAATGGATTTTATTTTTACATTTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta                                     <  2:106652/147‑1 (MQ=255)
                                                                                                        ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta                                     <  2:358/147‑1 (MQ=255)
                                                                                                             ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTTAAGTCGCGTTTTTCATATCCTGTATACAGCt                              >  1:126874/1‑149 (MQ=255)
                                                                                                              tattttattttACCATTTTTTTATATTAGGAAATACAAAAATCACTCCGCCATATGAACGACGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGa                            <  1:385684/149‑1 (MQ=255)
                                                                                                              gtttttttttttACATTTTTTTGTGTTTGTGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCt                              <  1:133664/146‑1 (MQ=255)
                                                                                                                  ttattttTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgcg                        <  1:301396/150‑1 (MQ=255)
                                                                                                                     tttttAAATTTTTTTGTTTTAAGGAAAAAAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGAc                     <  1:449415/150‑1 (MQ=255)
                                                                                                                      ttttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAAAACTTAGCAATCAACCATTAAAGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGa                      >  1:371341/1‑148 (MQ=255)
                                                                                                                                    gtgtTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAACAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAAAGCTGACGCGGACGGGCAATCTTCATAcc     >  2:29410/1‑151 (MQ=255)
                                                                                                                                       ttaaGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGt   <  2:166808/147‑1 (MQ=255)
                                                                                                                                        ttAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTc  >  1:335520/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
CATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATT‑GAAAGTATTGGGTTGCTG‑ATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTC  >  NC_000913/752982‑753265

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑