Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,326,925	(C)6→7	coding (817/1185 nt)	atoB →	acetyl‑CoA acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,326,919	1	.	C	97.1%	119.2 / ‑3.2	35	coding (811/1185 nt)	atoB	acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (16/18);  total (16/19)
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.93e-01

AAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑A‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGC  >  NC_000913/2326769‑2327064                                                                                                                                                         |                                                                                                                                                  aaaGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑a                                                                                                                                                    <  1:305284/151‑1 (MQ=255)      gAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑A‑ccccc                                                                                                                                              >  2:169432/1‑151 (MQ=255)      gAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑A‑ccccc                                                                                                                                              <  1:578680/151‑1 (MQ=255)       aaTTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑A‑cccccc                                                                                                                                             <  2:416949/151‑1 (MQ=255)       aaTTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑A‑cccccc                                                                                                                                             >  2:308918/1‑151 (MQ=255)          tCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑A‑ccccccctg                                                                                                                                          >  2:117298/1‑148 (MQ=255)          tCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑A‑ccccccct                                                                                                                                           <  2:407416/150‑3 (MQ=255)                           aGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTa                                                                                                                          <  1:337421/151‑1 (MQ=255)                             gTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATg                                                                                                                        <  1:800796/151‑1 (MQ=255)                                     cgcCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGt                                                                                                                >  2:236452/1‑151 (MQ=255)                                      gcCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTg                                                                                                               <  2:158591/151‑1 (MQ=255)                                        ccGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGc                                                                                                             >  2:752415/1‑151 (MQ=255)                                            ccTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGc                                                                                                         <  1:156925/151‑1 (MQ=255)                                              ttCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGccc                                                                                                       >  2:696514/1‑151 (MQ=255)                                                         aGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGAt                                                                                            >  1:91091/1‑151 (MQ=255)                                                                 tCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATgg                                                                                    >  2:338911/1‑151 (MQ=255)                                                                   aCCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATgggg                                                                                  >  1:95643/1‑151 (MQ=255)                                                                            aaCGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTAcc                                                                          >  2:702084/1‑150 (MQ=255)                                                                             aCGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTg                                                                        <  1:191242/151‑1 (MQ=255)                                                                                 gTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCAc                                                                    >  2:587040/1‑151 (MQ=255)                                                                                        aTTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCaaaa                                                              >  1:128897/1‑150 (MQ=255)                                                                                          tAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAg                                                            <  2:453611/150‑1 (MQ=255)                                                                                          tAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGc                                                           >  2:510507/1‑151 (MQ=255)                                                                                            acgacgGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGt                                                         >  2:820718/1‑151 (MQ=255)                                                                                                    gctgcCGCTCTGGTTATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACt                                                  <  1:644209/150‑1 (MQ=255)                                                                                                    gctgcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGGTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACt                                                  <  2:421616/150‑1 (MQ=255)                                                                                                       gcCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCg                                              <  2:666215/151‑1 (MQ=255)                                                                                                         cGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCggg                                            <  1:721440/151‑1 (MQ=255)                                                                                                         cGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCggg                                            <  1:595776/151‑1 (MQ=255)                                                                                                            tctGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACg                                                                   >  2:157093/1‑125 (MQ=255)                                                                                                            tctGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACg                                                                   <  1:157070/125‑1 (MQ=255)                                                                                                                  gATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACt                                   <  1:206935/151‑1 (MQ=255)                                                                                                                         gaagaaTCTGCGGCGCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGat                            <  1:251118/151‑1 (MQ=255)                                                                                                                                  gcggcgCTGGCAGCAGGCCTTGC‑CCCCCCCGGGTCGCATTAAAAGTTATGTCAGCGGTGTCGTGCCCACCGCATTGATGGGTATGGGGCGAGTTACGTCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGTCGGATATTGAt                      <  1:536056/148‑1 (MQ=255)                                                                                                                                  gcggcgCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGAtc                     <  1:637954/149‑1 (MQ=255)                                                                                                                                    ggcgCTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATCCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCCAAAAGCGTTACAACTGGCGGGGCTTCAACAGGCGGATATTGctct                    >  1:24828/1‑148 (MQ=255)                                                                                                                                       gcTGGCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTg               >  1:287756/1‑150 (MQ=255)                                                                                                                                          ggCAGCAGGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGGTATTGATCTCATTGAGGc           >  1:715402/1‑151 (MQ=255)                                                                                                                                              gcagGCCTT‑ACCCCCCCTGGCGCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTCCCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAg             <  1:378860/145‑1 (MQ=255)                                                                                                                                               cagGCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATg      <  1:476565/151‑1 (MQ=255)                                                                                                                                                  gCCTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAg   <  1:769838/151‑1 (MQ=255)                                                                                                                                                   ccTT‑ACCCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGc  >  1:237546/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                  AAAGCGAATTCAACGGCTGAAGCGTTAGGTGCATTGCGCCCGGCCTTCGATAAAGCAGGAACAGTCACCGCTGGGAACGCGTCTGGTATTAACGACGGTGCTGCCGCTCTGGTGATTATGGAAGAATCTGCGGCGCTGGCAGCAGGCCTT‑A‑CCCCCCTGGCTCGCATTAAAAGTTATGCCAGCGGTGGCGTGCCCCCCGCATTGATGGGTATGGGGCCAGTACCTGCCACGCAAAAAGCGTTACAACTGGCGGGGCTGCAACTGGCGGATATTGATCTCATTGAGGCTAATGAAGC  >  NC_000913/2326769‑2327064

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