Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	512,212	T→G	8.0%	V191G (GTG→GGG)	glsA →	glutaminase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	512,212	0	T	G	8.0%	167.3 / 7.4	50	V191G (GTG→GGG)	glsA	glutaminase 1
Reads supporting (aligned to +/- strand):  ref base T (26/20);  new base G (2/2);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

TCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATG  >  NC_000913/512065‑512360                                                                                                                                                    |                                                                                                                                                     tCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGTAATGGAAGCCTGTGACgtg                                                                                                                                                     <  3:269352/149‑1 (MQ=255)    aaCAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGtat                                                                                                                                                  >  2:72235/1‑149 (MQ=255)    aaCAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGtat                                                                                                                                                  >  1:268964/1‑149 (MQ=255)             tggctggcGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAg                                                                                                                                         >  3:27657/1‑149 (MQ=255)               gctggcGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTg                                                                                                                                       <  3:262596/149‑1 (MQ=255)                         aGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGctc                                                                                                                             >  1:371998/1‑149 (MQ=255)                               cgcTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAAt                                                                                                                       <  4:393809/149‑1 (MQ=255)                                gcTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATa                                                                                                                      <  4:427624/149‑1 (MQ=255)                                gcTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATa                                                                                                                      >  3:198380/1‑149 (MQ=255)                                gcTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATa                                                                                                                      <  3:418267/149‑1 (MQ=255)                                  tctctGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACt                                                                                                                    >  4:231424/1‑149 (MQ=255)                                  tctctGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACt                                                                                                                    >  3:49991/1‑149 (MQ=255)                                    tctGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                  <  2:268964/149‑1 (MQ=255)                                        aCGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGGGTATAc                                                                                                                                                >  4:158335/1‑115 (MQ=255)                                        aCGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGGGTATAc                                                                                                                                                <  3:158335/115‑1 (MQ=255)                                          gAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGagga                                                                                                            <  2:151653/149‑4 (MQ=255)                                            aGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTgg                                                                                                          <  3:231424/149‑1 (MQ=255)                                              tCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAtt                                                                                                                 >  1:462314/1‑140 (MQ=255)                                              tCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAtt                                                                                                                 <  2:462314/140‑1 (MQ=255)                                              tCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCa                                                                                                        >  3:142176/1‑149 (MQ=255)                                                       cGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGgc                                                                                               >  2:90537/1‑149 (MQ=255)                                                         gAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGgcgc                                                                                             >  4:334612/1‑149 (MQ=255)                                                          aacaaacaACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGcgcg                                                                                            >  2:354310/1‑149 (MQ=255)                                                          aacaaacaACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGcgcg                                                                                            >  2:238977/1‑149 (MQ=255)                                                            caaacaaCCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGAc                                                                                          >  1:152351/1‑149 (MQ=255)                                                             aaacaaCCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACg                                                                                         >  2:26642/1‑149 (MQ=255)                                                             aaacaaCCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACg                                                                                         >  2:107919/1‑149 (MQ=255)                                                              aacaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGc                                                                                        <  3:237716/149‑1 (MQ=255)                                                              aacaacCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGc                                                                                        >  4:237716/1‑149 (MQ=255)                                                                 aaccaacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTgg                                                                                     >  2:243590/1‑149 (MQ=255)                                                                           cATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtgg                                                                           <  1:243590/149‑1 (MQ=255)                                                                            aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                          <  4:117629/149‑1 (MQ=255)                                                                              aaCCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGgtgt                                                                        <  1:55430/149‑1 (MQ=255)                                                                                 cGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGaa                                                                     <  4:27657/149‑1 (MQ=255)                                                                                      cATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGt                                                                >  4:291659/1‑149 (MQ=255)                                                                                      cATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGt                                                                >  3:372070/1‑149 (MQ=255)                                                                                           ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACg                                                           >  4:181280/1‑149 (MQ=255)                                                                                                     gTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcgcg                                                 <  4:198380/149‑1 (MQ=255)                                                                                                       aCTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAAGCGtt                                               <  1:90537/149‑1 (MQ=255)                                                                                                         tCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGttct                                             >  4:232753/1‑149 (MQ=255)                                                                                                            gccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAACCGCGTTCTTCa                                          <  1:354310/149‑1 (MQ=255)                                                                                                                    tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACa                                  >  1:218589/1‑149 (MQ=255)                                                                                                                             tgATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGt                         <  3:334612/149‑1 (MQ=255)                                                                                                                              gATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTa                        <  4:372070/149‑1 (MQ=255)                                                                                                                                     tGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACa                      <  3:460995/144‑1 (MQ=255)                                                                                                                                     tGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACa                      >  4:460995/1‑144 (MQ=255)                                                                                                                                          gCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACTCGTTCTTCAGGGCGACAACGTGCCGTACATTCTGGCCGa            >  4:213148/1‑149 (MQ=255)                                                                                                                                            cactgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTg                 <  2:270417/139‑1 (MQ=255)                                                                                                                                            cactgACGGGTATACACGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTg                 >  1:270417/4‑142 (MQ=255)                                                                                                                                             tgtgACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAat         <  4:49991/149‑1 (MQ=255)                                                                                                                                               tgACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAatga       >  2:2515/1‑149 (MQ=255)                                                                                                                                                   gtgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGt                              >  4:236358/1‑122 (MQ=255)                                                                                                                                                   gtgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGt                              <  3:236358/122‑1 (MQ=255)                                                                                                                                                   gtgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAatgatgat   >  4:382142/1‑149 (MQ=255)                                                                                                                                                    tgtATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAatgatgatg  <  4:47462/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATG  >  NC_000913/512065‑512360

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