Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,338,494	A→T	11.3%	intergenic (+114/+28)	PP_3812 → / ← PP_5632	lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,338,494	0	A	T	11.3%	119.5 / 6.8	53	intergenic (+114/+28)	PP_3812/PP_5632	lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (19/28);  new base T (4/2);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

GCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATC  >  NC_002947/4338348‑4338638                                                                                                                                                      |                                                                                                                                                 gCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAAc                                                                                                                                                <  4:99947/149‑1 (MQ=255) gCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAAc                                                                                                                                                <  4:105356/149‑1 (MQ=255)   gggCAACGAACTTAAAGTGCATCTGCAGTGAGTCGATGAGCAATAGACCCGCATGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGGGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTTTAccc                                                                                                                                              <  8:343326/149‑1 (MQ=255)     gCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCcgc                                                                                                                                            <  4:223702/149‑1 (MQ=255)       aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑G‑‑GGGGGTACCcgcg                                                                                                                                           >  4:59835/1‑148 (MQ=255)         cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCAGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGaaa                                                                                                                                        <  1:166866/149‑1 (MQ=255)          gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGGCCCTTAG‑CCGGGGTAAGCCCGCGac                                                                                                                                         >  7:128933/1‑148 (MQ=255)            aCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAggg                                                                                                                                     <  4:187101/149‑1 (MQ=255)              ttAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCCTTAG‑CG‑GGGTAAACCCCGCAAAGGGc                                                                                                                                    >  1:369511/1‑149 (MQ=255)              ttAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTAG‑CG‑GGGTAAACCCCCAAAAggcc                                                                                                                                    >  3:211637/1‑147 (MQ=255)                aaaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAACCCCGCGAAAGGGcc                                                                                                                                   <  7:256097/147‑1 (MQ=255)                aaaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGcc                                                                                                                                   <  4:253071/147‑1 (MQ=255)                  aGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑CGGTAAACCCGCGAAAGGGCCGGc                                                                                                                                >  7:176076/1‑148 (MQ=255)                    tGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCaaa                                                                                                                             >  8:96379/1‑149 (MQ=255)                      cATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAAt                                                                                                                            <  2:272435/148‑1 (MQ=255)                       aTCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGcc                                                                                                                                   <  3:273129/140‑1 (MQ=255)                       aTCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTAG‑A‑‑GGGTTAACCCGCGAAAGGGCCgcgctcgg                                                                                                                           >  4:273129/1‑141 (MQ=255)                            cAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCt                                                                                                                      >  3:47812/1‑148 (MQ=255)                            cAGGGAGGCGATGAGCAATACACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGCAAACAGCCCGCAGGCTTCAATCAATCTGCCAGCGGGGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCt                                                                                                                         <  3:59835/145‑1 (MQ=255)                              ggcagccGATCAGCAATACACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGCAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                                                                                                   <  7:316826/143‑1 (MQ=255)                                         aGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAccc                                                                                                        >  6:324217/1‑149 (MQ=255)                                                gACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCcaca                                                                                                 <  7:96379/149‑1 (MQ=255)                                                 aCCCGCCTGGGCGGGTCTTTTGCTTCTGTTTGCGGAAACAGCCCGCAGGCTTCAATGAATCAGCCAGCTGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGGTCAAACCCTTCCACaa                                                                                                >  5:270459/1‑149 (MQ=255)                                                   ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCCTTAGCC‑‑GGGTAAACCCCCCCAAGGCCCGCCAAAACAT‑ATCAGGGGTCAAAaccagacaaa                                                                                                 >  4:6071/1‑139 (MQ=255)                                                             gggTATTTTTCTTATATTTGCTTAAATACCCCGCAGGATTAAATTAATCTCCCTGCGGTGCGGCTGGTCCGGCCCTTTCG‑C‑‑GGTAAACCCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACaa                                                                                                <  8:328275/137‑1 (MQ=255)                                                                     tggttatATTTGCGGAAACAGCCCGCAGGTTTCAATGAATCTGCCTGCGGGGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGTAAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGc                                                                             <  3:6071/142‑1 (MQ=255)                                                                        ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCG‑C‑‑GGGTAAACCCGCGAAAGGGGCGGCAAATCCTGCTCAGGGCGGAAACCCTTCCCCAATGATCACCTCACCCTTAGCCAcc                                                                         >  7:158913/1‑149 (MQ=255)                                                                                gCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATc                                                                                                                           <  5:182448/91‑1 (MQ=255)                                                                                gCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATc                                                                                                                           >  6:182448/1‑91 (MQ=255)                                                                                 cGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGa                                                                <  4:98863/149‑1 (MQ=255)                                                                                  ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCGCAAACCCTTCCACAATGATCACCTCAGCCTTAGACACCCCTTGGCGAt                                                               >  8:363858/1‑149 (MQ=255)                                                                                    aaaCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATgg                                                             <  7:363858/149‑1 (MQ=255)                                                                                        aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGc                                                         <  5:942/149‑1 (MQ=255)                                                                                           ccGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTAG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGGCAATCCAGCTCAGGGCTCAAACCCATGCACAATGATAACCGCCGCCATAGCCAACCCATGGCGATGGCTGCAtt                                                      >  4:304913/1‑148 (MQ=255)                                                                                                  cTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                               <  4:47812/149‑1 (MQ=255)                                                                                                   ttCAATGAATCTGCCTGGGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTgg                                              <  8:158913/149‑1 (MQ=255)                                                                                                         gAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTAG‑C‑‑GGGGTAACCCGCGAAAGGGGCGGCAAATCCTGGTCAGGGCGGAAACCCGTCCACAATGATCACCACCCCCTTCGCCCCCCCGTGGCGCTGGCTGCATGCCGCGTGGTagcgc                                        >  2:249978/1‑145 (MQ=255)                                                                                                               gCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTAAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAg                                  >  5:363953/1‑149 (MQ=255)                                                                                                                ccAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGc                                 <  4:211637/149‑1 (MQ=255)                                                                                                                  atcGGTTCGGGAGGTCCTGCCCGTTCG‑C‑‑GGGTTAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCgg                               <  6:270459/147‑1 (MQ=255)                                                                                                                       tGCGGCCGGTCCGCCCCTTTCG‑C‑‑GGGTAAACCCGCGAACGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACa                          <  8:128933/149‑1 (MQ=255)                                                                                                                       tGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACa                          <  4:171386/149‑1 (MQ=255)                                                                                                                         cGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCa                                                        <  1:330482/117‑1 (MQ=255)                                                                                                                              ttgccGGCCCTTTCG‑C‑‑GGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAAc                           >  6:76245/4‑141 (MQ=255)                                                                                                                              ttgccGGCCCTTTCG‑C‑‑GGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAAc                           <  5:76245/138‑1 (MQ=255)                                                                                                                                tCCGGCCCTTTAG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCGTCCACAATGATCACCACAGCCGTAGCCCCCCCGTGGCGAAGGCTGCATGCCGCATGGTACAGCGCGGGGCGGTAACAGGCGCCGGc                 >  8:143948/1‑149 (MQ=255)                                                                                                                                 ccGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGcc                >  1:4216/1‑149 (MQ=255)                                                                                                                                  cGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCt               >  6:347156/1‑149 (MQ=255)                                                                                                                                       cTTTCG‑C‑‑GGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCCCCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg          >  6:361851/1‑149 (MQ=255)                                                                                                                                        tttCG‑C‑‑GGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg          >  8:279389/1‑148 (MQ=255)                                                                                                                                          tCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa        <  3:304913/148‑1 (MQ=255)                                                                                                                                           cg‑c‑‑gGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTAc       <  4:193306/148‑1 (MQ=255)                                                                                                                                                  ggTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGCTCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATc  <  7:268127/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                 GCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCG‑C‑‑GGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATC  >  NC_002947/4338348‑4338638

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: ‑