Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	406,961	C→G	A353G (GCA→GGA)	atoB →	acetyl‑CoA C‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	406,961	0	C	G	98.2%	207.6 / NA	55	A353G (GCA→GGA)	atoB	acetyl‑CoA C‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (28/26);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAATTAAGCGATATTGATGTATTTGAATTAAATGAAGCATTTGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGCATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCAAAGTATAAATAATAAGAAAACA  >  CP000730/406821‑407095                                                                                                                                             |                                                                                                                                       aGAATTAAGCGATATTGATGTATTTGAATTAAATGAAGCATTTGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGt                                                                                                                                 <  1:357984/148‑1 (MQ=255) aGAATTAAGCGATATTGATGTATTTGAATTAAATGAAGCATTTGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGt                                                                                                                                 <  1:343120/148‑1 (MQ=255)     ttAAGCGATATTGATGTATTTGAATTAAATGAAGCATTTGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTa                                                                                                                             <  2:242278/148‑1 (MQ=255)           gATATTGATGTATTTGAATTAAATGAAGCATTTGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTAt                                                                                                                       <  1:41860/148‑1 (MQ=255)           gATATTGATGTATTTGAATTAAATGAAGCATTTGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTAt                                                                                                                       <  1:297439/148‑1 (MQ=255)           gATATTGATGTATTTGAATTAAATGAAGCATTTGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTAt                                                                                                                       <  1:95164/148‑1 (MQ=255)                    gTATTTGAATTAAATGAAGCATTTGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACat                                                                                                              >  1:14416/1‑148 (MQ=255)                              tAAATGAAGCATTTGCAGCACAATCATTAGCTGTTTATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGTTATTGCATTAGGACATCCTATTGGTGCATCTGGTGCTAGAGTATTAGTGACATTATTGCATCa                                                                                                    <  2:381474/148‑1 (MQ=255)                                  tGAAGCATTTGCAGCGCAATCATTAGCTGTGGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTg                                                                                                <  2:59574/148‑1 (MQ=255)                                  tGAAGCATTTGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTg                                                                                                <  1:376048/148‑1 (MQ=255)                                  tGAAGCATTTGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTg                                                                                                >  2:307514/1‑148 (MQ=255)                                           tGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGaa                                                                                       >  1:179122/1‑148 (MQ=255)                                           tGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGaa                                                                                       >  1:167733/1‑148 (MQ=255)                                                  caATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGaaa                                                                                >  1:214367/1‑148 (MQ=255)                                                  caATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGaaa                                                                                >  2:309985/1‑148 (MQ=255)                                                     tCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTg                                                                             >  1:221256/1‑148 (MQ=255)                                                      cATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTgg                                                                            <  2:159289/148‑1 (MQ=255)                                                       aTTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGt                                                                           >  2:376414/1‑148 (MQ=255)                                                       aTTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGt                                                                           >  2:274203/1‑148 (MQ=255)                                                       aTTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGt                                                                           >  2:236665/1‑148 (MQ=255)                                                       aTTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCAGTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGt                                                                           >  1:7485/1‑148 (MQ=255)                                                          aGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTa                                                                        <  2:199902/148‑1 (MQ=255)                                                          aGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTa                                                                        <  2:6302/148‑1 (MQ=255)                                                          aGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTa                                                                        <  2:6747/148‑1 (MQ=255)                                                          aGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTa                                                                        <  2:6366/148‑1 (MQ=255)                                                             tGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAAca                                                                     >  1:317052/1‑148 (MQ=255)                                                                       gAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTa                                                           <  1:9852/148‑1 (MQ=255)                                                                            aaaaTTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATtggt                                                      <  2:144931/148‑1 (MQ=255)                                                                            aaaaTTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATtggt                                                      <  2:217443/148‑1 (MQ=255)                                                                               aTTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGc                                                   >  1:334186/1‑148 (MQ=255)                                                                               aTTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTACATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGc                                                   >  1:254150/1‑148 (MQ=255)                                                                                ttACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCg                                                  <  2:169363/148‑1 (MQ=255)                                                                                ttACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCg                                                  <  2:112168/148‑1 (MQ=255)                                                                                   cctcctGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTc                                               >  1:108234/1‑148 (MQ=255)                                                                                        tGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCt                                          >  2:355783/1‑148 (MQ=255)                                                                                        tGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCt                                          >  2:300825/1‑148 (MQ=255)                                                                                        tGAAAAGGGGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCt                                          >  2:20430/1‑148 (MQ=255)                                                                                           aaaGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATc                                       >  2:359766/1‑148 (MQ=255)                                                                                           aaaGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATc                                       <  2:317052/148‑1 (MQ=255)                                                                                                  aaTGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAg                                <  2:9028/148‑1 (MQ=255)                                                                                                  aaTGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAg                                <  2:376736/148‑1 (MQ=255)                                                                                                          aGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCa                        <  1:299473/148‑1 (MQ=255)                                                                                                          aGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCa                        <  2:327733/148‑1 (MQ=255)                                                                                                          aGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCa                        <  2:185175/148‑1 (MQ=255)                                                                                                          aGGTGGCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCa                        >  2:141095/1‑148 (MQ=255)                                                                                                          aGGTGGCGCTATTGCATTAGGACATCCTATGGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCa                        >  2:89325/1‑148 (MQ=255)                                                                                                              ggCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCAAAGt                    <  1:225238/148‑1 (MQ=255)                                                                                                              ggCGCTATTGCATTAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCAAAGt                    <  2:244759/148‑1 (MQ=255)                                                                                                                           tAGGACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCAAAGTATAAATAATAAGa       >  1:1981/1‑148 (MQ=255)                                                                                                                              gACATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCAAAGTATAAATAATAAGaaaa    >  2:100491/1‑148 (MQ=255)                                                                                                                               aCATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCAAAGTATAAATAATAAGAAAAc   >  1:66889/1‑148 (MQ=255)                                                                                                                                cATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCAAAGTATAAATAATAAGAAAACa  >  2:60778/1‑148 (MQ=255)                                                                                                                                cATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCAAAGTATAAATAATAAGAAAACa  >  1:137796/1‑148 (MQ=255)                                                                                                                                cATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCAAAGTATAAATAATAAGAAAACa  >  2:230782/1‑148 (MQ=255)                                                                                                                                cATCCTATTGGTGGATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCAAAGTATAAATAATAAGAAAACa  >  2:130122/1‑148 (MQ=255)                                                                                                                                             |                                                                                                                                       AGAATTAAGCGATATTGATGTATTTGAATTAAATGAAGCATTTGCAGCACAATCATTAGCTGTTGATCGTGAATTAAAATTACCTCCTGAAAAGGTGAATGTTAAAGGTGGCGCTATTGCATTAGGACATCCTATTGGTGCATCTGGTGCTAGAGTATTAGTGACATTATTGCATCAACTGAATGATGAAGTTGAAACTGGTTTAACATCATTGTGTATTGGTGGCGGTCAAGCTATCGCTGCAGTTGTATCAAAGTATAAATAATAAGAAAACA  >  CP000730/406821‑407095

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