Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	2,781,639	T→A	6.9%	K112N (AAA→AAT)	U712_14340 ←	Primosomal protein DnaI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,781,639	0	T	A	6.9%	171.8 / 3.6	58	K112N (AAA→AAT)	U712_14340	Primosomal protein DnaI
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base A (2/2);  total (31/27)
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.09e-01

TCTCATTATAGCTTTTCAGGAAGTCTGTAACATGCTGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCTTCCGAACAATAGGAGCAA  >  CP006881/2781503‑2781777                                                                                                                                         |                                                                                                                                           tctgATTATAGCTTTTCAGGAAGTCAGTAACATCCTGAAACATCGCGAGTCTGCTTGCGTCACTGATATCAACTTGCTGGAACGTCGCTCAAAGAAGATCCTGCAGGATATACATGCATTTCATAAGAGACTgttgtta                                                                                                                                          >  6:95783/1‑138 (MQ=255)          aGCTTTTCAGGAAGTCTGTAACATGCTGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTctgctg                                                                                                                                 >  5:170177/1‑139 (MQ=255)             ttttCAGGAAGTCTTTAACATGCTGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTc                                                                                                                              <  7:137825/139‑1 (MQ=255)               ttCAGGAAGTCTGTAACATGCTGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTc                                                                                                                              >  1:177924/1‑137 (MQ=255)               ttCAGGAAGTCTGTAACATGCTGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTc                                                                                                                              >  3:152032/1‑137 (MQ=255)                                  gCTGCAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGATGGAACGTCGCTCCAAGAAGAGCCTGCTGGATATACATGCTTTTCATAAGAGACTGTGGGTTCGTCTGCGGGTCGAGCTTCCGTTTGACTGGACa                                                                                                         >  7:87374/1‑139 (MQ=255)                                  gCTGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACa                                                                                                         >  1:220997/1‑139 (MQ=255)                                  gCTGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACa                                                                                                         >  2:141014/1‑139 (MQ=255)                                  gCTGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACa                                                                                                         >  3:117865/1‑139 (MQ=255)                                  gCTGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACa                                                                                                         >  7:171474/1‑139 (MQ=255)                                  gCTGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCCCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACa                                                                                                         >  8:23572/1‑139 (MQ=255)                                   cTGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCCCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACAt                                                                                                        <  7:23572/139‑1 (MQ=255)                                    tGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACAt                                                                                                        >  1:181858/1‑138 (MQ=255)                                    tGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGATTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACAtt                                                                                                       >  3:210287/1‑139 (MQ=255)                                    tGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGATTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACAtt                                                                                                       >  3:206299/1‑139 (MQ=255)                                      aaaCATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCg                                                                                                     <  1:141014/139‑1 (MQ=255)                                      aaaCATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCGGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCg                                                                                                     <  7:155952/139‑1 (MQ=255)                                        aCATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGATTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTa                                                                                                   <  4:206299/139‑1 (MQ=255)                                        aCATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGATTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTa                                                                                                   <  4:210287/139‑1 (MQ=255)                                             gcgAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCTAGCTTCCGTTTGACTGGACATTCGTAATACt                                                                                              >  4:244330/1‑139 (MQ=255)                                                aGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGa                                                                                           <  4:49707/139‑1 (MQ=255)                                                aGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGa                                                                                           <  6:204324/139‑1 (MQ=255)                                                 gTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGat                                                                                          <  6:170177/139‑1 (MQ=255)                                                 gTCTGCGTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGat                                                                                          <  6:210489/139‑1 (MQ=255)                                                  tCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGata                                                                                         >  5:294416/1‑139 (MQ=255)                                                      cTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCaa                                                                                     <  2:181858/139‑1 (MQ=255)                                                           gTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGa                                                                                           >  4:285789/1‑128 (MQ=255)                                                           gTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGa                                                                                           <  3:285789/128‑1 (MQ=255)                                                           gTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGAc                                                                                <  6:294416/139‑1 (MQ=255)                                                           gTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTATGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGAc                                                                                >  2:202732/1‑139 (MQ=255)                                                              aCTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGAc                                                                                <  5:143683/136‑1 (MQ=255)                                                              aCTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGAc                                                                                >  6:143683/1‑136 (MQ=255)                                                                        aCTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGacaa                                                                   <  4:206461/139‑1 (MQ=255)                                                                            gCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAg                                                                                  >  8:219043/1‑120 (MQ=255)                                                                            gCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAg                                                                                  <  7:219043/120‑1 (MQ=255)                                                                            gCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAg                                                               >  3:117685/1‑139 (MQ=255)                                                                              tGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCt                                                             >  4:25808/1‑139 (MQ=255)                                                                                gAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTc                                                           >  4:283486/1‑139 (MQ=255)                                                                                 aaCGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCg                                                          >  1:235219/1‑139 (MQ=255)                                                                                    gTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGGTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGAt                                                       <  4:152032/139‑1 (MQ=255)                                                                                      cGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGAc                                                                                                          >  1:137908/1‑86 (MQ=255)                                                                                      cGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGAc                                                                                                          <  2:137908/86‑1 (MQ=255)                                                                                       gCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATgg                                                     >  8:224977/1‑138 (MQ=255)                                                                                                   tCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCcaaca                                        <  5:134010/139‑1 (MQ=255)                                                                                                     ctgctgGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACaaa                                      <  2:235219/139‑1 (MQ=255)                                                                                                       gctgGATATACATGCTTTTCATAAGAGACTGTTCTTTCATCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCAGTGACAACAAGCTGCGGATGGTAGCCCGCCAACAAAtt                                    >  6:98841/1‑139 (MQ=255)                                                                                                                aCATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGtt                           >  8:255640/1‑139 (MQ=255)                                                                                                                     cTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGAAATACTCGATATCAATAGCCCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCTACAAAATGTTGCAGGTGTCAt                      >  8:105835/1‑139 (MQ=255)                                                                                                                       tttCATAAGAGACTGTGGTGTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACAGTCGTATTACGCGATATCAATAGACCTCCCATTGACCACAAGCTTCGGATGGTAGCCCTCCCACAAATTGGTGCAGTGTTctgct                    >  2:1567/1‑135 (MQ=255)                                                                                                                        ttCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGCCCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCtt                   >  5:232916/1‑139 (MQ=255)                                                                                                                           aTAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAATCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCTTCCg                <  3:244330/139‑1 (MQ=255)                                                                                                                           aTAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCTTCCg                <  7:224977/139‑1 (MQ=255)                                                                                                                             aaGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCTTCCGaa              >  2:263967/1‑139 (MQ=255)                                                                                                                              agagACTGTTGTTTCTTCTGCGGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCTTCCGaa              <  2:177924/138‑1 (MQ=255)                                                                                                                                 gaCTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCTTCCGAACAAt          <  2:220997/139‑1 (MQ=255)                                                                                                                                 gaCTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCTTCCGAACAAt          <  8:171474/139‑1 (MQ=255)                                                                                                                                 gaCTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCTTCCGAACAAt          <  4:117865/139‑1 (MQ=255)                                                                                                                                  aCTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAAttgttg                                >  2:287938/1‑116 (MQ=255)                                                                                                                                  aCTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAAttgttg                                <  1:287938/116‑1 (MQ=255)                                                                                                                                    tgttgtTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCTTCCGAACAATAg        >  6:149404/1‑138 (MQ=255)                                                                                                                                      ttgttTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCTTCCGAACAATAGGAg     <  3:25808/139‑1 (MQ=255)                                                                                                                                         tttCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCTTCCGAACAATAGGAGCaa  <  7:255640/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                           TCTCATTATAGCTTTTCAGGAAGTCTGTAACATGCTGAAACATCGCGAGTCTGCTTGGGTCACTGATATCAACTTGCTGGAACGTCGCTCCAAGAAGATCCTGCTGGATATACATGCTTTTCATAAGAGACTGTTGTTTCTTCTGCTGGTCGAGCTTCCGTTTGACTGGACATTCGTAATACTCGATATCAATAGACCTCCCATTGACAACAAGCTTCGGATGGTAGCCCTCCAACAAATTGTTGCAGTTTTCATCTTCCGAACAATAGGAGCAA  >  CP006881/2781503‑2781777

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