Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	425,348	A→C	11.2%	intergenic (‑39/+18)	PP_0349 ← / ← PP_0350	membrane protein/ferrichrome‑iron receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	425,348	0	A	C	11.2%	133.8 / 8.8	53	intergenic (‑39/+18)	PP_0349/PP_0350	membrane protein/ferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base C (3/3);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

AAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTC  >  NC_002947/425208‑425494                                                                                                                                             |                                                                                                                                                     aaGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGcc                                                                                                                                              >  5:2726/1‑149 (MQ=255)  aGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGcc                                                                                                                                              <  6:292082/148‑1 (MQ=255)    gggCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCa                                                                                                                                           <  6:290808/149‑1 (MQ=255)    gggCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCa                                                                                                                                           <  8:20735/149‑1 (MQ=255)     ggccggccGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCGCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGAGTGTTGTCGGTTGGAACCCGGGGCCGCCTTGCGGccccgg                                                                                                                                          <  6:83895/149‑3 (MQ=255)             ggTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGctcactca                                                                                                                                  <  8:21693/149‑1 (MQ=255)             ggTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGctcactca                                                                                                                                  <  6:2726/149‑1 (MQ=255)                gACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGaa                                                                                                                               >  3:158614/1‑149 (MQ=255)                   tGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACtt                                                                                                                            >  1:81876/1‑149 (MQ=255)                    gCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTg                                                                                                                           >  5:229744/1‑149 (MQ=255)                     ccAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACtt                                                                                                                            <  6:155290/147‑1 (MQ=255)                      cAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTa                                                                                                                         <  6:220147/149‑1 (MQ=255)                                  gAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAggcggc                                                                                                                                               >  2:52239/1‑115 (MQ=255)                                  gAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAggcggc                                                                                                                                               <  1:52239/115‑1 (MQ=255)                                        ggTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAAc                                                                                                     >  5:149519/1‑149 (MQ=255)                                        ggTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAAc                                                                                                     >  5:158908/1‑149 (MQ=255)                                              gACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTgcg                                                                                               <  3:135070/149‑1 (MQ=255)                                              gACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGACGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTAGCTTGGAACCCGGGGCCGCAAGGCGGACACAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTgcg                                                                                               >  1:53410/1‑149 (MQ=255)                                                    gctgctCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTgcg                                                                                               <  1:171713/143‑1 (MQ=255)                                                            aGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCgtagt                                                                                 <  4:158614/149‑1 (MQ=255)                                                             ggTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCgtagta                                                                                <  2:81876/149‑1 (MQ=255)                                                               tgtgGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGtagtag                                                                               >  1:59417/1‑148 (MQ=255)                                                                   gATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagta                                                                          >  3:86492/1‑149 (MQ=255)                                                                      cAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCCGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtag                                                                         >  3:305491/1‑147 (MQ=255)                                                                       aaGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtagt                                                                        <  2:65966/147‑1 (MQ=255)                                                                       aaGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCagtagt                                                                        >  1:65966/1‑147 (MQ=255)                                                                       aaGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAgtagta                                                                       <  6:264871/148‑1 (MQ=255)                                                                       aaGGACCAGCGGCGGAGGGGGGGGCTCTTCATCACGGTACCGATTGGTGGCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGtagtag                                                                      <  6:41503/149‑1 (MQ=255)                                                                                    ggATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTag                                                         <  6:158908/149‑1 (MQ=255)                                                                                        ggtggGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagagag                                                     >  7:30787/1‑149 (MQ=255)                                                                                         gtggGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAgg                                                    <  6:149519/149‑1 (MQ=255)                                                                                                     atcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGtt                                         <  3:173550/148‑1 (MQ=255)                                                                                                     atcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCCCAAGGCGGCCCCCGCTCACACACAACATGTAAATTGCGC‑‑TGGGGACAACGCAGCGCGAATAGCCGGATTTCCAATAGTATGAGTCAAAAGTAGAAAAGTAATACtttgt                                         >  4:173550/1‑145 (MQ=255)                                                                                                      tcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGGCCCAGCTCACTCAGAACTTGTAAGTGGCGCCGTGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTagagag                                                     <  1:254213/137‑1 (MQ=255)                                                                                                      tcaCGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCGGCGCTGATGGCCGGAGTAGCAGTAGTAGGAGTCGCAAGTATAGAGGGAATCCTTGtt                                         >  4:93274/1‑147 (MQ=255)                                                                                                               cGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACGAGTAGAGAGGTAATCCTTGTTCAGCAGGTtgg                              >  5:3119/1‑149 (MQ=255)                                                                                                                 acaacaatCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTtggtg                            <  7:303376/142‑1 (MQ=255)                                                                                                                     ttgtCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACGAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCgg                        <  6:3119/149‑1 (MQ=255)                                                                                                                        tCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTGACTAAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATGGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGa                     >  5:259117/1‑149 (MQ=255)                                                                                                                        tCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTg                      <  6:229744/148‑1 (MQ=255)                                                                                                                        tCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGa                     <  4:86492/149‑1 (MQ=255)                                                                                                                               gAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcga               <  2:59417/148‑1 (MQ=255)                                                                                                                                aaCCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTcagc                                     >  7:252677/1‑125 (MQ=255)                                                                                                                                aaCCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTcagc                                     <  8:252677/125‑1 (MQ=255)                                                                                                                                   ctggggCCGCCTTGCGGCCCCGGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGtt                                >  6:259202/3‑127 (MQ=255)                                                                                                                                   ctggggCCGCCTTGCGGCCCCGGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGtt                                <  5:259202/125‑1 (MQ=255)                                                                                                                                   ctggggCCGCCTTGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgac              <  5:66763/143‑1 (MQ=255)                                                                                                                                   ctggggCCGCCTTGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGcgac              >  6:66763/3‑145 (MQ=255)                                                                                                                                     ggggCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTg                                           >  5:6949/1‑114 (MQ=255)                                                                                                                                     ggggCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTg                                           <  6:6949/114‑1 (MQ=255)                                                                                                                                      gggCCGCCTTGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGAcgcg       >  2:47615/1‑149 (MQ=255)                                                                                                                                         ccGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGccc    >  4:64965/1‑149 (MQ=255)                                                                                                                                           gCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTc  >  1:137100/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                     AAGGGGCAGGCCGGTGACTGCCAGCAACAGCAGGAACAGGGTACAGACCAGGCTGCTCCAGGTGTGGATCAAGGACCAGCGGCGGATGGTGGGGCTCTTCATCACGGTACCGATTGTTGTCGGTTGGAACCCGGGGCCGCAAGGCGGCCCCAGCTCACTCAGAACTTGTAAGTGGCGC‑‑TGGCGACAACGCTGCGCTCATCGCCGTAGTAGCAGTAGTAGGAGTCACAAGTAGAGAGGTAATCCTTGTTCAGCAGGTTGGTGGCGGTGACCGCGACCGACGCGCCCTC  >  NC_002947/425208‑425494

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