Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,020,166 C→T 14.4% intergenic (‑113/+244) ompA ← / ← sulA outer membrane porin A/cell division inhibitor SulA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,020,1660CT14.4% 52.6 / 4.0 28intergenic (‑113/+244)ompA/sulAouter membrane porin A/cell division inhibitor SulA
Reads supporting (aligned to +/- strand):  ref base C (11/13);  new base T (2/2);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

TCGCGATAGCTGTCTTTTTCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGAAAAAAAAGTCTTGTATAAGGTATGTTTAATCTTTTTTGTCAGCGACAATTTACAGAAG  >  NC_000913/1020032‑1020303
                                                                                                                                      |                                                                                                                                         
tCGCGATAGCTGTCTTTTTCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACGGGGCAATCTGCTGATTCATAACTCTACCGATATCTTCGGCTTATgggc                                                                                                                                   >  2:140104/1‑140 (MQ=255)
           gTCTTTTTCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTgg                                                                                                                        <  1:512367/143‑1 (MQ=255)
                     ttttttGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTCTGCCGAGCCCCCCTGGCGCTGTAAAg                                                                                                              >  1:219446/1‑143 (MQ=255)
                     tttgttGCGCCTCGTTCTCATACAAAAGACGCCATGAATAGCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGCGATCTGCTCAAGATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAg                                                                                                              <  1:282054/143‑1 (MQ=255)
                      tttttGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCCCCCCTGGCGATGTAAAGt                                                                                                             >  2:170782/1‑143 (MQ=255)
                          tGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTAc                                                                                                         <  1:168488/143‑1 (MQ=255)
                            ctccTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACaa                                                                                                       <  2:142343/141‑1 (MQ=255)
                                   ttATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTGCAACGTAGtt                                                                                                >  1:84370/1‑143 (MQ=255)
                                           cAACATACGCCATGTGTATCTCCATCGAGATAACACGGTTAAATCCGTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACtt                                                                                        <  2:399218/143‑1 (MQ=255)
                                           cAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTCGTTGCAAACtt                                                                                        >  1:339605/1‑143 (MQ=255)
                                                         aaTATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTc                                                                          >  1:323807/1‑143 (MQ=255)
                                                              ctcCTACGAGATAACACGGTTAAATCCTTCCCCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCCCCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCa                                                                     <  2:532621/143‑1 (MQ=255)
                                                                 cAACGAGATAACACGGTTAAATCCTTCCCCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGc                                                                  <  1:153767/143‑1 (MQ=255)
                                                                  tgcGAGATAACACGTTCAAATCCTCCCCCGGGGGCTCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCa                                                                 <  2:86414/141‑1 (MQ=255)
                                                                  aaCGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGGTGTAACGTCTACAACGTAGTTGAACACTTACCAGTGTGAACTCCGTCAGGCa                                                                 >  1:30210/1‑143 (MQ=255)
                                                                  aaCGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAACACTTACAAGTGTGAACTCCGTCAGGCa                                                                 >  2:506064/1‑143 (MQ=255)
                                                                         tAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGaaa                                                          <  1:233295/143‑1 (MQ=255)
                                                                                              cGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGAAAAAAAAGTCTTGTATAAGGTAt                                     >  1:362346/1‑143 (MQ=255)
                                                                                                    aTCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGAAAAAAAAGTCTTGTATAAGGTATGTTTaa                               <  1:13766/143‑1 (MQ=255)
                                                                                                        gCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGAAAAAAAAGTCTTGTATAAGGTATGTTTAATCtt                           <  1:170782/143‑1 (MQ=255)
                                                                                                         cTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGAAAAAAAAGTCTTGTATAAGGTATGTTTAATCttt                          >  1:549398/1‑143 (MQ=255)
                                                                                                              cagtAACTCTACGGATATCTGCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGAAAAAAAAGTCTTGTATAAGGTATGTTTAATCTTTTTTGt                     <  2:391769/140‑1 (MQ=255)
                                                                                                                     tctACCGATATCTTCGGTTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGaaaaaaaa                                                     >  1:49591/1‑104 (MQ=255)
                                                                                                                     tctACCGATATCTTCGGTTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGaaaaaaaa                                                     >  1:319085/1‑104 (MQ=255)
                                                                                                                     tctACCGATATCTTCGGTTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGaaaaaaaa                                                     <  2:319085/104‑1 (MQ=255)
                                                                                                                     tctACCGATATCTTCGGTTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGaaaaaaaa                                                     <  2:49591/104‑1 (MQ=255)
                                                                                                                                 ttCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGCACTCCGTCAGGCATATGAAAAAAAAGTCTTGTATAAGGTATGTTTAATCTTTTTTGTCAGCGACCATTTACagaag  >  1:213423/1‑143 (MQ=255)
                                                                                                                                 ttCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGAAAAAAAAGTCTTGTATAAGGTATGTTTAATCTTTTTTGTCAGCGACAATTTACagaag  <  2:323807/143‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                         
TCGCGATAGCTGTCTTTTTCATTTTTTGCGCCTCGTTATCATCCAAAATACGCCATGAATATCTCCAACGAGATAACACGGTTAAATCCTTCACCGGGGGATCTGCTCAATATTAACTCTACCGATATCTTCGGCTTATGCCGAGCACCCCTGGCGATGTAAAGTCTACAACGTAGTTGAAAACTTACAAGTGTGAACTCCGTCAGGCATATGAAAAAAAAGTCTTGTATAAGGTATGTTTAATCTTTTTTGTCAGCGACAATTTACAGAAG  >  NC_000913/1020032‑1020303

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