Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,436,713	C→T	8.9%	intergenic (+10/‑128)	yfbS → / → PP_3932	putative transporter/diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,436,713	0	C	T	8.9%	120.7 / 2.6	45	intergenic (+10/‑128)	yfbS/PP_3932	putative transporter/diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (24/17);  new base T (2/2);  total (26/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 = 5.36e-01

GTCATCACCGGTCAACACGCTGGTGCTTGGGCCCGGGCAATACCGCTTCGCTGACTTCGTCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGATTTCGTTTCCGCTTGCCGGTCGAGTTCATGGT  >  NC_002947/4436568‑4436860                                                                                                                                                  |                                                                                                                                                     ttcatcaCCGGGCAACACGGGGGTGGGTGGGCCCGGGCAATCCCGCTTCGCTGACTTCGTCAACTTTGGCGTGCCGTTCACCGTGCTTGTGATGCTGGTTCGCGGCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAACTGGCCAg                                                                                                                                                   <  6:266803/148‑1 (MQ=255)           gTCAACACGCTGGTGCTTGGGCCCGGGCAATACCGCTTCGCTGACTTCGTCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAgg                                                                                                                                          >  1:328102/1‑148 (MQ=255)               acacGCTGGTGCTTGGGCCCGGGCAATACCGCTTCGCTGACTTCGTCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAAtt                                                                                                                                     >  7:71914/1‑149 (MQ=255)                                 ccGGGCAATACCGCTTCGCTGACTTCGTCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCaa                                                                                                                   >  3:60398/1‑149 (MQ=255)                                  cGGGCAATACCGCTTCGCTGACTTCGTCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAAt                                                                                                                  <  8:267281/149‑1 (MQ=255)                                       aaTACCGCTTCGCTGACTTCGTCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTTACATTATCGCAACCAATCGCgg                                                                                                             >  8:268482/1‑149 (MQ=255)                                            cGCTTCGCTGACTTCGTCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGAt                                                                                                        >  6:236549/1‑149 (MQ=255)                                              cTTCGCTGACTTCGTCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGAtt                                                                                                       >  2:182178/1‑148 (MQ=255)                                                      aCTTCGTCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAac                                                                                              <  5:236549/149‑1 (MQ=255)                                                        ttCGTCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAacca                                                                                            >  3:262006/1‑149 (MQ=255)                                                            tCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGAtt                                                                                                       >  8:123763/1‑134 (MQ=255)                                                            tCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGAtt                                                                                                       <  7:123763/134‑1 (MQ=255)                                                            tCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCt                                                                                        <  4:60398/149‑1 (MQ=255)                                                              aaGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGa                                                                                      <  8:304895/149‑1 (MQ=255)                                                               aGATTGGCGTGCCGTTCACGGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTTGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAg                                                                                     <  7:150598/149‑1 (MQ=255)                                                                gATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGc                                                                                    >  2:109420/1‑149 (MQ=255)                                                                     gCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAac                                                                                <  4:262006/148‑1 (MQ=255)                                                                       gTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAaca‑ca                                                                            >  7:79163/1‑149 (MQ=255)                                                                                aCCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGTCAGAAATAGAGGCAATTTGACATTATc                                                                                                                           >  4:297103/1‑94 (MQ=255)                                                                                aCCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGTCAGAAATAGAGGCAATTTGACATTATc                                                                                                                           >  4:301843/1‑94 (MQ=255)                                                                                aCCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGTCAGAAATAGAGGCAATTTGACATTATc                                                                                                                           <  3:301843/94‑1 (MQ=255)                                                                                aCCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGTCAGAAATAGAGGCAATTTGACATTATc                                                                                                                           <  3:297103/94‑1 (MQ=255)                                                                                  cGTGCTGGCGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGAAATCAAACACCCTGAGCCAACA‑CACTGCAGGTaaa                                                                 >  5:219885/1‑148 (MQ=255)                                                                                          tGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAg                                                         <  8:71914/149‑1 (MQ=255)                                                                                            aTGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGAATGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCACCAgtt                                                       >  3:236719/1‑147 (MQ=255)                                                                                             tGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCaa                                                                                                                       >  6:63739/1‑85 (MQ=255)                                                                                             tGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCaa                                                                                                                       <  5:63739/85‑1 (MQ=255)                                                                                                  gTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACACCACTGCAGGTAATTGCCCCAGGTGctct                                                  <  1:182178/149‑1 (MQ=255)                                                                                                         tcgtGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGtt                                          >  6:49003/1‑149 (MQ=255)                                                                                                               tGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGGAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCt                                     <  8:79163/148‑1 (MQ=255)                                                                                                                  tGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGAt                                 >  1:167847/1‑149 (MQ=255)                                                                                                                     cGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGATTTc                              <  2:328102/149‑1 (MQ=255)                                                                                                                         gTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGATTTCGtt                           >  4:240378/1‑148 (MQ=255)                                                                                                                          ttcttcGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGATTTCGTTTc                         >  2:204877/1‑149 (MQ=255)                                                                                                                             ttcGTGCTGTGATGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGTTGCTCTATCGCGTTg                                         <  6:265624/130‑1 (MQ=255)                                                                                                                             ttcGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCTCAACCAATAGCGGCAGATTGACATCAACAACACTGAGCAAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATAGCGTTg                                         >  5:265624/1‑130 (MQ=255)                                                                                                                                   cgttctggAAAGTGGCCAGATATAGAGGCAATTTGTCATTATCGCTACCAATCGCGGCTGATTGACATCAACCACCCTGAGCCAACA‑CAGTGCAGGTAATTGCCCCAGGTGCTGTATCGCGTTGTCCTGGATTTCGTTTCCGCTTGcc                 <  4:72760/142‑1 (MQ=255)                                                                                                                                     gtgAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGATTTCGTTTCCGCTTGCCGGt              >  4:135064/1‑149 (MQ=255)                                                                                                                                       gAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGATTTCGTTTCCGCTTGCCGGTc             >  1:313479/1‑148 (MQ=255)                                                                                                                                        aGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGATTTCGTTTCCGCTTGCCGGTCGa           <  7:268482/149‑1 (MQ=255)                                                                                                                                          gAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATAAACCACCCTGAGCCACCA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGATTTCGTTTCCGCTTGcc                 >  3:72760/1‑141 (MQ=255)                                                                                                                                            aaGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGATTTCGTTTCCGCTTGCCGGTCGAGtt        >  1:18996/1‑148 (MQ=255)                                                                                                                                                ggCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGATTTCGTTTCCGCTTGCCGGTCGAGTTCATg    <  3:135064/148‑1 (MQ=255)                                                                                                                                                 gCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCACCA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCGGGATTTCGTTTCGTCTTGCCGGTCGAGTTCATgga  >  5:271720/1‑148 (MQ=255)                                                                                                                                                 gCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGATTTCGTTTCCGCTTGCCGGTCGAGTTCATgtt  >  7:127426/1‑147 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GTCATCACCGGTCAACACGCTGGTGCTTGGGCCCGGGCAATACCGCTTCGCTGACTTCGTCAAGATTGGCGTGCCGTTCACCGTGCTGGTGATGCTGGTTACCGTCGTGATGGTGCCGTGGTTCTTCGGGCTGTGAGGAAAGTGGCCAGAAATAGAGGCAATTTGACATTATCGCAACCAATCGCGGCAGATTGACATCAACCACCCTGAGCCAACA‑CACTGCAGGTAATTGCCCCAGGTGCTCTATCGCGTTGTCCTGGATTTCGTTTCCGCTTGCCGGTCGAGTTCATGGT  >  NC_002947/4436568‑4436860

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