Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,465,424	T→G	21.6%	A41A (GCA→GCC)	fcuA_1 ←	Ferrichrome receptor FcuA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,465,424	0	T	G	21.6%	82.2 / 4.0	48	A41A (GCA→GCC)	fcuA_1	Ferrichrome receptor FcuA
Reads supporting (aligned to +/- strand):  ref base T (11/25);  new base G (8/4);  total (19/29)
Fisher's exact test for biased strand distribution p-value = 4.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAAAGCCTTGCGCGATACTGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGCGTGGGTCC  >  ACB122/3465292‑3465556                                                                                                                                        |                                                                                                                                     ggCAAAGCCTTGCGCGATACTGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGc                                                                                                                                    <  3:149679/135‑1 (MQ=255) ggCAAAGCCTTGCGCGATACTGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGc                                                                                                                                    <  3:147910/135‑1 (MQ=255)   cAAAGCCTTGCGCGATACTGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCgg                                                                                                                                  <  5:69801/135‑1 (MQ=255)               cgATACTGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTG‑GGGGGCCAAGCGCGCCACTACCA‑CGGCCCCCCGGCC‑GGCCGCCCGGTGGGaac                                                                                                                      >  1:38984/1‑134 (MQ=255)                 aTACTGTAGCTGCCCTCCAGGCCTTGCGTTTGTACGCCCCTGCGCAGGGCGGGGTCTAATGACAACGCCACCCCGGGGTT‑GGCGGCAACGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTgc                                                                                                                   <  7:369174/136‑1 (MQ=255)                     tGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTTGGGCGGC‑AAGCCGGCCAGTACCA‑CGGCCCACCGGCC‑TGCGGCCCCGGGGTTAGTgggcgg                                                                                                               >  8:192204/1‑135 (MQ=255)                        aGCTGCCCGCCCGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTgcgcgc                                                                                                               <  2:88950/133‑1 (MQ=255)                          cTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTG‑GGGGGGCAACCGCGCCAGGACAT‑AGGCCAACGGGGC‑GGCGGCCAGGGGGGTAGTGGGCGCgggg                                                                                                           >  3:364662/1‑133 (MQ=255)                              ccccccggctttggctttGAACCCCCCGGCGCAGGGCCGGGTCGAATCACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTc                                                                                                      <  7:91946/130‑1 (MQ=255)                               cTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTT‑GGCGGGCACGCCCCCCCGGTCCG‑AGGCCAACCGGCC‑GGCGGCCCGGTGGGAAGGGGGCGCGGCGGGGTCg                                                                                                     >  8:206835/1‑136 (MQ=255)                               cTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAACCGCGCCCCGCCATCCGCCCCACGGCCG‑GGGGGCCCGGTGG‑AACTGCGCGCGGGGGGGTCg                                                                                                     >  8:118463/1‑136 (MQ=255)                                       cTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAACCGCGCCAGTTACT‑CGGCCCACCGGCC‑GGCGGCCCCGGGGTAAGGGCGCGGGGCGGTGTGGGCTGccca                                                                                             >  7:237003/1‑134 (MQ=255)                                           gCTTTGTACGCCCCGGAGCCGGGCCGGGGCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATcgc                                                                                          <  3:47807/135‑1 (MQ=255)                                                  aCGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑GGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTggg                                                                                    <  1:23735/134‑1 (MQ=255)                                                   cGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTT‑GGGGGGCAACCCCGCCCGGTCAA‑AGGCCCACCGGCCGGGGGGCCCGGGGGT‑AGTGGGCGGGGGGGGGTGGGCGGGCCTCGCgggggggg                                                                                  >  2:221006/1‑135 (MQ=255)                                                                ggCCGGGTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAACGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCAt                                                                     <  3:361277/135‑1 (MQ=255)                                                                gcccGGGTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCAt                                                                     <  6:213878/133‑1 (MQ=255)                                                                 gCCGGGCCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATg                                                                    <  1:49708/135‑1 (MQ=255)                                                                   cGGGTCGAATGACAACGCCACCCCGGCGTT‑TGGGGGAAAGCGCGGCAGTTCAT‑AGGGCAACGGGCC‑TTCGGGCAGGTGGTTAGTGCGCGCGGAGGTGTGGGCTTCCCTCGCGGGGGGGGGCAGGGCGGCCAttggg                                                                 >  8:113695/1‑132 (MQ=255)                                                                   cGGGTCGAATGACAACGCCACCCCGGCGTT‑GGGGGCCAAGCGCGCCCCTACAA‑AGGCCCAAGGGGC‑GGCGGCCAGGGGGTTAGTGGGCGCGGCGGTGTGGGCTGCCCTCGCCGGGGGGGGCCGGGCGGCCATGGCg                                                                 >  8:316301/1‑136 (MQ=255)                                                                     ggTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGAt                                                               <  7:113695/136‑1 (MQ=255)                                                                     ggTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGGCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGAt                                                               <  7:242570/136‑1 (MQ=255)                                                                      gccGAATCCCAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGAt                                                               <  3:273915/133‑1 (MQ=255)                                                                      gTCGAATGCCAACGCCACGCCGGCGTT‑GGCGGCAAAGCGCGCCAGTAGAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGAt                                                               <  5:126334/135‑1 (MQ=255)                                                                      gTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGGCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATg                                                              <  8:148016/136‑1 (MQ=255)                                                                      gTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑GGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATg                                                              <  8:284482/136‑1 (MQ=255)                                                                      gTCGAAGGCCAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGAt                                                               <  2:386940/135‑1 (MQ=255)                                                                          aaTGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTa                                                          <  8:265069/136‑1 (MQ=255)                                                                             gACAACGCCACCCCGGCGTT‑GGGGGCCAAGCGCGCCCGGACAT‑CGGGCCACGGGCCGGGGGGCCGGGGGGTATG‑GGGCGCGGGGGGGTGGGGTGCCCGAGGGGTGGGGGGCCGGGGGGGCCTGGCGCTggggtgg                                                        >  4:123025/1‑130 (MQ=255)                                                                                aaCGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCtt                                                     <  6:291543/135‑1 (MQ=255)                                                                                                             gcgcCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCt                         <  4:200857/135‑1 (MQ=255)                                                                                                                       aT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCg                >  3:259902/1‑134 (MQ=255)                                                                                                                           ggCCAACGGGCC‑GGCGGCCCGGAGGTATGTGCGCGGGGAGGTGTCGGCTGCCAGCGCGGTGGGGGCCAGGGCGGCGATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGAGTCAt            <  1:177029/135‑1 (MQ=255)                                                                                                                            gCCAACGGGCC‑TGCGGCCAGGTTGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCAt            <  5:187534/134‑1 (MQ=255)                                                                                                                            gCCAACGGGCC‑GGCGGCCAGGTGGGAAGTGCGGGCGGAGGGGTCGGCTGCCATCGCGGTGGGGGCCTGGGCGGCGAGGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATg           <  6:20587/135‑1 (MQ=255)                                                                                                                             ccAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATcgc                                                                                          >  7:243370/1‑55 (MQ=255)                                                                                                                             ccAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATcgc                                                                                          <  8:243370/55‑1 (MQ=255)                                                                                                                             ccAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCGATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGc                      <  5:391562/123‑1 (MQ=255)                                                                                                                             ccAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGc                      >  6:391562/1‑123 (MQ=255)                                                                                                                             ccAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGc          >  6:164983/1‑135 (MQ=255)                                                                                                                             ccAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGc          >  3:113367/1‑135 (MQ=255)                                                                                                                               aaCGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCCAGCGGTACGGCTGGCGAGGCGTCATGCGt        >  4:312001/1‑135 (MQ=255)                                                                                                                               aaCGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCCAGCGGTACGGCTGGCGAGGCGTCATGCGt        >  1:377453/1‑135 (MQ=255)                                                                                                                               aaCGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTc              >  6:282436/1‑129 (MQ=255)                                                                                                                               aaCGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTc              <  5:282436/129‑1 (MQ=255)                                                                                                                                  gggCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCAGCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGAGTCATGCGTggg     >  4:96220/1‑135 (MQ=255)                                                                                                                                     cc‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGCGTGGGTcc  <  4:259902/135‑1 (MQ=255)                                                                                                                                     cc‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGCGTGGGTcc  <  6:317832/135‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     GGCAAAGCCTTGCGCGATACTGTAGCTGCCCTCCAGGCCTTGGCTTTGTACGCCCCTGAGCAGGGCCGGGTCGAATGACAACGCCACCCCGGCGTT‑GGCGGCAAAGCGCGCCAGTACAT‑CGGCCAACGGGCC‑TGCGGCCAGGTGGTAAGTGCGCGCGGAGGTGTCGGCTGCCATCGCGGTGGGGGCCAGGGCGGCCATGGCGATGGCTAGGCTTAGGGCCTTCAAAGGCAAGCGGTACGGCTGGCGAGGCGTCATGCGTGGGTCC  >  ACB122/3465292‑3465556

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