Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,460,220	Δ1 bp	coding (131/195 nt)	bfd ←	bacterioferritin‑associated ferredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,460,218	0	C	.	92.1%	143.0 / 8.3	39	coding (133/195 nt)	bfd	bacterioferritin‑associated ferredoxin
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (17/18);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACA  >  NZ_CP009273/3460080‑3460349                                                                                                                                           |                                                                                                                                    ttCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGc                                                                                                                                     >  2:1219377/1‑139 (MQ=255)  tCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGTCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGc                                                                                                                                     <  1:1137371/138‑1 (MQ=255)       ttttGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGccggac                                                                                                                                <  2:981763/138‑5 (MQ=255)       ttttGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGccggac                                                                                                                                <  1:717233/138‑5 (MQ=255)                cTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCg                                                                                                                     >  2:86967/1‑139 (MQ=255)                          aCTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGc                                                                                                           >  2:28289/1‑139 (MQ=255)                                tGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGAc                                                                                                       <  2:645094/137‑1 (MQ=255)                                        aGAACGGCTTACGAGTAAGTCAAAGAACACACTCCGTTTATGCGGACTCCTTAAGCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTAACCGCACTGATTGCCGACCGGAATaact                                                                                             >  2:25723/1‑137 (MQ=255)                                               cTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAttttttt                                                                                       <  1:679193/139‑1 (MQ=255)                                                 tACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAAt                                                                                    >  1:1111921/1‑139 (MQ=255)                                                  aCGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTaa                                                                                     <  2:1051929/137‑1 (MQ=255)                                                         aGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGaa                                                                            <  1:28289/139‑1 (MQ=255)                                                         aGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGaa                                                                            >  2:1066913/1‑139 (MQ=255)                                                                     cacTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAgggg                                                                 >  2:755593/1‑138 (MQ=255)                                                                      acTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGaa                                                               <  2:1073862/139‑1 (MQ=255)                                                                      acTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGaa                                                               <  2:388403/139‑1 (MQ=255)                                                                       cTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGaa                                                               <  2:1190628/139‑1 (MQ=255)                                                                        tCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAttttttt                                                                                       >  2:338145/1‑113 (MQ=255)                                                                        tCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAttttttt                                                                                       <  1:338145/113‑1 (MQ=255)                                                                                  cGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTaa                                                                                     >  2:302171/1‑105 (MQ=255)                                                                                  cGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTaa                                                                                     <  1:302171/105‑1 (MQ=255)                                                                                  cGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGt                                                                       >  1:1181178/1‑119 (MQ=255)                                                                                  cGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGt                                                                       <  2:1181178/119‑1 (MQ=255)                                                                                    gACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGc                                                 >  1:1263423/1‑139 (MQ=255)                                                                                    gACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGc                                                 <  1:352142/139‑1 (MQ=255)                                                                                      cTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCt                                               >  1:1017106/1‑139 (MQ=255)                                                                                                   ggCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAtttttt                                   >  1:583952/1‑138 (MQ=255)                                                                                                    gCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTa                                 <  2:1253391/139‑1 (MQ=255)                                                                                                           cATTAATTCATCCTCCATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGa                          <  1:486502/139‑1 (MQ=255)                                                                                                            aTTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGa                          <  1:1005174/139‑1 (MQ=255)                                                                                                                      cctccATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACa                >  2:900039/1‑138 (MQ=255)                                                                                                                          cATCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGa            <  1:755593/138‑1 (MQ=255)                                                                                                                           aTCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACa          <  2:862046/139‑1 (MQ=255)                                                                                                                           aTCACTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACa          <  1:1066913/139‑1 (MQ=255)                                                                                                                           aTCACTTCACGCGC‑GTCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACa          <  2:735246/139‑1 (MQ=255)                                                                                                                               cTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACg      <  2:1153606/139‑1 (MQ=255)                                                                                                                               cTTCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACg      <  2:899177/139‑1 (MQ=255)                                                                                                                                ttCACGCGCGG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACg      >  2:100059/1‑138 (MQ=255)                                                                                                                                    cgcgcgG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAGCTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACa  >  2:850922/1‑138 (MQ=255)                                                                                                                                    cgcgcgG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACa          >  2:1096948/1‑130 (MQ=255)                                                                                                                                    cgcgcgG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACa          <  1:1096948/130‑1 (MQ=255)                                                                                                                                      cgcgG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTAc   <  1:727055/135‑1 (MQ=255)                                                                                                                                      cgcgG‑CCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTAc   >  2:727055/1‑135 (MQ=255)                                                                                                                                           |                                                                                                                                    TTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAAGAGTAGAACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTTTTTATCGCTGATACCATTACAAAGACAAACGTACA  >  NZ_CP009273/3460080‑3460349

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