Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	986,410	A→C	8.9%	S191S (TCT→TCG)	ompF ←	outer membrane porin F

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	986,410	0	A	C	8.9%	104.5 / 2.3	45	S191S (TCT→TCG)	ompF	outer membrane porin F
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base C (3/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

ACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAG  >  NC_000913/986277‑986552                                                                                                                                       |                                                                                                                                               aCTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTCTAACGGCga                                                                                                                                       <  1:687416/143‑2 (MQ=255) aCTGTTCAGCTTGTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGAACCATCAGCCACAACGGTACCACAGCCTTCGTATTGGTAGGCGATCGCACCGACCAAC‑CCGTGTCCGCTCGAacgacgc                                                                                                                                       >  2:338493/1‑142 (MQ=255)     ttCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACCGCGTGCAg                                                                                                                                   >  1:486496/1‑143 (MQ=255)           ttttttACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCAc                                                                                                                             >  1:544733/1‑143 (MQ=255)                      tGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTa                                                                                                                  <  1:579073/143‑1 (MQ=255)                        ccAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTCAcc                                                                                                                >  2:534255/1‑143 (MQ=255)                        ccAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTGCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGACCGGCGTGCAGTGTCACGCTCGTTTTTAcc                                                                                                                >  1:342496/1‑143 (MQ=255)                         cAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTAccc                                                                                                               >  1:575494/1‑143 (MQ=255)                               gTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTCTTTACCCAGGTAc                                                                                                         >  2:608921/1‑143 (MQ=255)                                        ctggaaGGTTGGTACGGTCAGCTGGACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGTTAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCg                                                                                                <  2:570081/138‑1 (MQ=255)                                         ttGCAGGTTGGTACGGTCAGATGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACACTCGTTTCTACGCAGGAACTGAACCGcaa                                                                                               >  2:241352/1‑141 (MQ=255)                                                ttGGTACGGTCAGCTGCACCATAAGCACCACCGCTACGTTAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGGTCAg                                                                                        >  1:150726/1‑143 (MQ=255)                                                 ttgTACGGGCCGCTGCACCATAAGCACCAACGATCCCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAgg                                                                                       <  2:241100/141‑1 (MQ=255)                                                  ggTACGGTCAGCTGCACCATAAGCACCAACGATTGAACAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                                                                      <  1:610920/143‑1 (MQ=255)                                                      cGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACACAGGTACTGAACAGCGAAGTTCAGGCCaca                                                                                  >  2:164543/1‑141 (MQ=255)                                                         tCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGGCATCAAc                                                                               <  2:443582/143‑1 (MQ=255)                                                           aGCTGCACCATCAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGCACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCa                                                                             <  2:490929/143‑1 (MQ=255)                                                             cTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CAAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGCTCAGGCCATCGACCAGa                                                                           >  1:314424/1‑143 (MQ=255)                                                              tGCACCATAAGCACCAACGATACCAAAGCTTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAc                                                                          <  2:114606/143‑1 (MQ=255)                                                               gCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGGTCGGCGAAGTTCAGGCCATCAACCAGAcc                                                                         >  1:398524/1‑143 (MQ=255)                                                                cACCATAAGCACCAGCGATACCAAAGCCTGCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGACCGGCCTCCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCa                                                                        <  2:726647/143‑1 (MQ=255)                                                                        gCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGtt                                                                <  1:49770/143‑1 (MQ=255)                                                                         cACCAACGACGCGATGGCCTTCGTATGCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTg                                                               <  1:696710/143‑1 (MQ=39)                                                                           cgaaCGAGACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGa                                                             <  1:435239/141‑1 (MQ=255)                                                                              aCGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTCCAGGCCATCAACCAGACCAAAGCAGTTggcgcc                                                          >  2:213818/1‑139 (MQ=255)                                                                                   cccACAGCCCTCTTACTCGTAGCTGACAGCACCG‑CCAACCCCGTCGCCGTGCGACCGTGGCGCAGGCTGACGCGCGGCTTTACCCACGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGAt                                                     <  1:20175/142‑1 (MQ=19)                                                                                      aaaGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAgg                                                  <  2:472608/143‑1 (MQ=255)                                                                                      aaaGCCCTAGTACACGTAGCTGATAGACCCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAgg                                                  <  1:247257/143‑1 (MQ=38)                                                                                       aaGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCTTTTTTCCCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGt                                                 <  1:28834/143‑1 (MQ=255)                                                                                             gcgtatacgtaGCGGATAGAACCG‑CCAACACCGTCGCCGTTAGAGCGGCGTCCAGTCTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAc                                           <  1:178548/142‑1 (MQ=38)                                                                                                      taGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGGCGGAGTTACGATAGGTAGCAACGCCGCCAac                                  >  1:249281/1‑143 (MQ=255)                                                                                                          tGATAGAACCT‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGa                              <  1:173716/143‑1 (MQ=255)                                                                                                          tGATAGAACCG‑CCAACACCGTCGCCGTTAGACCGCCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGCACAGCGAAGTTCAGGCCATCAACCAGACCAATGAAGTTGGAGTTACGATAGGTAGCAACGCCACCAACACGa                              >  1:281924/1‑143 (MQ=255)                                                                                                          gaatggcaCCG‑CCGACACCGTCGCCGGTAGAACGGCGCGCAGCGGCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTAAGATAGGTAGCAACGCCGCCAACACGa                              <  1:629977/136‑1 (MQ=37)                                                                                                            aTAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGAcc                            >  1:90257/1‑143 (MQ=255)                                                                                                                 aCCG‑CCAACACCGTCGCCGTTAGAACGACGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGGAGCAACGCAGCCAACACGACCAACga                       >  2:446878/1‑143 (MQ=255)                                                                                                                  ccG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCCAAGAAGTTGGAGTTGCGATAGGTAGCAACGCCGCCAACACGACCAACgaa                      >  2:572660/1‑143 (MQ=255)                                                                                                                  ccG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACgaa                      >  1:151426/1‑143 (MQ=255)                                                                                                                        aaCACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGt                 >  2:218002/1‑143 (MQ=255)                                                                                                                            tgcacGCCGTTCTAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATc             >  2:572869/5‑143 (MQ=255)                                                                                                                              cacGCCGTTCTAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGATGGAGTCACGATAGGTAGCAACGCCACCAACACGACCAACGAAGAAGTCATCGc           >  2:271862/3‑143 (MQ=38)                                                                                                                               tCGCCGTGAGGACGGCGTGCAGTGTCACGCTCGGTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACATGACCAACGAAGAAGTCATCGCt          <  2:212984/143‑1 (MQ=255)                                                                                                                                  ccGTGAGAACGGAGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTa       <  2:425776/143‑1 (MQ=255)                                                                                                                                   ctttAGAACGGCGTGCAGTGTCACGGTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCATCCCGTCCAAAGTAGGTGGAGTTACGTTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTAt      <  2:197388/141‑1 (MQ=255)                                                                                                                                     ttAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGCAGAAGTCATCCCTGTAtcc    >  1:206741/1‑141 (MQ=255)                                                                                                                                       aGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGCTAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAg  <  1:748999/143‑1 (MQ=255)                                                                                                                                       |                                                                                                                                               ACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCG‑CCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAG  >  NC_000913/986277‑986552

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