Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,369,520 T→C V43A (GTT→GCT)  pspD → peripheral inner membrane phage‑shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,369,5200TC95.1% 138.2 / ‑6.3 41V43A (GTT→GCT) pspDperipheral inner membrane phage‑shock protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (9/30);  minor base A (2/0);  total (11/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGATGGCAATTT  >  NC_000913/1369393‑1369666
                                                                                                                                |                                                                                                                                                    
aTGAATAATCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTAGGGGATTAAAACAGACGATCGcccccc                                                                                                                                          >  1:253672/1‑136 (MQ=255)
          gCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGCGGGGGGTAAAAACAACTGCTCGCCGACCCCCGAAAAtttt                                                                                                                              >  2:661002/1‑142 (MQ=255)
             ggCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCtgtg                                                                                                                    <  2:160152/149‑1 (MQ=255)
                      ccGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAAcc                                                                                                       <  1:473939/151‑1 (MQ=255)
                        gggCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGc                                                                                                     <  1:202148/151‑1 (MQ=255)
                        gggCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGc                                                                                                     <  2:336181/151‑1 (MQ=255)
                            aaaaGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGCGGCGGGCTTGGTGCTTAAATTACCTGCTCCCCGCCCCCGGAAAATTTTGCCGGCTGGGGCGCTGGaaactccgtt                                                                                                 >  2:63697/1‑143 (MQ=255)
                                 gTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTc                                                                                            <  2:404152/151‑1 (MQ=255)
                                      gCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGgct                                                                                       <  1:506760/151‑1 (MQ=255)
                                          ggTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGgctgct                                                                                    <  1:553833/150‑1 (MQ=255)
                                               cAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTGGGGGGTAAAAATAACCGCTTGCCCACCCCTGGAAAAATTGGTGGGTGGGGCGGCTGGACCCCCGTTTAATCGGGGTGGtaataa                                                                               >  2:42915/1‑150 (MQ=255)
                                               cAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCtaataa                                                                               <  2:14478/150‑1 (MQ=255)
                                                        aGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGc                                                                    <  1:111047/151‑1 (MQ=255)
                                                         ggCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGc                                                                    <  1:209132/151‑1 (MQ=255)
                                                           cAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTGGGCGTTTAAATTAGATTCCAGCCGACCACTGAAAAAATTTCTTGCTTTGGGGGCGGAACCCCTGTTTAGTCGGGCCGCTaaaaaaactgaac                                                                       >  2:781273/1‑139 (MQ=255)
                                                            aaGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGtt                                                                 <  1:411905/151‑1 (MQ=255)
                                                              gCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTtat                                                               <  2:309644/151‑1 (MQ=255)
                                                                      ttcttACCGCACTGCGCTATGGCCCGGCGGGGGTGGCGGGGTTGGGGGTAAAAAAAGACGGCTCCCCGCCCCCGGAAAAGGTTGTGGGGGGGGGGGCGGGAACCGCGGTTAATCGGGGCGGCTATAAAATGGGCCCGCGGTTTaaaaaggg                                                       >  1:354704/1‑147 (MQ=255)
                                                                              gCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGGTAAAAACAGCTGCTCGCCCACCGCCGGAAAAGGTGGTGGGGGTGGCGCTGGGACCGCTGTTTAATCTGGGTGCTAATAAAACGGCACAGCGGTATAAAAAgggggggg                                                  >  1:240895/1‑148 (MQ=255)
                                                                                   gcgcTATGGCCCGGCGGGTGTGGCGGGGTTGGCGGTAAAAATAACGGCTCGCCCACCGCCGGAAAAGTTTGTGGGTGGGGGGCTGGGAACGCCCGTTAATCGGGCGGCCTATAAAATGGGACAGCGGTTTAAAAAgggggggggggtt                                             >  2:217488/1‑143 (MQ=255)
                                                                                   gcgcTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTc                                            <  1:136533/149‑1 (MQ=255)
                                                                                   gcgcTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCg                                           <  1:168261/150‑1 (MQ=255)
                                                                                        aTGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCaaa                                       <  2:283571/149‑1 (MQ=255)
                                                                                          ggCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAAtt                                   <  2:11337/151‑1 (MQ=255)
                                                                                            cccGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGt                                 <  1:619278/151‑1 (MQ=255)
                                                                                                      gtgGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCa                        <  1:229538/150‑1 (MQ=255)
                                                                                                      gtgGCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAg                       <  2:410819/151‑1 (MQ=255)
                                                                                                         gCGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCg                    <  2:776202/151‑1 (MQ=255)
                                                                                                          cGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGt                   <  1:781217/151‑1 (MQ=255)
                                                                                                          cGGGCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGt                   <  2:730730/151‑1 (MQ=255)
                                                                                                            ggCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGtat                 >  1:150633/1‑151 (MQ=255)
                                                                                                            ggCTGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGtat                 <  2:95560/151‑1 (MQ=255)
                                                                                                               tGGGCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAat              >  2:426239/1‑151 (MQ=255)
                                                                                                                gggCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGGTGTGGCGCTGGAACCGCTGTTAAGGCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGATGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATaa               <  1:15427/149‑1 (MQ=255)
                                                                                                                gggCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATaa               <  1:190201/149‑1 (MQ=255)
                                                                                                                gggCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAatg             <  2:210369/151‑1 (MQ=255)
                                                                                                                gggCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAat              <  1:14392/150‑1 (MQ=255)
                                                                                                                  gCGATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAatgg            >  1:390802/1‑150 (MQ=255)
                                                                                                                    gATAAAATCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAatggatg         <  2:240914/151‑1 (MQ=255)
                                                                                                                          aTCAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGATGGCAAtt   <  1:187829/151‑1 (MQ=255)
                                                                                                                            cAGCTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGATGGCAAttt  <  1:597388/150‑1 (MQ=255)
                                                                                                                                |                                                                                                                                                    
ATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGATAAAATCAGTTGCTCGCCGACCGCTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGATGGCAATTT  >  NC_000913/1369393‑1369666

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