Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	44,310	A→G	N44S (AAC→AGC)	fixC →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	44,310	0	A	G	88.6%	95.0 / 5.7	35	N44S (AAC→AGC)	fixC	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (11/20);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

GCTGACAGGGAGCTCTTATGTCCGAAGATATCTTTGACGCCATCATCGTCGGTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAACGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGAAAAGTCAGCGATGACTATGGACTACTGCAATGG  >  NC_000913/44163‑44457                                                                                                                                                    |                                                                                                                                                    gCTGACAGGGAGCTCTTATGTCCGAAGATATCTTTGACGCCATCATCGTCGGTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGt                                                                                                                                                  <  2:628003/151‑1 (MQ=255)           aGCTCTTATGTCCGAAGATATCTTTGACGCCATCATCGTCGGTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCggcggg                                                                                                                                        <  2:376522/151‑1 (MQ=255)              tctTATGTCCGAAGATATCTTTGACGCCATCATCGTCGGTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCggcggg                                                                                                                                        >  2:57521/1‑148 (MQ=255)                ttATGTCCGAAGATATCTTTGACGCCATCATCGTCGGTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCggcggg                                                                                                                                        >  2:399523/1‑146 (MQ=255)                 tATGTCCGAAGATATCTTTGACGCCATCATCGTCGGTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGtct                                                                                                                                   >  1:298860/1‑150 (MQ=255)                  aTGTCCGAAGATATCTTTGACGCCATCATCGTCGGTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGtct                                                                                                                                   >  2:319979/1‑149 (MQ=255)                   tGTCCGAAGATATCTTTGACGCCATCATCGTCGGTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGGGcctc                                                                                                                                  >  2:323649/1‑149 (MQ=255)                      ccGAAGATATCTTTGACGCCATCATCGTCGGTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAACGTCACCGGCGGGCGTCTCTAt                                                                                                                               >  1:276131/1‑149 (MQ=255)                                   tGACGCCATCATCGTCGGTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGa                                                                                                                 <  1:345501/150‑1 (MQ=255)                                         catcatCGTCGGTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACatt                                                                                                          <  2:106691/151‑1 (MQ=255)                                                    gTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGttt                                                                                               <  1:525276/151‑1 (MQ=255)                                                               gCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCg                                                                                    <  2:316806/151‑1 (MQ=255)                                                                       ggttGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGGGTCTCTATGGCCACAGTCTGGGACACATTATTCCTGGTTTTGCCGGCTcccccccc                                                                               >  1:47605/1‑148 (MQ=255)                                                                              gCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGcc                                                                     <  1:256260/151‑1 (MQ=255)                                                                                                  aGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTc                                                 <  1:594160/151‑1 (MQ=255)                                                                                                    gTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAACGTCACCGGCGGGGGTCTCTATGCCCACAGTCTGGAAAACATTATTCCTGGTTTCGCCGACTCCCCCCCCGTAGAACGCCTGATCACCCATGAAAAACTcgc                                               >  2:557899/1‑151 (MQ=255)                                                                                                           aGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATg                                        <  1:151397/151‑1 (MQ=255)                                                                                                            gtgtTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGa                                       <  2:184354/151‑1 (MQ=255)                                                                                                                tagttaTCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTAGGCCCACCATCGGGGACAAAATATTTCTCGGTTCGCCGGCTTCGCCCCCCTAGCAAGCCCCGACACCCAAGAAAAAAACGCGATTTTg                                        >  2:68884/1‑146 (MQ=255)                                                                                                                 agttaTCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGa                                   <  1:385334/150‑1 (MQ=255)                                                                                                                    taTCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTTGCCGACTCCGCCCCCGTAGAAAGCCTGATCACCCATGAAAAAATCGCGGTTTTGACGGgaaa                                >  1:131062/1‑150 (MQ=255)                                                                                                                     aTCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGAAAAg                               <  2:381915/150‑1 (MQ=255)                                                                                                                     aTCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGAAAAGt                              <  1:401370/151‑1 (MQ=255)                                                                                                                     aTCGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGGTTCGCCGACTCCCCCCCCCTAGAACGCCTGATCACCCATGAAAAAATCGCGGTTTTGACGGgaa                                 >  2:83769/1‑148 (MQ=255)                                                                                                                     aTCGAGCGTGGCAATTCCGCAGGTGCCAAGAACGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGAAAAg                               <  1:557910/150‑1 (MQ=255)                                                                                                                       cGAGCGTGGCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGAAAAGTCa                            <  1:245720/151‑1 (MQ=255)                                                                                                                              ggCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCGCCCATGAAAAACTCGCGTTTATGACGGAAAAGTCAGCGATg                      <  2:294667/150‑1 (MQ=255)                                                                                                                               gCAATTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGAAAAGTCAGCGATGAc                    <  2:170981/151‑1 (MQ=255)                                                                                                                                 aaTTCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCCCAGTCTGGGACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGGAGAAAGCCTTATCACCCCTGAAAAACTCGCGTTTATGACGGAAAAGCCAg                           >  2:91431/1‑142 (MQ=255)                                                                                                                                    tCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGAAAAGTCAGCGATGACTATg                <  2:464160/150‑1 (MQ=255)                                                                                                                                    tCCGCAGGTGCCAAGAGCGTCACCGGCGGGCGGCTCTATTCCCACAGTCTGGGAAACATTATTCCTGGGTTTGCCGAATCCGCCCCCGTAGAACGCCTGATCACCCCTGAAAAACTCGCGTTTATGACGGGAAAATCAGCGATGACt                   >  2:19498/1‑147 (MQ=255)                                                                                                                                      cGCAGGTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGAAAAGTCAGCGATGACTATGGac             <  1:7882/151‑1 (MQ=255)                                                                                                                                           gTGCCAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGAAAAGTCAGCGATGACTATGGACTACTg        <  2:368734/151‑1 (MQ=255)                                                                                                                                              ccAAGAGCGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGAAAAGTCAGCGATGACTATGGACTACTGCa      <  1:520081/150‑1 (MQ=255)                                                                                                                                                 agagCGTCACCGGCGGGGGTCTCTATGCCCCCCGTCTGGAACACATTATTTCTGGGTTTGCCGACTCCGCCCCCGTAGAAAGCCCGATCACCCATGAAAAACTCGCGGTTATGACGGAAAAGTCAGCGATGACTATGGACCACTTCCATgg  >  2:153334/1‑151 (MQ=255)                                                                                                                                                 agagCGTCACCGGCGGGCGGCTCTATGCCCACAGTCTGGAACACATTATTTCTGGGTTCGCCGACTCCCCCCCCGCAGAACGCCTGATCACCCCTGAAAAACTCGCGGTTTTGACGGAAAAGTCAGCGATGAATATGGACTACTg        >  2:426404/1‑145 (MQ=255)                                                                                                                                                 aGAACGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGAAAAGTCAGCGATGACTATGGACTACTGCAATgg  <  2:276124/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCTGACAGGGAGCTCTTATGTCCGAAGATATCTTTGACGCCATCATCGTCGGTGCAGGGCTTGCCGGTTCGGTTGCCGCACTGGTGCTCGCCCGCGAAGGTGCGCAAGTGTTAGTTATCGAGCGTGGCAATTCCGCAGGTGCCAAGAACGTCACCGGCGGGCGTCTCTATGCCCACAGTCTGGAACACATTATTCCTGGTTTCGCCGACTCCGCCCCCGTAGAACGCCTGATCACCCATGAAAAACTCGCGTTTATGACGGAAAAGTCAGCGATGACTATGGACTACTGCAATGG  >  NC_000913/44163‑44457

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