Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	7,925	+A	coding (35/1431 nt)	yaaJ ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	7,922	1	.	A	94.7%	113.0 / ‑0.5	38	coding (38/1431 nt)	yaaJ	putative transporter
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (15/21);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑C‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTG  >  NC_000913/7793‑8056                                                                                                                                     |                                                                                                                                       ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACg                                                                                                                                <  1:252544/139‑1 (MQ=255) ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACg                                                                                                                                <  1:48194/139‑1 (MQ=255)  cTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGc                                                                                                                               >  2:728087/1‑139 (MQ=255)   tGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCt                                                                                                                              <  2:781342/139‑1 (MQ=255)         tgGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAAt                                                                                                                        >  2:110262/1‑139 (MQ=255)           gATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGa                                                                                                                      <  2:494712/139‑1 (MQ=255)                     aTTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGaaaaa                                                                                                            >  2:704596/1‑139 (MQ=255)                     aTTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGaaaaa                                                                                                            >  2:671538/1‑139 (MQ=255)                      tttttAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGaaaaaa                                                                                                           >  2:473386/1‑139 (MQ=255)                                  ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑C‑aaaaggac                                                                                                                                <  1:163493/105‑6 (MQ=255)                                  ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑C‑aaaaggac                                                                                                                                >  2:163493/1‑100 (MQ=255)                                      cAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                           <  1:801820/139‑1 (MQ=255)                                      cAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                           <  2:38594/139‑1 (MQ=255)                                             gCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                                                                <  1:705339/127‑1 (MQ=255)                                             gCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGCGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGccatgca                                                                                                >  2:705339/1‑127 (MQ=255)                                                       cGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑C‑aaaagga                                                                                                                                 <  1:430794/83‑5 (MQ=255)                                                       cGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑C‑aaaagga                                                                                                                                 >  2:430794/1‑79 (MQ=255)                                                               ctcaAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACtt                                                                  >  1:480188/3‑139 (MQ=255)                                                                 caAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                                <  1:38836/139‑1 (MQ=255)                                                                    aTCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCtg                                                             >  1:730807/1‑139 (MQ=255)                                                                      ccGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCtgtg                                                           <  2:730807/139‑1 (MQ=255)                                                                       cGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTgtgt                                                          >  2:547703/1‑139 (MQ=255)                                                                         gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCTGTCGCGATGACTTTCCTGTTTaa                                                        >  1:138551/1‑139 (MQ=255)                                                                         gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa                                                        <  1:73899/139‑1 (MQ=255)                                                                         gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa                                                        >  1:95147/1‑139 (MQ=255)                                                                         gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa                                                        <  2:25352/139‑1 (MQ=255)                                                                           gCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAAc                                                      <  2:88906/139‑1 (MQ=255)                                                                              aaaaGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACAGATCCAC‑A‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGtt                                                   >  1:183815/1‑139 (MQ=255)                                                                                 aGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTAcc                                                <  1:728087/139‑1 (MQ=255)                                                                                           caacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTATCCTGTGTAAGCGTTACCAATTGTTTaa                                      <  1:83273/139‑1 (MQ=255)                                                                                             acCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaaga                                    <  1:547703/139‑1 (MQ=255)                                                                                             acCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaaga                                    <  2:183815/139‑1 (MQ=255)                                                                                                 gcgtCGAAGAGCAGGTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCAGATCCCTCTTATTGCCGGACGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGta                                >  1:515342/1‑139 (MQ=255)                                                                                                               gTAAATCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAgg                  <  1:671538/139‑1 (MQ=255)                                                                                                                  aaTCATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTAc               <  1:177302/139‑1 (MQ=255)                                                                                                                   aTCATTACCGATCC‑CAA‑AAAGGACGCTGTTAATGAAGGAGAAGAAATCTGGCATGCATATCCTTCGTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTAACAATTGAGTAAGAAGTATATACGCTACGAGCTACt              >  1:545178/1‑139 (MQ=255)                                                                                                                     cATTACCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTg            <  1:168311/139‑1 (MQ=255)                                                                                                                         aCCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATaa        >  2:536099/1‑139 (MQ=255)                                                                                                                         aCCGATCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATaa        >  2:721273/1‑139 (MQ=255)                                                                                                                            atccc‑c‑a‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTc    >  1:540876/4‑139 (MQ=255)                                                                                                                             tccc‑c‑a‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaa                                      >  1:95656/3‑104 (MQ=255)                                                                                                                             aTCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACtt     <  2:315402/138‑1 (MQ=255)                                                                                                                             aTCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACtt     >  1:315402/1‑138 (MQ=255)                                                                                                                              tCC‑C‑CAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaa                                      <  2:95656/104‑1 (MQ=255)                                                                                                                              ccc‑c‑a‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACg                     <  2:180436/119‑1 (MQ=255)                                                                                                                              ccc‑c‑a‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTg  >  2:260809/2‑139 (MQ=255)                                                                                                                                     |                                                                                                                                       CCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑C‑C‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTG  >  NC_000913/7793‑8056

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑