Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	380,786	G→A	A19A (GCC→GCT)	yaiO ←	outer membrane protein YaiO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	380,786	0	G	A	97.1%	115.3 / NA	34	A19A (GCC→GCT)	yaiO	outer membrane protein YaiO
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (11/22);  minor base T (0/1);  total (11/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCGCAAAGGCAATACCCGTTCGGGTTGTCAGCCAGTTATTCCATTTATACCAGACCGCACCCTGACCGCGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGG  >  NC_000913/380546‑381025                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                atCCGCAAAGGCAATACCCGTTCGGGTTGTCAGCCAGTTATTCCATTTATACCAGACCGCACCCTGACCGCGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCTGCGGTAATGGAAGTTAACCCAGCATAAGCg                                                                                                                                                                                                                                        <  1:74749/250‑1 (MQ=255)    cGCAAAGGCAATACCCGTTCGGGTTGTCAGCCAGTTATTCCATTTATACCAGACCGCACCCTGACCGCGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCa                                                                                                                                                                                                                                    <  1:146412/251‑1 (MQ=255)             aaTACCCGTTCGGGTTGTCAGCCAGTTATTCCATTTATACCAGACCGCACCCTGACCGCGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTaa                                                                                                                                                                                                                           <  1:317296/251‑1 (MQ=255)             aaTACCCGTTCGGGTTGTCAGCCAGTTATTCCATTTATACCAGACCGCACCCTGACCGCGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTaa                                                                                                                                                                                                                           <  2:130341/251‑1 (MQ=255)                     tacGGGTTGTCAGCCAGTTATTCCATTTATACCAGACCGCACCCTGACCGCGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCTACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGgcc                                                                                                                                                                                                                   >  2:265448/3‑251 (MQ=255)                        gggTTGTCAGCCAGTTATTCCATTTATACCAGACCGCACCCTGACCGCGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGTCAATGCGCTAAAAATGgccgcc                                                                                                                                                                                                                <  2:339343/251‑1 (MQ=255)                             gTCAGCCAGTTATTCCATTTATACCAGACCGCACCCTGACCGCGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCta                                                                                                                                                                                                            >  2:174015/1‑250 (MQ=255)                                aGCCAGTTATTCCATTTATACCAGACCGCACCCTGAACGAGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCATGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTCAAACTGGCCGCCGCtaata                                                                                                                                                                                                         >  2:33226/1‑250 (MQ=255)                                  ccAGTTATTCCATTTATACCAGACCGCACCCTGACCGCGAGTGGCATTGAAATGTTCATTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGTCAATGCGCTAAAAATGGCCGCCGCTAATAACg                                                                                                                                                                                                      <  1:17916/251‑1 (MQ=255)                                     ggtATTCCATTTATACCAGACCGCACCCTGACCGCGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTcc                                                                                                                                                                                                   <  1:293096/249‑1 (MQ=255)                                        atgccATTTATACCAGCCCGCACCCTGACCGCGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGGAAAAGCGTAGCCCGCGGTAATGGAAGGGAACCCAGCATAAGCGGTCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGt                                                                                                                                                                                                 <  1:241975/247‑1 (MQ=255)                                                    ccAGACCGCACACTGACCGCGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTa                                                                                                                                                                                     >  2:19651/1‑250 (MQ=255)                                                            cACCCTGCCCGCGCGTGGCATTGACATGTTCAGTTTCATAATCACGACGCCCATGGGAAAGATTCAAAAGCAGTGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATCCGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCACTGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCATGt                                                                                                                                                                            >  2:371768/1‑251 (MQ=255)                                                                                                           cgCCCCTGGGAAAGATTTAAAAGCAGGTTTGAGTTTTCAACTTGCGCAACCAGTTCAGCAAACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCTTAAGCGTGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTGTATCATTTTAATCTCCTATACGGATAAGTGCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTAc                                                                                                                             <  2:873/251‑1 (MQ=255)                                                                                                                    gAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACa                                                                                                                     >  1:135164/1‑250 (MQ=255)                                                                                                                             aaaaGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCAcgcgcg                                                                                                           <  1:174015/251‑1 (MQ=255)                                                                                                                                      gttgcgttTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGtt                                                                                                  >  2:79461/1‑251 (MQ=255)                                                                                                                                        ggcgttTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGc                                                                                                <  2:166032/250‑1 (MQ=255)                                                                                                                                         gcgttTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGc                                                                                                <  1:297467/250‑1 (MQ=255)                                                                                                                                               tCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTc                                                                                          <  1:180514/250‑1 (MQ=255)                                                                                                                                                      tCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAaat                                                                                  >  1:113060/1‑251 (MQ=255)                                                                                                                                                                ttCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCgt                                                                                              >  1:62962/1‑229 (MQ=255)                                                                                                                                                                ttCAGCATACGCTAAATTACGGGTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCgt                                                                                              <  2:62962/229‑1 (MQ=255)                                                                                                                                                                                 tACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCTGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTTTTCGGGGAAATAATTCGCCCAGGCAAAATAATAAcgt                                                        <  1:362962/250‑3 (MQ=255)                                                                                                                                                                                 tACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCGCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTATGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTTCGTGTTCTGGGAAATAATTCGCCCAGGCAAAATAATAACtt                                                        <  1:362979/250‑1 (MQ=255)                                                                                                                                                                                         ccATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCaa                                               <  1:67365/251‑1 (MQ=255)                                                                                                                                                                                                                   aTCGTATCCCGCGGTAATGGACGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAATTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGCTCGGGGAAATACTTCGCCCAGGCAAAATAATAACTTGTACGCCATGGCGTGAGCGCTGCCCGATATaaca                      >  2:61869/1‑250 (MQ=255)                                                                                                                                                                                                                   aTCGTAGCCCGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATaacaa                     <  1:314786/251‑1 (MQ=255)                                                                                                                                                                                                                           ccGCGGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAAttt             <  1:79461/251‑1 (MQ=255)                                                                                                                                                                                                                           ccGCGGTAATGGAACTTAACACAGCCTAAGAGCGGAATGCGCTAAAAATGGCAGCCACTAACAACGTCCGTTTTATCATTCTAATCTCCTGTAAGGATAAGCTCTTCTCGGAGTGAATAAGTTCGGGTAGGTAAATTACGTTATACAATCACGCGCCGGACAAGCTGCGTGTTCGGGAAAATAATACGCCCAGGCAAAATAATAACTTGTACGACAAGGCGTGAACGCTGCCCGATATAACAACCTAAtat             >  2:163791/1‑249 (MQ=255)                                                                                                                                                                                                                             gcgGTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTAt           <  1:125547/251‑1 (MQ=255)                                                                                                                                                                                                                               ggTAATGGAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCg         <  2:177543/251‑1 (MQ=255)                                                                                                                                                                                                                                   aTGGAAGTTAACCCAGCATAAGCGGGAAATTCGCTAAAAAGGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATc     <  2:94556/251‑1 (MQ=255)                                                                                                                                                                                                                                      gAAGTTAACCCAGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCgg  <  1:179574/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                ATCCGCAAAGGCAATACCCGTTCGGGTTGTCAGCCAGTTATTCCATTTATACCAGACCGCACCCTGACCGCGAGTGGCATTGAAATGTTCAGTTTCATAATCACGACGCCCCTGGGAAAGATTAAAAAGCAGCGTTGCGTTTTCAACTTTCGCCACCAGTTCAGCATACGCTAAATTACGGTTGCCATGATCGCCAGAATAATCGGTAAAATCGTAGCCCGCGGTAATGGAAGTTAACCCGGCATAAGCGGGCAATGCGCTAAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGG  >  NC_000913/380546‑381025

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