Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	698,669	C→T	8.3%	R170H (CGC→CAC)	asnB ←	asparagine synthetase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	698,669	0	C	T	8.3%	143.0 / 3.5	48	R170H (CGC→CAC)	asnB	asparagine synthetase B
Reads supporting (aligned to +/- strand):  ref base C (21/23);  new base T (2/2);  total (23/25)
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.41e-01

TGAATCTTCCAGTGCCTGACGCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGC  >  NC_000913/698521‑698811                                                                                                                                                     |                                                                                                                                               tGAATCTTCCAGTGCCTGACGCAGCTCGTTTTTGTCGGTCACGTTATCTTTAACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCACCGCCGGGAGCTCTTTAATCGTGc                                                                                                                                                >  6:84905/1‑149 (MQ=255)   aaTCTTCCAGTGCCTGACGCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCgg                                                                                                                                              >  6:12440/1‑149 (MQ=255)   aaTCTTCCAGTGCCTGACGCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCgg                                                                                                                                              >  6:299736/1‑149 (MQ=255)   aaTCTTCCAGTGCCTGACGCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCg                                                                                                                                               <  3:83728/148‑1 (MQ=255)             tGCCTGACGCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCa                                                                                                                                    <  1:154583/149‑1 (MQ=255)                 tGACGCAGCTCGTTTTTGTCGGTCACGTTATATTTCACCGAATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCATTGCTCCACAAATAGCTCCCCGCCGGGAAATCTTTAATCGTCCGGCAAAATGGCAc                                                                                                                                   >  5:134301/1‑146 (MQ=255)                       aGCTCGGTTTTGTCGGTCCCGGTATCTTTCGCCGCATCGGAGTCGAACCAGTCGCGATGCTAGTAAGAACGGATTTCGCCGTCCTGGCTCGCCAAATTGCTCCCCGCCGGGCACTCTTGGATCGTGCGGCTACCTGGCACCAGCGCgtt                                                                                                                          >  2:237293/1‑149 (MQ=255)                        gCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGAGTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCTGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTc                                                                                                                         <  1:83387/149‑1 (MQ=255)                         cTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTAGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGAGCGGCAAACTGGCACCAGCGCTTTCa                                                                                                                        >  5:134210/1‑149 (MQ=255)                          acGTTTTTGTCGGTCACGTTATCTTTCACAGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACCGATTTCGCCGTCCTGGATCCACAAATAGCTCCCCGCCTGGAACTCTTTAATCGTCCGGCAAACTGGCACCAGCGCTTTCAt                                                                                                                       >  2:1269/2‑149 (MQ=255)                           cGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCAtt                                                                                                                      >  7:239745/1‑149 (MQ=255)                            gTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCAttt                                                                                                                     <  8:270701/149‑1 (MQ=255)                            gTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCAtt                                                                                                                      <  7:14978/148‑1 (MQ=255)                            gTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCAt                                                                                                                       <  4:292744/147‑1 (MQ=255)                            gTTTTTGTCGGCCACGTTATCTTTCACCGCATCGTAGGCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCAttt                                                                                                                     <  7:309634/149‑1 (MQ=255)                               tttGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCt                                                                                                                   <  7:142995/148‑1 (MQ=255)                                  gTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAg                                                                                                                <  5:12440/148‑1 (MQ=255)                                         aCGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACAt                                                                                                         <  3:176019/148‑1 (MQ=255)                                         aCGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATa                                                                                                        <  5:299736/149‑1 (MQ=255)                                                             gTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGTGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCa                                                                                     >  5:263917/1‑148 (MQ=255)                                                             gTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCATTAATCGTGTGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCa                                                                                     >  3:204119/1‑148 (MQ=255)                                                                 tCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATAccc                                                                                >  5:183418/1‑149 (MQ=255)                                                                   gatccAGTCGCGATGATAGTAAGAAAGGATTTCGCCGTCCTGGCTCCACAAATAGCT‑CCAGCCGGTATCTCTTTAATCGTGCGGCAAACTGGCACCATCGCTTTCATTTCTGAGGCCACATACAGCTGACGGTGTTCGTCATACCCCat                                                                             <  4:65675/146‑1 (MQ=255)                                                                     aCCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCata                                                                            <  7:96129/149‑1 (MQ=255)                                                                        aGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACAAACAGCTCACCGTGTTCGTCATACACCAtata                                                                          >  4:98039/1‑148 (MQ=255)                                                                         gTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACa                                                                        <  6:134210/149‑1 (MQ=255)                                                                                tgatAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGgat                                                                 <  1:128114/149‑1 (MQ=255)                                                                                 gatAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGgat                                                                 >  3:125204/1‑148 (MQ=255)                                                                                                              cTCCACAAATCGCTCCCCGCCGGGAACTCTTTAATCGGGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTa                                   <  5:84905/149‑1 (MQ=255)                                                                                                                 cacaAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAgg                                 >  5:85428/1‑148 (MQ=255)                                                                                                                 cacaAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAgg                                 >  2:210939/1‑148 (MQ=255)                                                                                                                 cacaAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGc                                <  5:139363/149‑1 (MQ=255)                                                                                                                    aaaTAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATc                             >  1:78245/1‑149 (MQ=255)                                                                                                                    aaaTAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATc                             >  1:340541/1‑149 (MQ=255)                                                                                                                       tAGCTCCCCGCCGGGAACTCTTTAATTGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCCCATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCTCGACCAATCAGGTAGGCATCttt                          <  5:118430/149‑1 (MQ=255)                                                                                                                         gCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCttttt                        <  6:183418/149‑1 (MQ=255)                                                                                                                          cTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTc                       >  7:57592/1‑149 (MQ=255)                                                                                                                             cccGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGc                     <  8:239745/148‑1 (MQ=255)                                                                                                                             cccGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCGTTTTCGCt                    >  4:200308/1‑149 (MQ=255)                                                                                                                              ccgccgGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGCTCCCCAGATGGACGCGACCAATCAGGTAGGCATCTTTTTCGCTg                   >  6:298278/1‑149 (MQ=255)                                                                                                                                   gggAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTa              <  8:111176/149‑1 (MQ=255)                                                                                                                                    ggAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTgg                                                                    >  7:312357/1‑94 (MQ=255)                                                                                                                                    ggAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTgg                                                                    <  8:312357/94‑1 (MQ=255)                                                                                                                                     gAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTa              >  4:271644/1‑147 (MQ=255)                                                                                                                                              aaTCGTGTGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAg    <  6:263917/148‑1 (MQ=255)                                                                                                                                              aaTCGTGTGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACACCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAg    <  4:204119/148‑1 (MQ=255)                                                                                                                                               aTCGTGCGGCAATCTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGc  <  3:98039/149‑1 (MQ=255)                                                                                                                                               aTCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATCCCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAgg   >  5:188917/1‑148 (MQ=255)                                                                                                                                                tCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGCTGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGc  >  8:33859/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                               TGAATCTTCCAGTGCCTGACGCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGATTTCGCCGTCCTGGCTCCACAAATAGCTCCCCGCCGGGAACTCTTTAATCGTGCGGCAAACTGGCACCAGCGCTTTCATTTCTGAGGCCACATACAGCTGACCGTGTTCGTCATACCCCATATACAGTGGGATGATCCCCAGATGGTCGCGACCAATCAGGTAGGCATCTTTTTCGCTGTCGTACAGTGCAAAGGC  >  NC_000913/698521‑698811

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