Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,289,257	A→C	V78G (GTA→GGA)	ybbP ←	hypothetical membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,289,257	0	A	C	97.1%	112.1 / NA	35	V78G (GTA→GGA)	ybbP	hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (11/23);  total (11/24)
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

TATTGCACAGCCTTTGAAACCGATTGAATCAATTTCTCTTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAG  >  CP000730/2289110‑2289399                                                                                                                                                    |                                                                                                                                               tatTGCACAGCCTTTGAAACCGATTGAATCAATTTCTCTTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATcccc                                                                                                                                             >  2:214342/1‑151 (MQ=255)     gCACAGCCTTTGAAACCGATTGAATCAATTTCTCTTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCa                                                                                                                                         >  1:120934/1‑151 (MQ=255)     gCACAGCCTTTGAAACCGATTGAATCAATTTCTCTTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCGACCACGACGATTGCGGTCTAACGCACTTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCa                                                                                                                                         >  1:270873/1‑151 (MQ=255)       acaGCCTTTGAAACCGATTGAATCAATTTCTCTTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCAtt                                                                                                                                       <  2:277631/151‑1 (MQ=255)        caGCCTTTGAAACCGATTGAATCAATTTCTCTTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTg                                                                                                                                      >  2:199023/1‑151 (MQ=255)              ttGAAACCGATTGAATCAATTTCTCTTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAAc                                                                                                                                <  1:57576/151‑1 (MQ=255)              ttGAAACCGATTGAATCAATTTCTCTTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAAc                                                                                                                                <  1:92782/151‑1 (MQ=255)                      gATTGAATCAATTTCTCTTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGa                                                                                                                        >  2:440514/1‑151 (MQ=255)                        ttGAATCAATTTCTCTTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGaat                                                                                                                      <  2:303553/151‑1 (MQ=255)                              cAATTTCTCTTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTa                                                                                                                >  1:72296/1‑151 (MQ=255)                                 tttCTTTTCATTATTAATATAAGTATTAGAAGTATAGCGTTATAAAAAGCTACCTCTACCAAGTTGTTCTAACGCAAATCTAATTTCTTGTTAGAATATTAATATTAAAGTTAATAACACCAATTGAATAACGATATCGAATAATTTAGAt                                                                                                             <  2:92733/151‑1 (MQ=255)                                    ctctTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGtt                                                                                                          >  2:412542/1‑151 (MQ=255)                                           tCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCa                                                                                                   >  2:235790/1‑151 (MQ=255)                                           tCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCa                                                                                                   <  2:120934/151‑1 (MQ=255)                                            cTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCaa                                                                                                  >  1:387860/1‑151 (MQ=255)                                             tttACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAg                                                                                                 <  1:425849/151‑1 (MQ=255)                                             tttACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAg                                                                                                 <  2:113937/151‑1 (MQ=255)                                               tACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGtt                                                                                               <  1:358301/151‑1 (MQ=255)                                               tACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGtt                                                                                               <  1:358323/151‑1 (MQ=255)                                                      cGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATc                                                                                        >  2:74734/1‑151 (MQ=255)                                                       gTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCGCCCCCATTGAATAACGATATCGACTAATTTAGATGTTGCAGTCAAGTTCAATATCa                                                                                       <  2:388563/151‑1 (MQ=255)                                                        tATTAGAAGTATAGCGTTTTAAAAAGCTAATTCTACCAAGTTGTTCTAACGAACGTCTAATTTCTGGTTGGACTATTACTATTAAAGCTCATCCCCCCCATTGAATAACGATTTCGAATAATTTAGATGTTGAAATGAAGTTCAATATCAt                                                                                      <  2:193116/151‑1 (MQ=255)                                                                       gTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGCATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGAcc                                                                       <  2:865/151‑1 (MQ=255)                                                                         tttAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCaa                                                                     <  2:37223/151‑1 (MQ=255)                                                                             aaaaGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGCTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAAc                                                                 <  2:266218/151‑1 (MQ=255)                                                                                      ctctACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTaa                                                        <  2:14476/151‑1 (MQ=255)                                                                                           ccAAGTTGTTCTAACGCACGTCTAATNTCTGGTTGNAATNTTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATcc                                                   >  1:421020/1‑151 (MQ=255)                                                                                              aGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCtt                                                <  1:262155/151‑1 (MQ=255)                                                                                               gttgttCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCttt                                               <  1:62886/151‑1 (MQ=255)                                                                                                    tCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTa                                          >  1:194105/1‑151 (MQ=255)                                                                                                                cTAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCtt                              <  1:396510/151‑1 (MQ=255)                                                                                                                 tAATTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCttt                             <  1:412542/151‑1 (MQ=255)                                                                                                                  aaTTTCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTa                            <  2:171161/151‑1 (MQ=255)                                                                                                                      tCTGGTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTc                        <  2:320477/151‑1 (MQ=255)                                                                                                                         ggTTGGAATATTACTATTAAAGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCt                     <  2:419783/151‑1 (MQ=255)                                                                                                                                            aaGCTAATCCCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGagaag  <  1:301484/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               TATTGCACAGCCTTTGAAACCGATTGAATCAATTTCTCTTCATCTTTACTATACGTATTAGAAGTATAGCGTTTTAAAAAGCTACCTCTACCAAGTTGTTCTAACGCACGTCTAATTTCTGGTTGGAATATTACTATTAAAGCTAATACCCCCCATTGAATAACGATATCGAATAATTTAGATGTTGCAGTCAAGTTCAATATCATACTTATCTGCTGACCAATAACAATTACTAATATCCCTTTAAGTAATTGTATCGCTTTAGTTCCCTTAAAGACCGTGATGAGAAG  >  CP000730/2289110‑2289399

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