Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	340,444	(A)5→4	coding (558/816 nt)	BW25113_RS01680 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	340,440	0	A	.	98.8%	363.7 / ‑1.7	82	coding (554/816 nt)	BW25113_RS01680	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (36/45);  minor base G (1/0);  total (37/45)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

GTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTT  >  NZ_CP009273/340352‑340524                                                                                         |                                                                                     gTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGaa                                                                                     >  1:455165/1‑90 (MQ=255) gTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGaa                                                                                     >  1:2166030/1‑90 (MQ=255)     gAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGaaaatg                                                                                 >  1:2720792/1‑88 (MQ=255)         aGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCtt                                                                            <  1:2985563/90‑1 (MQ=255)         aGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCtt                                                                            <  2:283410/90‑1 (MQ=255)           aTTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGTTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTcg                                                                          >  2:743653/1‑90 (MQ=255)           aTTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTcg                                                                          >  2:2511048/1‑90 (MQ=255)           aTTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTcg                                                                          >  2:391611/1‑90 (MQ=255)           aTTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTcg                                                                          >  1:1338643/1‑90 (MQ=255)           aTTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTcg                                                                          >  1:2419543/1‑90 (MQ=255)           aTTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTcg                                                                          >  1:69294/1‑90 (MQ=255)             tGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTcgcg                                                                        <  2:1200991/90‑1 (MQ=255)              gcgaaTGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGa                                                                       <  2:2283259/87‑1 (MQ=255)              gCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGa                                                                       <  2:2245327/90‑1 (MQ=255)                aaaTGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGaaa                                                                     <  2:1957780/90‑1 (MQ=255)                aaaTGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGaaa                                                                     <  1:2741824/90‑1 (MQ=255)                 aaTGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTG‑‑AAAATGGCTTCGCGAAATg                                                                   <  1:2300508/90‑1 (MQ=255)                  aTGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATg                                                                   <  1:1774397/90‑1 (MQ=255)                    gTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtc                                                                 >  1:2030361/1‑90 (MQ=255)                     tATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtcg                                                                <  1:11795/90‑1 (MQ=255)                      aTTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtcgt                                                               >  2:1836071/1‑90 (MQ=255)                      aTTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtcgt                                                               >  1:2748866/1‑90 (MQ=255)                         cttcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTa                                                            >  2:2170492/1‑90 (MQ=255)                         cttcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTa                                                            >  2:33082/1‑90 (MQ=255)                         cttcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTa                                                            <  2:2419543/90‑1 (MQ=255)                         cttcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTa                                                            <  2:1338643/90‑1 (MQ=255)                              tGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCg                                                       <  2:2537119/90‑1 (MQ=255)                              tGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCg                                                       >  1:1215911/1‑90 (MQ=255)                                   aCCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGagac                                                  >  1:1811809/1‑90 (MQ=255)                                   aCCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATG‑GGGAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagac                                                  >  2:1090949/1‑90 (MQ=255)                                      atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                               >  1:3008866/1‑90 (MQ=255)                                      atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                               >  2:2730227/1‑90 (MQ=255)                                      atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                               >  1:1073618/1‑90 (MQ=255)                                      atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                               >  2:920149/1‑90 (MQ=255)                                       tatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacagac                                              <  2:2497321/90‑1 (MQ=255)                                         tGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgt                                                                  <  2:2080603/68‑1 (MQ=255)                                         tGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgt                                                                  >  1:2080603/1‑68 (MQ=255)                                         tGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgc                                            <  2:1836197/90‑1 (MQ=255)                                          gCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                           <  1:1384084/90‑1 (MQ=255)                                           ccTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACgcgc                                          <  2:1666606/90‑1 (MQ=255)                                            cTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGcc                                         >  2:1211695/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgc                                            <  1:1452114/85‑1 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgc                                            >  2:1452114/1‑85 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                       >  2:2599740/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                       >  2:2706414/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                       >  1:1048955/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                       >  1:2889249/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                       >  2:1065671/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                       >  1:1221503/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                       >  1:1302612/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                       >  1:679544/1‑90 (MQ=255)                                               tGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                                      <  1:2511048/90‑1 (MQ=255)                                               tGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                                      <  2:1073618/90‑1 (MQ=255)                                               tGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                                      <  1:2607367/90‑1 (MQ=255)                                                     gaATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCa                                >  1:2178913/1‑90 (MQ=255)                                                     gaATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCa                                >  1:232519/1‑90 (MQ=255)                                                      aaTTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAg                               <  1:743653/90‑1 (MQ=255)                                                      aaTTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAg                               <  1:939852/90‑1 (MQ=255)                                                      aaTTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAg                               <  1:2308547/90‑1 (MQ=255)                                                      aaTTTTCCCTAAAGGGAAAGAGTGTAGGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAg                               <  1:1353743/90‑1 (MQ=255)                                                           tCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCt                          <  2:273512/90‑1 (MQ=255)                                                                 aaGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa                    <  2:50732/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                   <  1:2600980/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                   <  1:1472997/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                   <  2:2720792/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                   <  1:2735913/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                   <  1:1028707/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                   <  1:2828315/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                   <  2:1164742/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                   <  1:490542/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                   <  2:1236575/90‑1 (MQ=255)                                                                     gAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGa                <  2:1221503/90‑1 (MQ=255)                                                                     gAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGa                <  2:455165/90‑1 (MQ=255)                                                                      aaaGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGat               >  1:900914/1‑90 (MQ=255)                                                                       aaGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt              <  2:1004970/90‑1 (MQ=255)                                                                               aTGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTac      <  2:2461083/90‑1 (MQ=255)                                                                               aTGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTac      <  2:2178913/90‑1 (MQ=255)                                                                                  cATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACt   <  1:2742703/90‑1 (MQ=255)                                                                                   aTGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa                    >  2:623781/1‑72 (MQ=255)                                                                                   aTGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCa                    <  1:623781/72‑1 (MQ=255)                                                                                   aTGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt              <  1:2677804/78‑1 (MQ=255)                                                                                   aTGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt              >  2:2677804/1‑78 (MQ=255)                                                                                   aTGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACtt  >  2:1938787/1‑90 (MQ=255)                                                                                   aTGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACtt  >  2:2036903/1‑90 (MQ=255)                                                                                   aTGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACtt  >  2:2074901/1‑90 (MQ=255)                                                                                         |                                                                                     GTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTT  >  NZ_CP009273/340352‑340524

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