Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,652,331 T→C intergenic (+2346/+397) ydfD → / ← ynfP lysis protein/protein YnfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,652,3310TC96.7% 96.7 / ‑5.6 30intergenic (+2346/+397)ydfD/ynfPlysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (11/18);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCC  >  NC_000913/1652087‑1652571
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
aTGGTGTGTGTCTGTTTGTACTGGATGTCGCNCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAAGTTAGCAAAAGACCGTTTCTCGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTACGGAGTCAGATGACATTAACGAGACTGAGCGAAGAGCATATCATTGATTTTC‑ACCTCagag                                                                                                                                                                                                                                            <  2:218884/251‑1 (MQ=255)
               ttATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGACAGCAATCCACGCCCTCCGGAGTCAGATGAAATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCACAGAGAATGGcagaaca                                                                                                                                                                                                                              >  1:420117/1‑250 (MQ=255)
                 ataCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGcagaacaga                                                                                                                                                                                                                            <  2:349273/250‑1 (MQ=255)
                      ggAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGACAAATGGAACAATGATGATCAGGAGAAATTTCGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGACCAATAGGGGCATTCATCTTATTAAGCCAGCCATCGACGCCCTCCGGCGTCAGAAGACCATAAAGAGACTGAGAAAAGAGCATATAATTGACGGTT‑ACCTCCGAGAGTATGGCAGAACAGaaacac                                                                                                                                                                                                                       >  1:41711/1‑247 (MQ=255)
                                                aCTTGGATGGGAGGACATAGATCTGAAAAATNGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTg                                                                                                                                                                                            <  2:225179/251‑1 (MQ=255)
                                                                         aaaaaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCCAACGGGGCATTTATCTTATTAAGCCAGCAATCGACGCTCTCCCGCGTAAGATGACATTAACGAGACTGAGCAAAAAGCATATCATTGATGTTC‑ACCCCCGAGAGCATGGCACAACCGAAAAAAAAAAATGCACCCTTTTTTTTTAACCTGCCATGTCAGCGAGAGTCAAAAAttt                                                                                                                                                                   >  2:170918/1‑250 (MQ=255)
                                                                            aaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATgg                                                                                                                                                                <  2:307575/251‑1 (MQ=255)
                                                                                    aaTGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATTCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGGGAccctt                                                                                                                                                        >  1:358444/1‑248 (MQ=255)
                                                                                                    aaTTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAAt                                                                                                                                        <  1:406524/251‑1 (MQ=255)
                                                                                                                 aCCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTgg                                                                                                                           <  2:156039/251‑1 (MQ=255)
                                                                                                                    gTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGTGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGAc                                                                                                                        <  1:60311/251‑1 (MQ=255)
                                                                                                                      ttCACGGTACCAAAAACACAGGCGGGGACCAATCGGGTCCTTCATCTTATTAAAACAGCAATCGAACCTTTCCGGAGTCAGATGACATTAACGAGACTGAGCCAAGAGAATATTATTGATGTTCCACT‑CAGAGGGGATGGCAGAAAAAATAAAAAAAAATGCACCTTTGTTTTTTAACCTGAAGGGTCAGCGAGAGTAAAAAATTATGGTGACCCTTTTACCGTTGACTCAAAAAAGCAGATGTggggagg                                                                                                                      >  1:464346/1‑247 (MQ=255)
                                                                                                                           ggTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTTTTAAGCCAGCAATCGACGCTCTCCGGAGTCAGAAGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCCGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGGGAACATTTTACCGTTTACTCAATAAGGCAGATGTGGGACGCcgcg                                                                                                                  >  1:57635/1‑250 (MQ=255)
                                                                                                                            gTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGAt                                                                                                                <  2:352010/251‑1 (MQ=255)
                                                                                                                             tACCAAAAACACAGGCGGGAACCCATCGGGGCCTTCATCTTATTTAACCAGCAATCGACGCCCTCCGGAGGCAGATGACATTAAAGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGGATGGCAGAACAAAAAAAAAAAAATGAACCTTTGTTTTTCAAACTGAAGTGTCAGCGAGAGTAAAAAATTATTGTGACCATTTTTCCGTTTTCTCAAAAAGGCAGATGTTGGACGccacgat                                                                                                                >  2:298944/1‑244 (MQ=255)
                                                                                                                              aCCAAAAACACAGGCGGGAACCAATCGGGTNATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATa                                                                                                               <  2:228506/250‑1 (MQ=255)
                                                                                                                                cAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGTCGCAGCGATaaaa                                                                                                            <  2:159228/251‑1 (MQ=255)
                                                                                                                                      cacaGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGcc                                                                                                      <  1:250596/251‑1 (MQ=255)
                                                                                                                                           gCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGAATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACt                                                                                                 >  1:124953/1‑251 (MQ=255)
                                                                                                                                                         ggTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAAt                                                                                   >  1:147351/1‑251 (MQ=255)
                                                                                                                                                                     tattaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGATCATATCATTTATGTTC‑ACCTCAGAGAGGATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAATTCATATCAGTc                                                                         <  2:464346/249‑1 (MQ=255)
                                                                                                                                                                                        gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgctg                                                    >  2:147237/1‑251 (MQ=255)
                                                                                                                                                                                               ggAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTg                                             <  1:225947/251‑1 (MQ=255)
                                                                                                                                                                                                          aCATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGc                                  <  1:170918/251‑1 (MQ=255)
                                                                                                                                                                                                                aCGAGACTGAGCAAAGAGCATATCATTGATGNTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAccc                            <  2:216332/251‑1 (MQ=255)
                                                                                                                                                                                                                 cGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCg                           <  1:298944/251‑1 (MQ=255)
                                                                                                                                                                                                                        gatcAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATTCCCCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTTCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAtt                      <  2:337706/246‑1 (MQ=255)
                                                                                                                                                                                                                           cAAAGAGCATATCATTGATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTtat                   >  1:32401/1‑249 (MQ=255)
                                                                                                                                                                                                                                           gATGTTC‑ACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCGTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACGCCGCCATCGAAAATCATATCAGTAGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGGGCTAACCCGGCATTTATAGCAAACCAGATGGGcc  <  2:124953/250‑1 (MQ=255)
                                                                                                                                                                                                                                           gATGTGC‑ACCTCAGAGAGCATTTCATAACAGAAAACCAAAAATGTACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAATTTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTCGGACGCAGCGATAACACGTGCCGGACTCCGCCATCGAAAGACATATCCGTCGATACATACTTATGCCTGCTTGTCGATGACAGCAGGTGCTAACCCGGCATATATGGCAATCAAGTTGGGcc  >  2:116463/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
ATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTC‑ACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCC  >  NC_000913/1652087‑1652571

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