Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	538,630	C→A	10.8%	A333D (GCC→GAC)	ybbW →	putative allantoin transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	538,630	0	C	A	10.8%	92.5 / 2.8	37	A333D (GCC→GAC)	ybbW	putative allantoin transporter
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base A (2/2);  total (19/18)
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.49e-01

CCGGAGCCAGTATTCACTACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGATATTA  >  NC_000913/538483‑538770                                                                                                                                                    |                                                                                                                                             ccGGAGCCAGTATTCACTACGGCGCTGATACCTGGAACGTGCTCGATATTGTTCAGCGTTGGGACAGCCTGGAGGCATCGGTCTTTGCGGTACAGGATATTATTGTGTGAGTTATCGCAACTCACGCGACCGGCAATCTAATTGCAGcg                                                                                                                                             >  2:511849/1‑148 (MQ=255)   ggggcgcgTATTCACTACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCgg                                                                                                                                           <  2:623183/142‑1 (MQ=255)        cAGTATTCACTACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATc                                                                                                                                      >  2:151316/1‑149 (MQ=255)             ttCACTACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGGTGACAACTATCTCCACTAACGCGGCCGGTAATTTTATTCCAGCCGGTTTTcgggtt                                                                                                                                 >  2:319396/1‑144 (MQ=255)             ttCACTACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGACGGTTATCAGAtt                                                                                                                                 >  1:314748/1‑149 (MQ=255)             ttCACTACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGACGGTTATCAGAtt                                                                                                                                 >  1:395381/1‑149 (MQ=255)                 cTACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTgccg                                                                                                                             >  1:469193/1‑149 (MQ=255)                  tACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTgccgc                                                                                                                            >  1:769091/1‑149 (MQ=255)                    cGGCGCTGATACCTGGAACGGGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCa                                                                                                                          <  1:81859/149‑1 (MQ=255)                     ggCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCAt                                                                                                                         >  1:206316/1‑149 (MQ=255)                               ccTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATGGCACCGCCa                                                                                                               >  2:383580/1‑149 (MQ=255)                                    caCGTGTTGGTTATTGTTCATCTTTGGGCCAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCAAGCGGTTTATCAGATTGCCGCCATTGCACCGACAAAACt                                                                                                          <  2:362741/148‑1 (MQ=255)                                          cTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACACAACTGACCta                                                                                                    >  1:246382/1‑149 (MQ=255)                                          cTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTGCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCCCCGACAAACCTGCCCtt                                                                                                    >  1:81185/1‑148 (MQ=255)                                           tGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCtat                                                                                                   >  1:265921/1‑149 (MQ=255)                                                  tGTCCAGCTTTGGGACGGAACGTACGCGTCTGTCTTGGCGGTACTGGTTATTCTTATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCATCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACg                                                                                            <  2:745881/149‑1 (MQ=255)                                                      cAGCGTTGGGACAGCCTGTTCTCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCCGCCGGTTTGTTGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGt                                                                                        <  1:154253/149‑1 (MQ=255)                                                       cgCGTGGGGACCGCCTGTCCCCCTCGTTCTCTGCGGTACTGGTTATTCTGATGACAACTATCTCCACCAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTa                                                                                       <  2:591495/148‑1 (MQ=255)                                                       aGCGGTGGGAGAGGCTGTTCGCATCGGTCCTTGCGGTAGCGGTTATGCGGATGACAACGCTCTCCACTAAGGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTa                                                                                       <  2:591383/149‑1 (MQ=255)                                                         cGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACt                                                                                     <  2:206316/149‑1 (MQ=255)                                                                             cTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCt                                                                 >  1:58067/1‑149 (MQ=255)                                                                               cGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTg                                                               >  1:9729/1‑149 (MQ=255)                                                                                  tCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTg                                                            <  2:769091/149‑1 (MQ=255)                                                                                     ttGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATc                                                         <  1:723742/149‑1 (MQ=255)                                                                                        cGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGc                                                      >  2:630567/1‑149 (MQ=255)                                                                                        cGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGACGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGc                                                      <  2:314748/149‑1 (MQ=255)                                                                                        cGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGACGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGc                                                      <  2:395381/149‑1 (MQ=255)                                                                                        cGGTACCGGTTATTCCGATGACAACCATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGc                                                      <  1:319396/149‑1 (MQ=255)                                                                                           tACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCg                                                   <  2:469193/149‑1 (MQ=255)                                                                                           tACTGGTTATTCTGATGAAGACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCACTATTATCAGCTTGCTGATCTGCCCg                                                   >  1:574851/1‑149 (MQ=255)                                                                                              tGGTCATTCTGATGACAACTACCTCCACGAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTgg                                                <  1:242297/149‑1 (MQ=255)                                                                                                             cAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTACCGCCATAGCACCTACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGATTGCTGAACTGCCCGTGGGAATTAATGGAAAAt                                 >  1:477260/1‑149 (MQ=255)                                                                                                              aaCTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCACGGGGAACTTAATGGAAAAta                                >  2:605702/1‑148 (MQ=255)                                                                                                                                  ccGGTAATAGTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTc            <  1:627754/149‑1 (MQ=255)                                                                                                                                           ttattCCAGCCGGTTATCAGATTGCCGCCATTGCCCCGACAAAAATGAACCATAAACACGCCGTACTGCTTCCAACAATTCTTAGCCTGCTGACGAGCTCTGGGAAGATAATGCAGGATCGGGTCAGCAGTTAGCTTTTCCTCGAtatt   >  2:379248/1‑149 (MQ=255)                                                                                                                                           ttattCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCCTTTCCTCGAcatt   >  2:404028/1‑149 (MQ=255)                                                                                                                                            tcttCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGAtatta  <  2:575518/147‑1 (MQ=255)                                                                                                                                            tagtCCAGCCGGTTATCAGGTTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGAtatta  <  2:131653/146‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                             CCGGAGCCAGTATTCACTACGGCGCTGATACCTGGAACGTGCTGGATATTGTTCAGCGTTGGGACAGCCTGTTCGCCTCGTTCTTTGCGGTACTGGTTATTCTGATGACAACTATCTCCACTAACGCGACCGGTAATATTATTCCAGCCGGTTATCAGATTGCCGCCATTGCACCGACAAAACTGACCTATAAAAACGGCGTACTGATTGCCAGTATTATCAGCTTGCTGATCTGCCCGTGGAAATTAATGGAAAATCAGGACAGCATTTATCTTTTCCTCGATATTA  >  NC_000913/538483‑538770

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: ‑