Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 747,900 G→C 12.7% intergenic (‑131/+21) abrB ← / ← ybgO regulator of aidB expression; inner membrane protein/putative fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913747,9000GC12.7% 173.7 / 14.6 55intergenic (‑131/+21)abrB/ybgOregulator of aidB expression; inner membrane protein/putative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (24/24);  new base C (3/4);  total (27/28)
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.49e-01

AAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACAT  >  NC_000913/747760‑748031
                                                                                                                                            |                                                                                                                                   
aaaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                              >  4:799064/1‑148 (MQ=255)
                                              aTACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAg                                                                                                                               <  2:1066682/101‑1 (MQ=255)
                                              aTACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAg                                                                                                                               >  1:1066682/1‑101 (MQ=255)
                                                   gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGa                                                                           <  3:799064/148‑1 (MQ=255)
                                                    cATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCaa                                                                                                                                 >  4:728137/1‑93 (MQ=255)
                                                    cATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCaa                                                                                                                                 <  3:728137/93‑1 (MQ=255)
                                                           aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATg                                                                                                                      <  2:509413/97‑1 (MQ=255)
                                                           aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATg                                                                                                                      >  1:509413/1‑97 (MQ=255)
                                                           aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                           <  1:186716/124‑1 (MQ=255)
                                                           aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                           >  2:186716/1‑124 (MQ=255)
                                                                         aaaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGc                                                                                                                            >  3:679965/1‑77 (MQ=37)
                                                                         aaaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGc                                                                                                                            <  4:679965/77‑1 (MQ=37)
                                                                         aaaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAggg                                                                                                                             >  1:735154/1‑76 (MQ=255)
                                                                         aaaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAggg                                                                                                                             <  2:735154/76‑1 (MQ=255)
                                                                                 aCAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgg                                                         <  3:810613/136‑1 (MQ=255)
                                                                                 aCAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgg                                                         >  4:810613/1‑136 (MQ=255)
                                                                                    gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATc                                                                                                               <  2:991130/79‑1 (MQ=255)
                                                                                    gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATc                                                                                                               >  1:991130/1‑79 (MQ=255)
                                                                                    gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                            >  1:306382/1‑82 (MQ=255)
                                                                                    gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                            <  2:306382/82‑1 (MQ=255)
                                                                                    gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                        <  1:932131/134‑1 (MQ=255)
                                                                                    gCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                        >  2:932131/1‑134 (MQ=255)
                                                                                      aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcaccac                                                                                   <  2:461465/105‑1 (MQ=255)
                                                                                      aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcaccac                                                                                   >  1:461465/1‑105 (MQ=255)
                                                                                        tCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                          >  4:109085/1‑112 (MQ=255)
                                                                                        tCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                          <  3:109085/112‑1 (MQ=255)
                                                                                         cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATgg                                                                                                                     >  1:862312/1‑68 (MQ=255)
                                                                                         cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATgg                                                                                                                     <  2:862312/68‑1 (MQ=255)
                                                                                           aGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATg                                                          <  4:436919/125‑1 (MQ=255)
                                                                                           aGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATg                                                          >  3:436919/1‑125 (MQ=255)
                                                                                             gAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGa                                                            <  2:912091/121‑1 (MQ=255)
                                                                                             gAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGa                                                            >  1:912091/1‑121 (MQ=255)
                                                                                                 aGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                 >  1:538301/1‑128 (MQ=255)
                                                                                                 aGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                 <  2:538301/128‑1 (MQ=255)
                                                                                                      aaaaTCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                         <  4:569128/99‑1 (MQ=255)
                                                                                                      aaaaTCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                         >  3:569128/1‑99 (MQ=255)
                                                                                                             gCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                         <  4:790658/92‑1 (MQ=255)
                                                                                                             gCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                         >  3:790658/1‑92 (MQ=255)
                                                                                                                            agaTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAg                                                                                         >  3:149913/1‑61 (MQ=255)
                                                                                                                            agaTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAg                                                                                         <  4:149913/61‑1 (MQ=255)
                                                                                                                                gCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGa                                                            >  3:836753/1‑86 (MQ=38)
                                                                                                                                gCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGa                                                            <  4:836753/86‑1 (MQ=38)
                                                                                                                                 cccTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                        <  4:212358/89‑1 (MQ=39)
                                                                                                                                 cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAACGCACCACTGTATCGAAGGGGCCAGCCTCGATggggg                                                      >  1:984104/1‑91 (MQ=255)
                                                                                                                                 cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggggg                                                      <  2:984104/91‑1 (MQ=255)
                                                                                                                                 cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                          <  2:640374/103‑1 (MQ=255)
                                                                                                                                 cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                          >  1:640374/1‑103 (MQ=255)
                                                                                                                                  ccTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtat                        <  4:342309/120‑1 (MQ=255)
                                                                                                                                  ccTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtat                        >  3:342309/1‑120 (MQ=255)
                                                                                                                                  ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                        >  3:63591/1‑88 (MQ=255)
                                                                                                                                  ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggg                                                        <  4:63591/88‑1 (MQ=255)
                                                                                                                                  ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaa                                     <  1:288351/107‑1 (MQ=255)
                                                                                                                                  ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTAGCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaa                                     >  2:288351/1‑107 (MQ=255)
                                                                                                                                   cTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACAt  <  3:943323/141‑1 (MQ=255)
                                                                                                                                   cTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACAt  >  4:943323/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                   
AAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACAT  >  NC_000913/747760‑748031

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