Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	512,167	G→T	12.1%	A43A (GCG→GCT)	fetB →	iron exporter permease subunit, ABC transporter FetAB; peroxide resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	512,167	0	G	T	12.1%	99.3 / 2.4	33	A43A (GCG→GCT)	fetB	iron exporter permease subunit, ABC transporter FetAB; peroxide resistance protein
Reads supporting (aligned to +/- strand):  ref base G (16/13);  new base T (2/2);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

AAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGA  >  CP009273/511967‑512423                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                 aaGACGAAATTAATCATGCGGAGAAAGTGATTACACTNCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAAGATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAATACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTAAGGCGTGAGATAATTCAGCTGATTATTGt                                                                                                                                                                                                                                              <  1:353194/221‑1 (MQ=255)         aTTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGTGCGAGCGATAATTCAg                                                                                                                                                                                                                                                         >  2:1910/1‑202 (MQ=255)             aTCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGa                                                                                                                                                                                                                                                     <  1:1910/202‑1 (MQ=255)                          aGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGCATGATGCCAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAAt                                                                                                                                                               <  1:463255/275‑1 (MQ=255)                                             cATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGATGgt                                                                                                                                                                                  >  1:408526/1‑237 (MQ=255)                                                                     cGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCTATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGATGgtg                                                                                                                                                                                 >  2:289430/1‑214 (MQ=255)                                                                         aTGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCTATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGATGGTGttat                                                                                                                                                                             <  1:289438/214‑1 (MQ=255)                                                                                cGCATAATATTACTAACGAATCATTAGCGCTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTCTCAGCGTGGATGATGCCANCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCa                                                                                                         >  1:363905/1‑275 (MQ=255)                                                                                 gCATAATATTACTAACGAATCATTAGCACTGTCATTAATGCTGGTGGTGGGGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATGCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGGTGATGCCAGCCTGa                                                                                                                                                                                               <  1:141830/188‑1 (MQ=255)                                                                                  cATAATATTACTAACGAATCATTAGCACTGGCATTGATGCTGGTGGAGGTGGCAATCTTAATTAGACATAAAGAAAAACTGGCGCTGGAGAAAGATATGATCTGGAGCGTCGGGCGAGAGATAATTAAGCTGAGTAATGTCGGCTATGTGCTGACGTAGATTTTCAGCGTGGATGATGATAGCCTGACATTATTTATGGAGGTATTTAACTGCTTTAATGCGGCGTGGAACGGGGAAACAGGCAGGAAAAATAATTCTAAAgggtacctctcacc                                                                                                       >  2:76954/1‑262 (MQ=255)                                                                                  cATAATATTACTAACGAATCATTAGCACTGGCATCAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCNTGACATTANTGATGGTGTTATTTATCTGCTTTAATGCGGCGTGGAACGCGCAAAAACg                                                                                                                                        >  1:327682/1‑242 (MQ=255)                                                                                                   aTCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGATGAAGTATATTTTCAGCGTGGATGATGCCAGCCt                                                                                                                                                                                                 >  1:6722/1‑168 (MQ=255)                                                                                                       ttAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACa                                                                                                                                                                                             <  2:6722/168‑1 (MQ=255)                                                                                                           cACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTc                                                                                                                                                                                                                                             >  2:205249/1‑116 (MQ=255)                                                                                                               gtacttAATGTTGTTGGGGGTGTAAATTTTAATTAGCCATAATGAACAACTGGCGCTGGAGAAAGTTATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAATTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGaa                                                                          <  2:140046/271‑1 (MQ=255)                                                                                                               ggCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCt                                                                                                                                                                                                                                         <  1:205255/116‑1 (MQ=255)                                                                                                                            gtggtggtggCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCg                                                                              >  1:140052/1‑258 (MQ=255)                                                                                                                                 ggtggCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGATAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTTTGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGGCATTATTGATGGTGTTATTTATCTGCTTTAATGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTAATCGTTTATTGCCATTACGTTCGGGGCGGGAATTACCCTGGAgggnttgn                                                        >  2:562/1‑269 (MQ=255)                                                                                                                                    ggCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATgc                                                                                                                                                             >  2:53597/1‑171 (MQ=255)                                                                                                                                        aTCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATgcggcg                                                                                                                                                         <  1:53599/171‑1 (MQ=255)                                                                                                                                                   ccATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCTATAATTCAGCTGATTATTGTc                                                                                                                                                                                                                                             >  1:171775/1‑76 (MQ=255)                                                                                                                                                      tAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCGGAAGTATATTTTCAGCGTGGATGATGCAAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATGCAGCATGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGTTTTTCTCTCAACTATTAATGCCATTAACGTCCGAGGAGGAATTAAGCTGGCGTTGCTGAATTTCTCacactcgttgtaat                                   >  2:96107/1‑262 (MQ=255)                                                                                                                                                       aaaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCTATAATTCAGCTGATTATTGTCGGCt                                                                                                                                                                                                                                         <  2:171769/76‑1 (MQ=255)                                                                                                                                                           aaaaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGa                                                                                                                                                                                      <  1:47836/123‑1 (MQ=255)                                                                                                                                                             aaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGatta                                                                                                                                                                                                                                                  >  2:161293/1‑61 (MQ=255)                                                                                                                                                                 tGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGt                                                                                                                                                                                                                                              <  1:161299/61‑1 (MQ=255)                                                                                                                                                                    cgcTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATTCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCCCATCGTTTATTGCCATTACGGTAGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCCNNGTCGATTTATTNNTTccccttttt                      >  2:231534/1‑268 (MQ=255)                                                                                                                                                                                   ttCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATACCAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCgg                                                                                  >  1:396527/1‑199 (MQ=255)                                                                                                                                                                                     ctctGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat                                                                                                                                                                                                            >  1:124264/1‑75 (MQ=255)                                                                                                                                                                                       ctGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATACCAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCg                                                                              <  2:396514/199‑1 (MQ=255)                                                                                                                                                                                        tGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGga  <  1:96109/274‑1 (MQ=255)                                                                                                                                                                                         ggAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGatgatg                                                                                                                                                                                                        <  2:124258/75‑1 (MQ=255)                                                                                                                                                                                          gAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTTCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGc                                                           >  1:172131/1‑215 (MQ=255)                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                 AAGACGAAATTAATCATGCGGATAAAGTGATTACACTGCAACCGCATGCCGGAGAAATGCAGGAAGCACGCTATGAACTCGCATAATATTACTAACGAATCATTAGCACTGGCATTAATGCTGGTGGTGGTGGCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCTGACATTATTGATGGTGTTATTTATCTGCTTTAATGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGA  >  CP009273/511967‑512423

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