Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,564,565	T→C	T673A (ACG→GCG)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,564,565	0	T	C	97.6%	143.1 / NA	42	T673A (ACG→GCG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (24/17);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACTGCGCCAGGCTAACGCTGTAGTT  >  NC_000913/3564438‑3564697                                                                                                                                |                                                                                                                                     tcGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAg                                                                                                                          >  2:591091/1‑140 (MQ=255)    tagtaTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAgcg                                                                                                                        <  1:688118/139‑1 (MQ=255)    tagtaTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAgcg                                                                                                                        >  2:688118/1‑139 (MQ=255)     agtaTTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCa                                                                                                                      >  1:138089/1‑140 (MQ=255)         ttCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACTTTCGCACCGTCAAACGCACCGATAGTCAGCGCACCGt                                                                                                                  >  1:288148/1‑140 (MQ=255)         ttCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGt                                                                                                                  >  2:533751/1‑140 (MQ=255)            aCGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTGGCACCGTTCAAAGCGCCGATAGTCAGCGCACCGTTaa                                                                                                               >  1:215137/1‑140 (MQ=255)             cgcgGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAg                                                                                                              >  1:576849/1‑140 (MQ=255)                    tGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGACAGCACACCGTTAAGCGCAAAc                                                                                                       >  2:274680/1‑140 (MQ=255)                        gCCCTGACGACGCAGTTCTTCCTCTTCTTCCCCTGAGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCa                                                                                                   <  2:288148/140‑1 (MQ=255)                             gacgacGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTa                                                                                              <  2:372595/140‑1 (MQ=255)                             gacgacGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTa                                                                                              <  1:156110/140‑1 (MQ=255)                                        tcttcCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGa                                                                                   <  1:35196/140‑1 (MQ=255)                                          ttcCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAg                                                                                 >  1:30763/1‑140 (MQ=255)                                             cACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCtt                                                                              <  2:152995/140‑1 (MQ=255)                                                 tcttcCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGt                                                                          >  1:618947/1‑140 (MQ=255)                                                  cttcCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTCAGCGCAACCTTCATGTAACTGGTGCCGGAAGCTTCCGTc                                                                         >  2:308694/1‑140 (MQ=255)                                                     ccGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCACAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCt                                                                      >  2:730409/1‑140 (MQ=255)                                                       gCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGc                                                                    <  2:465948/140‑1 (MQ=255)                                                           tgtTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGc                                                                >  2:709730/1‑140 (MQ=255)                                                             ttACCAAAAATAAAGATATTGTCAGAACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGa                                                              <  1:253018/140‑1 (MQ=255)                                                               acagaaaaTAAAGATATCGTCAACACCGACATGATCCAGAATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGaaa                                                            <  2:193056/136‑1 (MQ=255)                                                               aCCAAAAATAAAGATGTTGTCACCATCGACATGATCCCGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGaaa                                                            <  1:274680/140‑1 (MQ=255)                                                                ccAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAAt                                                           <  1:623765/140‑1 (MQ=255)                                                                     aaTAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGtt                                                      >  2:608812/1‑140 (MQ=255)                                                                       tAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTca                                                    <  2:138089/140‑1 (MQ=255)                                                                              atTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAGGCTTCCGGCCCAGCGAGCGAATTCTGTTcagacctgt                                             >  1:193056/1‑136 (MQ=255)                                                                                       cccGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGccgccg                                    <  2:487208/139‑1 (MQ=255)                                                                                       aCCGACATGCTCCAGCATCTCGACATTCGCACCGTCAAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGGCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGccgccg                                    <  2:734015/140‑1 (MQ=255)                                                                                                   cAGCATCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACt                        >  2:138307/1‑140 (MQ=255)                                                                                                       aTCTCGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACTgcgc                    >  1:97683/1‑140 (MQ=255)                                                                                                         ctcGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACTGCGCCa                  >  2:611043/1‑140 (MQ=255)                                                                                                          tcGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACTGCGCCAg                 >  1:601942/1‑140 (MQ=255)                                                                                                           cGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTTAGACAGACCTGGCGCCGGAATCATCAACTGCGCCAgg                >  1:275396/1‑140 (MQ=255)                                                                                                           cGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACTGCGCCAgg                >  2:44444/1‑140 (MQ=255)                                                                                                           cGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACTGCGCCAgg                >  2:582148/1‑140 (MQ=255)                                                                                                           cGACATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACTGCGCCAgg                >  2:60330/1‑140 (MQ=255)                                                                                                            gACATTCGCACCGTCAAACGCAACGCTAGCCAGCGCAGCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACTGCGCCAGGc               <  1:273443/140‑1 (MQ=255)                                                                                                              cATTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACTGCGCCAGGCTa             <  2:30763/140‑1 (MQ=255)                                                                                                               aTTCGCACCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACTGCGCCAGGCTaa            <  2:601942/140‑1 (MQ=255)                                                                                                                     aCCGTCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACTGCGCCAGGCTAACGCTGt      <  1:533751/140‑1 (MQ=255)                                                                                                                        gcccAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTCACTGGTGCCGGAGGCTTCCGCCCCTGCCAGCGAAATCTGTTCATACAGATCTCCCGGCGGCATGATCAAGTGCGCCAGGCTAACGCTGTAGt   >  2:293221/3‑140 (MQ=255)                                                                                                                         tCCAACGCACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGGCGCCGGAGTGATCAACGGCGCCAGGCTAACGCTGTAGtt  >  2:277291/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                     TCGTAGTATTCACGCGGTTTGTAGCCCTGACGACGCAGTTCTTCCACTTCTTCCGCTGTGTTACCAAAAATAAAGATATTGTCAGCACCGACATGATCCAGCATCTCGACATTCGCACCGTCCAACGTACCGATAGTCAGCGCACCGTTAAGCGCAAACTTCATGTTACTGGTGCCGGAAGCTTCCGTCCCTGCCAGCGAAATCTGTTCAGACAGATCTGCCGCCGGAATGATCAACTGCGCCAGGCTAACGCTGTAGTT  >  NC_000913/3564438‑3564697

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