Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,157,282	A→C	D169A (GAT→GCT)	ycfH →	putative DNase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,157,282	0	A	C	86.4%	119.5 / 0.1	45	D169A (GAT→GCT)	ycfH	putative DNase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (11/27);  minor base G (6/0);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 1.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAA  >  NC_000913/1157144‑1157425                                                                                                                                           |                                                                                                                                                 tgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTAATGGCCCT‑Cggggttt                                                                                                                                      >  2:743738/1‑145 (MQ=255)  gaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTAc                                                                                                                                    <  2:840192/151‑1 (MQ=255)       aGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACAtctc                                                                                                                               <  2:646308/151‑1 (MQ=255)                 cGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCgg                                                                                                                      <  1:59937/150‑1 (MQ=255)                     cATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATtg                                                                                                                 <  2:875136/151‑1 (MQ=255)                     cATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATtg                                                                                                                 <  1:819798/151‑1 (MQ=255)                      aTACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATtgt                                                                                                                <  2:137519/151‑1 (MQ=255)                             tGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAAATACTGGCTCT‑CGGGTTTTATATCTCCTTTTTCGGCCTTGTGACCtt                                                                                                          >  2:815240/1‑150 (MQ=255)                             tGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACAGGTCCT‑CGGTTTTTAAATCTTCTTTTTCGCCCTTTTGAccct                                                                                                          >  1:869861/1‑148 (MQ=255)                                  cccGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTaa                                                                                                    <  2:926586/151‑1 (MQ=255)                                   ccGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAAt                                                                                                   <  2:305705/151‑1 (MQ=255)                                   ccGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAAt                                                                                                   <  2:602092/151‑1 (MQ=255)                                       gCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCg                                                                                                <  1:107123/150‑1 (MQ=255)                                             aCGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTATTGGCTCT‑CGGATTTTAAATCTCCTTTTTCGGCACTGTGACCTTTCGTAAAGCGGGGCaa                                                                                          >  1:18027/1‑150 (MQ=255)                                             aCGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTTGCTCT‑CCGATTTTAAATCTCCTTTTTCCGCCTTGTGACCTTCCCTAATTCgggggaa                                                                                          >  2:580125/1‑145 (MQ=255)                                              cGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACt                                                                                        <  1:740550/151‑1 (MQ=255)                                              cGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTTGCCCT‑CGGGTTTTTTATCTCCTTTTTCCGCCCTGTGGCCCTCCCCCAAGCGGAGc                                                                                            >  2:881007/1‑147 (MQ=255)                                              cGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATGGGTCT‑CGGGTTTTAAAATTCCTTTTTCCGCCTTGTGGACTTCCGTAAAGGGGGGCaa                                                                                          >  1:161430/1‑149 (MQ=255)                                               gCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATAATTGCTCT‑CCGGTTTTAAATCTCCCTTTTCGGCATTTTGACCTTCCGTAATGCGGAGCAACTg                                                                                       >  2:109484/1‑151 (MQ=255)                                                 tGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAAATTCTGGGGCT‑CGGGTTTTACAACTCCTTTTTCGGCATTGTGACCTTTCGTAATGCGGAGCaaatccg                                                                                     >  2:807750/1‑146 (MQ=255)                                                          tGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGACTCCGGCTTTTAAAATCCCTTTTCCGCTCTTTTTA‑CTTTCGTGATGGGGAGAGACTGcgcggatct                                                                              >  1:244773/1‑144 (MQ=255)                                                           gcgcGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGAtgctgc                                                                            <  2:640733/150‑1 (MQ=255)                                                             gcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGcgcg                                                                         <  2:352593/151‑1 (MQ=255)                                                               gaagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGtt                                                                       <  1:658735/151‑1 (MQ=255)                                                                aagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGGCCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTa                                                                      <  2:167607/151‑1 (MQ=255)                                                                aagaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTa                                                                      <  1:226398/151‑1 (MQ=255)                                                                  gaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATg                                                                    <  2:609811/151‑1 (MQ=255)                                                                         tGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTCT‑CGGGTTTTAAATCTCCTTTTTCGGCATTGTGACCTTCCGTAATTCGGGGCAACTGCGCGGTGCTGCGCGTTATGTccccc                                                              >  2:89604/1‑150 (MQ=255)                                                                          gACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCt                                                             <  1:208121/150‑1 (MQ=255)                                                                             ggATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTAATGGGTTT‑CGGGTTTTAAATCTCCCTTTTTGGCATTTTGTACTTCCCTCATGCGGGGCAACTGCGCGATGCTGCGCGTTATGTgcccc                                                              >  2:521756/1‑146 (MQ=255)                                                                                 tGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCgg                                                      <  2:599786/150‑1 (MQ=255)                                                                                   cGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTTAAATTCTTGCTTT‑CGGATTTTTAATCTCCTTTTTCGGCAATGTGGCCTTTCGCAATACGGGGCAAACGCGCGCTGGTGCGCGGTATGTTCCCCtcggcccggtt                                                   >  2:359553/1‑141 (MQ=255)                                                                                      tGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAAATACTTGCTCT‑CCGGGTTTAAATCTCCTTTTTCCGCATTTTGACCTTTCGTAATGCGGGGCAACCGCGCGGTGCTGCGCGGTATGTTCCCCCGGGTCgggttct                                                 >  2:85913/1‑145 (MQ=255)                                                                                        gCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACtggt                                              <  1:483680/151‑1 (MQ=255)                                                                                            actacACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTTGCTCT‑CGGGTTTTAAATCTCCTTTTTCCGCATTGTGACCCTCCGTAATTCGGAGCAACTGCGCGATGCTGCGCGGTTTGTTCCCCCGGGTCGGTTACtggtgg                                            >  2:685979/1‑149 (MQ=255)                                                                                            actacACTGTTTTACAGAGGACAGAGAAACGGCGGGGAAATTACTGGGTCT‑CGGGTTTTACATCTCCCTTTTCTGCATTGTGAACTTCCGTAATGGGGAGGAACTGCGCGATGCTGCGCGTTATGTTCCCCTCGGACGGTTACtgggg                                             >  2:368195/1‑146 (MQ=255)                                                                                                           agagGACAGAGAAACGGCGGGGAAAATACTTGCTCT‑CGGGTTTTAAATCTCCTTTTTCGGCATTGTGACCTTTCCTAATTCGGGGCAAATGCGCGGTGCTGCGCGTTATGTCCCCCCGGGTCGGTTACTGGTGGAAAATGACTCacctt                             >  2:340017/1‑149 (MQ=255)                                                                                                              ggACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCaccttacctt                        <  1:333998/151‑1 (MQ=255)                                                                                                               gACAGAGAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTg                       <  1:932841/151‑1 (MQ=255)                                                                                                                  agagAAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTgcgc                    <  2:778205/151‑1 (MQ=255)                                                                                                                    agaAACGGCGGGTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCg                  <  2:157146/151‑1 (MQ=255)                                                                                                                              ggTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATc        <  1:348472/151‑1 (MQ=255)                                                                                                                               gTAAATTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCg       >  1:531670/1‑151 (MQ=255)                                                                                                                                   aTTACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGa   <  2:18028/151‑1 (MQ=255)                                                                                                                                     tACTGGCTCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGaa  <  1:822412/150‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                 TGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGTAAATTACTGGATCT‑CGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAA  >  NC_000913/1157144‑1157425

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