Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	21,250	T→C	S24P (TCC→CCC)	yaaY →	DUF2575 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	21,250	0	T	C	75.4%	90.3 / 41.4	57	S24P (TCC→CCC)	yaaY	DUF2575 family protein
Reads supporting (aligned to +/- strand):  ref base T (7/7);  new base C (20/23);  total (27/30)
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.02e-01

GCCGAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACG  >  NZ_CP009273/21119‑21383                                                                                                                                    |                                                                                                                                      gCCGAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGAt                                                                                                                                <  2:579906/139‑1 (MQ=255)    gAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATaa                                                                                                                              >  1:625230/1‑138 (MQ=255)    gAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATaa                                                                                                                              >  2:837314/1‑138 (MQ=255)    gAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATaa                                                                                                                              >  1:776341/1‑138 (MQ=255)      ggAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATaaaaa                                                                                                                           <  2:67617/139‑1 (MQ=255)      ggAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATaaaaa                                                                                                                           <  2:200833/139‑1 (MQ=255)                  tttAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCg                                                                                                               >  1:117121/1‑139 (MQ=255)                  tttAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCg                                                                                                               >  1:359902/1‑139 (MQ=255)                     aGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGt                                                                                                             <  2:359902/138‑1 (MQ=255)                           cttTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAgg                                                                                                                 <  1:129894/128‑1 (MQ=255)                           cttTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAgg                                                                                                                 >  2:129894/1‑128 (MQ=255)                              ttGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTttc                                                                                                      <  1:837314/136‑1 (MQ=255)                                       gATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGAtt                                                                                           <  2:117121/138‑1 (MQ=255)                                             tgCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGcc                                                                                     <  1:481632/138‑1 (MQ=255)                                                  aTAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCa                                                                               <  2:502618/139‑1 (MQ=255)                                                      gcgcCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGa                                                                           <  1:783613/139‑1 (MQ=255)                                                      gcgcCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATg                                                                            >  1:676418/1‑138 (MQ=255)                                                       cgcCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGa                                                                           >  1:282345/1‑138 (MQ=255)                                                       cgcCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGa                                                                           >  1:323509/1‑138 (MQ=255)                                                            ttAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCCGCCTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGc                                                                      >  1:339804/1‑138 (MQ=255)                                                                    ccGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCGCGTGATTTGCGCCATGGCAATGAAAAGCCACttcttt                                                             >  2:273777/1‑139 (MQ=255)                                                                    ccGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACttcttt                                                             >  2:342096/1‑139 (MQ=255)                                                                    ccGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACttcttt                                                             >  1:791528/1‑139 (MQ=255)                                                                        cACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGa                                                         <  1:765887/139‑1 (MQ=255)                                                                                     gTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAg                                                                       <  1:160731/112‑1 (MQ=255)                                                                                     gTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAg                                                                       >  2:160731/1‑112 (MQ=255)                                                                                      tGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAAt                                           <  2:791528/139‑1 (MQ=255)                                                                                       gATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACttctttct                                                           <  2:473172/122‑1 (MQ=255)                                                                                       gATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACttctttct                                                           >  1:473172/1‑122 (MQ=255)                                                                                       gATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATc                                          >  1:315083/1‑139 (MQ=255)                                                                                        aTGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCg                                         >  1:172903/1‑139 (MQ=255)                                                                                        aTGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAgg                                                                                                                 >  1:177355/1‑67 (MQ=255)                                                                                        aTGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAgg                                                                                                                 <  2:177355/67‑1 (MQ=255)                                                                                        aTGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATc                                          <  2:653817/138‑1 (MQ=255)                                                                                             gcAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGt                                    <  1:360671/139‑1 (MQ=255)                                                                                              cAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGt                                    >  2:780099/1‑138 (MQ=255)                                                                                              cAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGt                                    >  2:360671/1‑138 (MQ=255)                                                                                               aGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTa                                  <  1:768237/139‑1 (MQ=255)                                                                                                ggATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTaa                                 >  1:254209/1‑139 (MQ=255)                                                                                                ggATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTaa                                 >  1:771933/1‑139 (MQ=255)                                                                                                        ccAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGAc                          <  1:97205/138‑1 (MQ=255)                                                                                                         cAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCg                        <  2:68311/139‑1 (MQ=255)                                                                                                         cAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCg                        <  2:282345/139‑1 (MQ=255)                                                                                                         cAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCg                        <  1:294468/139‑1 (MQ=255)                                                                                                         cAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCg                        <  2:339804/139‑1 (MQ=255)                                                                                                         cAGCTTGCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCg                        <  2:822500/139‑1 (MQ=255)                                                                                                               gCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTa                   <  1:342096/138‑1 (MQ=255)                                                                                                               gCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTa                   <  1:695093/138‑1 (MQ=255)                                                                                                               gCGGGGTGTGAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTa                   <  1:742711/138‑1 (MQ=255)                                                                                                                cGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTAc                  >  2:546703/1‑138 (MQ=255)                                                                                                                     tgtgAATACAGCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGta            >  2:26350/1‑139 (MQ=255)                                                                                                                             cAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTccc                                                                                                 <  2:485599/46‑1 (MQ=255)                                                                                                                             cAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTccc                                                                                                 >  1:485599/1‑46 (MQ=255)                                                                                                                               gCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGAcc                         >  2:368684/1‑116 (MQ=255)                                                                                                                               gCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGAcc                         <  1:368684/116‑1 (MQ=255)                                                                                                                               gCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGAc   >  1:407019/1‑138 (MQ=255)                                                                                                                               gCTTTCCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACg  <  2:315083/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      GCCGAGGAATATGCCCTTTTAGCCTTCTTTTGTCAATGGATTTGTGCAAATAAGCGCCGTTAATGTGCCGGCACTCGTTACGTAGTGATGGCGCAGGATTCTACCAGCTTGCGGGGTGTGAATACAGCTTTTCCGCGATAAAAATTGCAGCAGGCGGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACG  >  NZ_CP009273/21119‑21383

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