Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,792,767	A→T	18.0%	intergenic (‑216/+59)	kbl ← / ← yibB	glycine C‑acetyltransferase/YibB family protein, function unknown

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,792,767	0	A	T	18.0%	125.8 / 22.3	50	intergenic (‑216/+59)	kbl/yibB	glycine C‑acetyltransferase/YibB family protein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base T (4/5);  total (25/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 = 8.69e-01

CGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATC  >  NC_000913/3792620‑3792903                                                                                                                                                    |                                                                                                                                         cgTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGGTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaa                                                                                                                                          <  3:1093268/148‑1 (MQ=255)          ttCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGGAAAAATGATGCTTTGGCATCAtt                                                                                                                                 <  3:546381/148‑1 (MQ=255)                                aaCGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAg                                                                                                                    >  4:1071683/1‑139 (MQ=255)                                aaCGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAg                                                                                                                    <  3:1071683/139‑1 (MQ=255)                                           tGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                 >  1:125992/1‑115 (MQ=255)                                           tGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                 <  2:125992/115‑1 (MQ=255)                                                tctctcAACGCCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAGTGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca                                                                                                                                      >  4:739232/1‑105 (MQ=255)                                                tctctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                 <  3:924451/110‑1 (MQ=255)                                                tctctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                   <  1:1067939/124‑1 (MQ=255)                                                tctctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                   >  2:1067939/1‑124 (MQ=255)                                                tctctcAACACCAACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca                                                                                                                                      <  3:739232/105‑1 (MQ=255)                                                 ctctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaa                                                                                                                                          <  1:291747/100‑1 (MQ=255)                                                 ctctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaa                                                                                                                                          >  2:291747/1‑100 (MQ=255)                                                   ctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTAt                                                                                                             >  1:93771/1‑127 (MQ=255)                                                   ctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTAt                                                                                                             <  2:93771/127‑1 (MQ=255)                                                       acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                <  1:5389/104‑1 (MQ=255)                                                       acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                >  2:5389/1‑104 (MQ=255)                                                       acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatc                                                                                                                                    <  1:211799/100‑1 (MQ=255)                                                       acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatc                                                                                                                                    >  2:211799/1‑100 (MQ=255)                                                               aGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaa                                                                                                                                          <  3:797379/86‑1 (MQ=255)                                                               aGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaa                                                                                                                                          >  4:797379/1‑86 (MQ=255)                                                                 aTAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                   >  4:982843/1‑107 (MQ=255)                                                                 aTAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                   <  3:982843/107‑1 (MQ=255)                                                                      agAAAAACTGTACAGTTTTGGCTGTGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaa                                                                                                                                          <  3:587156/79‑1 (MQ=255)                                                                      agAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGCTTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaa                                                                                                                                          >  4:587156/1‑79 (MQ=255)                                                                       gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaa                                                                                                                                          >  2:27178/1‑78 (MQ=255)                                                                       gAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCaa                                                                                                                                          <  1:27178/78‑1 (MQ=255)                                                                        aaaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTcc                                                                                                                            >  3:644204/1‑91 (MQ=255)                                                                        aaaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTcc                                                                                                                            <  4:644204/91‑1 (MQ=255)                                                                         aaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                 >  3:424388/1‑85 (MQ=39)                                                                               gTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAg                                                                                                                                        <  4:25986/72‑1 (MQ=255)                                                                               gTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAg                                                                                                                                        >  3:25986/1‑72 (MQ=255)                                                                                                      gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAg                                                                                      <  2:475409/99‑1 (MQ=255)                                                                                                      gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAg                                                                                      >  1:475409/1‑99 (MQ=255)                                                                                                       atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaa                                                                                  >  1:946134/1‑102 (MQ=255)                                                                                                       atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaa                                                                                  <  2:946134/102‑1 (MQ=255)                                                                                                       atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                       >  4:837206/1‑145 (MQ=255)                                                                                                       atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                       <  3:837206/145‑1 (MQ=255)                                                                                                                               taAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       <  2:914553/89‑1 (MQ=255)                                                                                                                               taAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       >  1:914553/1‑89 (MQ=255)                                                                                                                                 aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtg                            <  2:304532/130‑1 (MQ=255)                                                                                                                                 aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtg                            >  1:304532/1‑130 (MQ=255)                                                                                                                                 aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTCGaa                                                                   <  2:1013613/91‑1 (MQ=255)                                                                                                                                 aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTCGaa                                                                   >  1:1013613/1‑91 (MQ=255)                                                                                                                                  aaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaa                                                                   <  2:940836/90‑1 (MQ=255)                                                                                                                                  aaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaa                                                                   >  1:940836/1‑90 (MQ=255)                                                                                                                                   aaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       <  1:165574/85‑1 (MQ=255)                                                                                                                                   aaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       >  2:165574/1‑85 (MQ=255)                                                                                                                                         aaaTGATGCTTTGGCATCATTTTTCCTATCTCAGTTGGTATCATAACGCTACTCACTAATGGAGGGAAACTTATTTTCTTGaa                                                                   >  2:229890/1‑83 (MQ=37)                                                                                                                                         aaaTGATGCTTTGGCATCATTTTTCCTATCTCAGTTGGTATCATAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaa                                                                   <  1:229890/83‑1 (MQ=38)                                                                                                                                         aaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                       <  4:802002/111‑1 (MQ=255)                                                                                                                                         aaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                       >  3:802002/1‑111 (MQ=255)                                                                                                                                         aaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAGATc  >  4:219191/1‑148 (MQ=255)                                                                                                                                          aaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       >  2:431772/1‑78 (MQ=255)                                                                                                                                          aaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       <  1:431772/78‑1 (MQ=255)                                                                                                                                              atgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCt                                                                                                    >  2:395811/1‑45 (MQ=255)                                                                                                                                              atgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCt                                                                                                    <  1:395811/45‑1 (MQ=255)                                                                                                                                                gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                 >  4:192038/1‑110 (MQ=255)                                                                                                                                                gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                 <  3:192038/110‑1 (MQ=255)                                                                                                                                                   aaaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       >  1:783164/1‑69 (MQ=255)                                                                                                                                                   aaaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       <  2:783164/69‑1 (MQ=255)                                                                                                                                                    aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATtc                                                         <  4:754307/82‑1 (MQ=255)                                                                                                                                                    aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATtc                                                         >  3:754307/1‑82 (MQ=255)                                                                                                                                                    aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAAt   >  3:324479/1‑136 (MQ=255)                                                                                                                                                    aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTCTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCACTAAAGTGTTACTCCTGAAGCAGCGAAACaatt   <  4:324479/136‑3 (MQ=255)                                                                                                                                                    |                                                                                                                                         CGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATC  >  NC_000913/3792620‑3792903

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