Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	264,098	G→T	L50L (CTC→CTA)	ykfH ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	264,098	0	G	T	97.7%	130.7 / NA	44	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (24/19);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAT  >  NC_000913_3_pae_tpiA/263952‑264241                                                                                                                                                   |                                                                                                                                                cGGGCTCCCAGGGCGAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAg                                                                                                                                               <  1:209515/149‑1 (MQ=255)  gggCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAgg                                                                                                                                              >  1:376513/1‑149 (MQ=255)        ccAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCCGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCg                                                                                                                                        >  1:429514/1‑149 (MQ=255)        ccAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCg                                                                                                                                        >  1:429516/1‑149 (MQ=255)                       tctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCACTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACg                                                                                                                         >  1:174316/1‑149 (MQ=255)                                    gCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtgt                                                                                                            >  2:337491/1‑149 (MQ=255)                                      cccTGGTAGGGTTGCTCATTACAGTCCTCCACTATTGTCATTTCAGGGCGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                          <  1:195340/149‑1 (MQ=255)                                       ccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                         <  2:429516/149‑1 (MQ=255)                                       ccTGGGAGTGTTGCTCATTACAGTTCTCCACTATTTTCATTTCAGGGTGAGGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                         <  2:429514/149‑1 (MQ=255)                                         tGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGa                                                                                                       >  2:310123/1‑149 (MQ=255)                                                   gCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCg                                                                                                                              >  1:108732/1‑116 (MQ=255)                                                   gCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCg                                                                                                                              <  2:108732/116‑1 (MQ=255)                                                   gCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTAc                                                                                             >  1:231782/1‑149 (MQ=255)                                                     tCATTACAGTTCTCCACTATTGTCAGATCCGCTTAATCGTCAGCAGGCAGGGATAGAGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTTCGGCGATATCCGGAACGTCCCTGCCGGGGTAATGACAGACGcggccg                                                                                           >  1:132876/1‑144 (MQ=255)                                                            aGTTCTCCACTATGGTCAGTTCAGGGAGATGCTCATCAGGCAGGTATAGGTCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACt                                                                                    >  1:450967/1‑149 (MQ=255)                                                              ttCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGCCGTATACCGCTAGGACTCCGGCGGTATCCGGAACGTCCCTGCCGGTGTAATGACAGCCGCTCCCGTGcccccgg                                                                                  >  1:78476/1‑144 (MQ=21)                                                                  ccACTATTGTCAGTTCAGGGTGATGCGCATCAGGCAGGTATATGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAtt                                                                              <  2:129713/149‑1 (MQ=255)                                                                   cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGc                                                                                        >  1:312409/1‑138 (MQ=255)                                                                   cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGc                                                                                        <  2:312409/138‑1 (MQ=255)                                                                       aTTTTCAGTTCAGGGTGATGCGCAACAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCg                                                                         <  2:252676/149‑1 (MQ=255)                                                                       aTTGTCAGTTCAGGGTGATGCTCAGCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCg                                                                         <  2:432596/149‑1 (MQ=255)                                                                           tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGc                                                                     >  1:389525/1‑149 (MQ=255)                                                                                          gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTc                                                      <  1:207541/149‑1 (MQ=255)                                                                                           cTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                                     >  2:317975/1‑149 (MQ=255)                                                                                                  ggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCa                                              <  1:494952/149‑1 (MQ=255)                                                                                                       gTATAGGGCCCATTGCGGTCCTTGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATg                                         >  2:145027/1‑149 (MQ=255)                                                                                                         ataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCCCTCGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCt                                       >  2:304900/1‑149 (MQ=255)                                                                                                          taGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                                                                                                    <  1:445672/71‑1 (MQ=255)                                                                                                          taGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                                                                                                    >  2:445672/1‑71 (MQ=255)                                                                                                          taGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTg                                      <  2:9035/149‑1 (MQ=255)                                                                                                             ggCCCAATGCGGTCCTGGCAGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGCAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAAATTCGGGACTCAGGATGCTGGCg                                   <  2:517683/149‑1 (MQ=17)                                                                                                                          ccTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGc                                                                                        >  2:167762/1‑83 (MQ=255)                                                                                                                          ccTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGc                                                                                        <  1:167762/83‑1 (MQ=255)                                                                                                                           cTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTg                     >  1:283940/1‑149 (MQ=255)                                                                                                                                ggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTg                >  1:373146/1‑149 (MQ=255)                                                                                                                                 gcgTTCGGCGTATACCGCTAGGACCACCACACTATCCTGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGAAGTAGCGAAAGATTCCAGACTCCAGGTGCTGGGGGTATATCGTCCTTTCCCTGc               >  2:172683/1‑149 (MQ=14)                                                                                                                                    ttCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGtt            <  2:376513/149‑1 (MQ=255)                                                                                                                                       gtcGTAGACCGCGAGCAGTCCGGCGATAGCGGGAACGTCCCTGCCGGTGTATTGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATAGCGTCATTGCCCTGCGTTTAc         <  2:476454/147‑1 (MQ=17)                                                                                                                                        gcctataCCGCTAGGACTCCGGCGCTATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGTTATCCGCCGGTCCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACt        <  1:299214/146‑1 (MQ=14)                                                                                                                                        gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAg                                                                   >  2:103495/1‑90 (MQ=255)                                                                                                                                        gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAg                                                                   <  1:103495/90‑1 (MQ=255)                                                                                                                                          gTATACCGCTAGGTCTCCGGCGATATCCGGAACGTCCCTGCTGGTGTAATGACAGACGCTACCGTGTCACTGGTATTTGCCGGTCCAGTAGCGTAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGa      <  2:314257/149‑1 (MQ=17)                                                                                                                                              aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAt  >  1:443487/1‑149 (MQ=255)                                                                                                                                              aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGAGTCGGGACTCAGGATGCTGGCGGTAGATCGTAATTCCCCTTCGTTTACTCATagt  >  1:473156/1‑147 (MQ=17)                                                                                                                                                   |                                                                                                                                                CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAT  >  NC_000913_3_pae_tpiA/263952‑264241

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