Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,126	A→C	11.7%	I419S (ATT→AGT)	trmF ←	tRNA uridine 5‑carboxymethylaminomethyl modification protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,126	0	A	C	11.7%	137.2 / 9.1	50	I419S (ATT→AGT)	trmF	tRNA uridine 5‑carboxymethylaminomethyl modification protein
Reads supporting (aligned to +/- strand):  ref base A (23/21);  new base C (3/3);  total (26/24)
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.85e-01

GAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGCC  >  NC_002947/2981‑3273                                                                                                                                                   |                                                                                                                                                    gAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCa                                                                                                                                                  <  6:25697/149‑1 (MQ=255)   aTGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAgg                                                                                                                                                <  1:293529/149‑1 (MQ=255)    tGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACACTGCGGTACGGCTCCTGAGTACCCAAGGTAATCAgg                                                                                                                                                <  1:54789/149‑1 (MQ=255)    tGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGATTACCTTGGGTACTCAgga                                                                                                                                               <  1:139370/149‑2 (MQ=255)    tGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGgt                                                                                                                                               >  5:124666/1‑149 (MQ=255)             gcgcTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTccacca                                                                                                                                      <  3:128534/149‑1 (MQ=255)                   gTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGc                                                                                                                                >  5:43684/1‑149 (MQ=255)                     cATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCa                                                                                                                                                  >  7:90943/1‑129 (MQ=255)                     cATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCa                                                                                                                                                  <  8:90943/129‑1 (MQ=255)                      aTCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGTTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTcca                                                                                                                                         <  6:2904/137‑1 (MQ=255)                      aTCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTcca                                                                                                                                         >  5:2904/1‑137 (MQ=255)                         gATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGt                                                                                                                          <  5:333585/149‑1 (MQ=255)                              ggCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTAt                                                                                                                        >  5:313045/1‑146 (MQ=255)                              ggCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTAt                                                                                                                        <  6:313045/146‑1 (MQ=255)                              ggCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGcc                                                                                                                     <  8:115920/149‑1 (MQ=255)                                cccAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCtc                                                                                                                   <  3:297399/149‑1 (MQ=255)                                   aGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCtcgtc                                                                                                                >  3:131455/1‑149 (MQ=255)                                        aCGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggc                                                                                                           >  8:148824/1‑149 (MQ=255)                                          ggCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggcgc                                                                                                         >  3:89140/1‑149 (MQ=255)                                           gCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggcgcg                                                                                                        <  3:73480/149‑1 (MQ=255)                                             ccTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggcgcg                                                                                                        >  3:165550/1‑147 (MQ=255)                                             ccTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggcgcg                                                                                                        <  4:165550/147‑1 (MQ=255)                                                    ggTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGc                                                                                               <  6:43684/149‑1 (MQ=255)                                                          gcgcAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCAc                                                                                         >  6:216190/1‑149 (MQ=255)                                                            gcAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACgg                                                                                       <  6:124666/149‑1 (MQ=255)                                                                               gcgcAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCtcgtcg                                                                                                               >  4:95995/1‑106 (MQ=255)                                                                                cacAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTgcggc                                                                   >  5:313633/3‑149 (MQ=255)                                                                                    gTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTa                                                               <  4:131455/149‑1 (MQ=255)                                                                                         aGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTccc                                                          >  1:24443/1‑149 (MQ=255)                                                                                          gCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCg                                                         >  8:48240/1‑149 (MQ=255)                                                                                                    ggCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAgg                                                <  5:216190/148‑1 (MQ=255)                                                                                                    ggCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGc                                               >  6:292801/1‑149 (MQ=255)                                                                                                                   cA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGct                                 >  8:303628/1‑148 (MQ=255)                                                                                                                    a‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGaa                                                                          <  5:315105/106‑1 (MQ=255)                                                                                                                    a‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGaa                                                                          >  6:315105/1‑106 (MQ=255)                                                                                                                      tGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGcttct                              <  7:148824/149‑1 (MQ=255)                                                                                                                        cGGTACGGCTCCTGAGTACCCAATGTAATCAGGTCGTCGACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTGACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCttcttc                            <  4:89140/149‑1 (MQ=255)                                                                                                                         ggTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCg                           >  2:80294/1‑149 (MQ=255)                                                                                                                                gaccTGATTACCTTGGGTACTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTgg                                          >  7:84559/3‑127 (MQ=255)                                                                                                                                gaccTGATTACCTTGGGTACTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTgg                                          >  7:89708/3‑127 (MQ=255)                                                                                                                                gaccTGATTACCTTGGGTACTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTgg                                          <  8:89708/125‑1 (MQ=255)                                                                                                                                gaccTGATTACCTTGGGTACTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTgg                                          <  8:84559/125‑1 (MQ=255)                                                                                                                                gaccTGATTACCTTGGGTACTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCgg                    >  1:221708/3‑149 (MQ=255)                                                                                                                                     gAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGgtggtg               >  3:221757/1‑149 (MQ=255)                                                                                                                                          cccAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGtt          >  3:371011/1‑149 (MQ=255)                                                                                                                                          cccAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGt           <  1:80294/148‑1 (MQ=255)                                                                                                                                           ccAAGGTAATCTGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGtt          >  1:286216/1‑148 (MQ=255)                                                                                                                                           ccAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGtt          >  4:254315/1‑148 (MQ=255)                                                                                                                                                gTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTg    >  6:351112/1‑149 (MQ=255)                                                                                                                                                 tAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGc   <  5:292801/149‑1 (MQ=255)                                                                                                                                                  aaTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGcc  >  5:206018/1‑149 (MQ=255)                                                                                                                                                  aaTCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGcc  >  2:107630/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACA‑TGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAGGCCCTGGGCGCCGGCTTCTTCGTACCCGGTGGTGCCGTTGATCTGCC  >  NC_002947/2981‑3273

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: ‑