Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	2,179,015	C→G	12.2%	T178T (ACC→ACG)	U712_10905 →	UPF0750 membrane protein ypjC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,179,015	0	C	G	12.2%	115.0 / 8.1	41	T178T (ACC→ACG)	U712_10905	UPF0750 membrane protein ypjC
Reads supporting (aligned to +/- strand):  ref base C (24/12);  new base G (3/2);  total (27/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

GAGTTGATATCATCGCCCGTTTGGTCAATAAGTATTTTGGGATTCCGATGGGCAGAACGATGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCTTGAAGAAATGGAACGAGGCGT  >  CP006881/2178882‑2179152                                                                                                                                      |                                                                                                                                          gAGTTGATATCATCGCCCGTTTGGTCAATAAGTATTTTGGGATTCCGATGGGCAGAACGATGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCtggt                                                                                                                                      <  7:39267/139‑1 (MQ=255)      gATATCATCGCCCGTTTGGTCAATAAGTATTTTGGGATTCCGATGGGCAGAACGATGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCt                                                                                                                                  <  1:94354/138‑1 (MQ=255)                   gTTTGGTCAATAAGTATTTTGGGATTCCGATGGGCAGAACGATGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGc                                                                                                                    <  4:226925/139‑1 (MQ=255)                                          aTTCCGATGGGCAGAACGATGTTTGCTTGCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAgg                                                                                             <  3:336342/139‑1 (MQ=255)                                           ttCCGATGGGCAGAACGATGTTTGCTTTCGATGCCTGCGTCATCATTTTGTCTTTGCTTACTTATCTCTCCTATAAAGAAGCGATGTATACGCTGGTAGCTGTATTTGTAGGGGCAAGATTAATCGATTTTATTCAgg                                                                                             >  2:319607/1‑138 (MQ=255)                                           ttCCGATGGGCAGAACGATGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAgg                                                                                             >  1:221486/1‑138 (MQ=255)                                           ttCCGATGGGCAGAACGATGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAgg                                                                                             >  5:117917/1‑138 (MQ=255)                                           ttCCGATGGGCAGAACGATGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAgg                                                                                             >  3:264023/1‑138 (MQ=255)                                                        aaCGATGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATTCCCTGGTGGCTGTTTTTGTAGCTGCAAGATTTATTGATTTTATTGAGGAAGGCGGATACgct                                                                               >  1:305599/1‑139 (MQ=255)                                                        aaCGATGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACgct                                                                               >  1:108110/1‑139 (MQ=255)                                                        aaCGATGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACgct                                                                               >  7:19033/1‑139 (MQ=255)                                                             tGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGAtt                                                                                                      <  1:330633/111‑1 (MQ=255)                                                             tGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGAtt                                                                                                      >  2:330633/1‑111 (MQ=255)                                                                      tCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACgctgct                                                                            >  1:42944/1‑128 (MQ=255)                                                                      tCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACgctgct                                                                            <  2:42944/128‑1 (MQ=255)                                                                         aTGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGc                                                                                        >  1:186835/1‑113 (MQ=255)                                                                         aTGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGc                                                                                        <  2:186835/113‑1 (MQ=255)                                                                                     tcatTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGa                                                   >  3:126404/1‑138 (MQ=255)                                                                                     tcatTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGa                                                   >  5:288379/1‑138 (MQ=255)                                                                                                ttGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTCTATTTGTAGCTGCAAGATTAATATATTTTATTAAGGAAGTCTGAAACGCTGCTAAAGGCGCTACGATTATCTATTCGAAAAACGAGCTGa                                       >  3:192374/1‑139 (MQ=255)                                                                                                 tGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGAt                                           >  1:288807/1‑134 (MQ=255)                                                                                                 tGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGAt                                           <  2:288807/134‑1 (MQ=255)                                                                                                          aTCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGaaa                               <  1:207998/137‑1 (MQ=255)                                                                                                           tctctcCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGAtt                                                                                                      <  3:262069/65‑1 (MQ=255)                                                                                                           tctctcCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGAtt                                                                                                      >  4:262069/1‑65 (MQ=255)                                                                                                                  tACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAATATTCtt                      >  6:324069/1‑138 (MQ=255)                                                                                                                  tACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCtt                      >  3:318221/1‑138 (MQ=255)                                                                                                                    cAAAGAAGCGATGTATACCCTTGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGCAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCTTg                     <  4:248359/137‑1 (MQ=255)                                                                                                                    cAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCTTgaa                   >  5:12411/1‑139 (MQ=255)                                                                                                                         aagCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCTTGAAGAAAt              >  2:204818/1‑139 (MQ=255)                                                                                                                         aagCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCTTGAAGAAAt              >  5:251789/1‑139 (MQ=255)                                                                                                                         aagCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCTTGAAGAAAt              >  6:150788/1‑139 (MQ=255)                                                                                                                         aagCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCTTGAAGAAAt              >  7:137743/1‑139 (MQ=255)                                                                                                                            cGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCTTGAAGAAATgg            >  4:104191/1‑138 (MQ=255)                                                                                                                              aTGTATACGCTGGTAGCTGTATTTGTAGCGGCAAGATTAATCGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATCATCTCTTCGaaaaa                                               <  7:269374/101‑1 (MQ=255)                                                                                                                              aTGTATACGCTGGTAGCTGTATTTGTAGCGGCAAGATTAATCGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATCATCTCTTCGaaaaa                                               >  8:269374/1‑101 (MQ=255)                                                                                                                              aTGTATACGCTGGTAGCTGTATTTGTAGCGGCAAGATTAATCGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATCATCTCTTCGaaaaa                                               >  6:237641/1‑101 (MQ=255)                                                                                                                              aTGTATACGCTGGTAGCTGTATTTGTAGCGGCAAGATTAATCGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATCATCTCTTCGaaaaa                                               <  5:237641/101‑1 (MQ=255)                                                                                                                              aTGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCTTGAAGAAATGGAAc         <  6:251789/139‑1 (MQ=255)                                                                                                                                   taCCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCTTGAAGAAATGGAACGAGGc    >  8:115138/1‑139 (MQ=255)                                                                                                                                   taCCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCTTGAAGAAATGGAACGAGGc    >  3:263545/1‑139 (MQ=255)                                                                                                                                      ccTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGAGGATACGCTGCTAAAGGCGCTATGTTTATCTCTTCGAAATACGATCTGATTCAGAAAAAGATTCTTGAAGAAATGGAACGAgtgtt  <  4:264023/138‑5 (MQ=255)                                                                                                                                      |                                                                                                                                          GAGTTGATATCATCGCCCGTTTGGTCAATAAGTATTTTGGGATTCCGATGGGCAGAACGATGTTTGCTTTCGATGCCTGCGTCATCATTTTATCTTTGCTTACTTATCTCTCCTACAAAGAAGCGATGTATACCCTGGTGGCTGTATTTGTAGCTGCAAGATTAATTGATTTTATTCAGGAAGGCGGATACGCTGCTAAAGGCGCTACGATTATCTCTTCGAAAAACGATCTGATTCAGAAAAAGATTCTTGAAGAAATGGAACGAGGCGT  >  CP006881/2178882‑2179152

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