Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	445,105	A→T	14.9%	intergenic (+78/‑210)	PP_0367 → / → PP_0368	hypothetical protein/acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	445,105	0	A	T	14.9%	110.3 / 11.9	47	intergenic (+78/‑210)	PP_0367/PP_0368	hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (22/18);  new base T (4/3);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

GAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAAC  >  NC_002947/444959‑445252                                                                                                                                                   |                                                                                                                                                     gAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGCCACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCATGCCCGCGGATGCGTCCGCTCAGCcaacaaat                                                                                                                                                    >  5:280021/1‑147 (MQ=255)   gCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                  <  6:151498/149‑1 (MQ=255)   gCTCGGCGCAAAGGTCGTCAACGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                  <  6:15610/149‑1 (MQ=255)       ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                              <  8:207338/149‑1 (MQ=255)       ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCt                                                                                                                                              >  6:302220/1‑149 (MQ=37)         cgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGGCCCCGAATGAGTCCGCTCAGACAACAAATTTGTGTgg                                                                                                                                            >  2:279571/1‑149 (MQ=255)         cgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCAGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTgg                                                                                                                                            <  5:183617/149‑1 (MQ=255)          gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTgat                                                                                                                                           <  7:60047/149‑3 (MQ=25)             aaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                                        <  7:292475/149‑1 (MQ=255)                  cgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctc                                                                                                                                   <  5:131458/149‑1 (MQ=255)                   gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctct                                                                                                                                  <  1:334800/149‑1 (MQ=255)                        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGGCCGCGAATGCGTCCCATAAGAAAACAATTTTTTATGGTGCGGCCTCTTagcg                                                                                                                             >  3:301257/1‑149 (MQ=255)                         gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAATATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                            <  1:168365/149‑1 (MQ=18)                         gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                            <  5:302220/149‑1 (MQ=18)                           aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcg                                                                                                                             <  8:121305/146‑1 (MQ=255)                           aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcg                                                                                                                             >  7:121305/1‑146 (MQ=255)                             ccGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATGCGTCCGATCAGACAACAATTTTGTGTGGTCGGGCCTCTTAGCGGG‑Caa                                                                                                                       >  5:375033/1‑149 (MQ=255)                                  gaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCg                                                                                                                  <  2:332510/149‑1 (MQ=255)                                    aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCATTCCCCCGGAATCGTGCCCTCACACAACAAATTTTTAAGGTGGGGCCACTTAGCCGG‑CCAGcccccg                                                                                                                >  7:2557/1‑146 (MQ=255)                                       gCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCCCGAATTCGTCCGATAAGACAAAAATTTTTTCTGGGCGGGCCGATTAGCGGGCCAAGccgcg                                                                                                                 >  8:152349/1‑143 (MQ=255)                                             ccGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCCGCGCCTGCCTTGTCGGCGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAgg                                                                                                       <  1:258129/149‑1 (MQ=255)                                                                    gAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTTACCCCGAAAACGGACCAAAAGACAAAAATTTTTTGTGGTAGGGGCGATAAGCGGG‑CAAAGCCCCGCCAACAAAGGTTTGTGgcatagcatgcaagat                                                                                >  4:86640/1‑143 (MQ=255)                                                                        aCTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                                                 >  2:144749/1‑97 (MQ=255)                                                                        aCTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                                                 <  1:144749/97‑1 (MQ=255)                                                                                 aTCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATGCGTCCGATCAGACAACAATTTTTTCTGGTCGGGCCTCTTAGCGGG‑CAAGCCCGCTCCGACAACGGTTTGTGGCATGGCATGACTGATTTGCCAtttttag                                                                   >  7:93497/1‑147 (MQ=255)                                                                                  tCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGCGTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CTAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGt                                                                  <  5:278511/149‑1 (MQ=255)                                                                                    cccTATTTCCGCCCCGGCCTTGTAGGAGCGGGCTTGCCCGCTAATGCGTCCTATCACACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTg                                                                 <  4:50832/148‑1 (MQ=38)                                                                                                gcctgcctTGTAGGAGCGGGCTGGCCCGCTAATGCGTCCTATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCAc                                                    <  1:279571/149‑1 (MQ=255)                                                                                                            ttcgCGGGCAAGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                                        >  7:246133/4‑149 (MQ=255)                                                                                                                cGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                    <  7:251288/149‑1 (MQ=255)                                                                                                                   gCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                    >  7:89959/1‑146 (MQ=255)                                                                                                                        cccGCTAATGCGTCCGTTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc                            <  3:86640/149‑1 (MQ=255)                                                                                                                        cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc                            >  5:212092/1‑149 (MQ=18)                                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGAc                                                                                    >  2:4866/1‑92 (MQ=255)                                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                           >  1:278994/1‑149 (MQ=255)                                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                           >  2:243089/1‑149 (MQ=255)                                                                                                                          cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACaa                                                                                                         <  7:232935/70‑1 (MQ=255)                                                                                                                          cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACaa                                                                                                         >  8:232935/1‑70 (MQ=255)                                                                                                                          cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                           >  2:13328/1‑148 (MQ=255)                                                                                                                           gcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                           >  4:330305/1‑147 (MQ=255)                                                                                                                              aaTGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg                        <  1:229941/147‑1 (MQ=255)                                                                                                                              aaTGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg                        >  2:229941/1‑147 (MQ=255)                                                                                                                                   gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACg                 >  2:172231/1‑149 (MQ=255)                                                                                                                                   gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACg                 >  5:159226/1‑149 (MQ=255)                                                                                                                                    tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACg                 >  3:320823/1‑148 (MQ=255)                                                                                                                                     ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                >  5:316564/1‑148 (MQ=38)                                                                                                                                          cAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGttt          >  2:61283/1‑149 (MQ=37)                                                                                                                                             acaacaATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCa        >  3:152137/1‑148 (MQ=255)                                                                                                                                              caacaaTTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa      <  6:159226/149‑1 (MQ=255)                                                                                                                                               aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCGCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACACAAAGCTGAGATAAACTTAACTTTCaaa      >  8:288766/1‑148 (MQ=255)                                                                                                                                               aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc         >  5:21619/1‑145 (MQ=255)                                                                                                                                               aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc         <  6:21619/145‑1 (MQ=255)                                                                                                                                               aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa      >  3:134284/1‑148 (MQ=255)                                                                                                                                                  aaTTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac  <  8:93497/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGG‑CAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAAC  >  NC_002947/444959‑445252

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