Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,021,487	G→C	9.5%	intergenic (+24/‑415)	PP_1798 → / → gmd	outer membrane efflux protein/GDP‑mannose 4,6‑dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,021,487	0	G	C	9.5%	118.4 / 4.2	42	intergenic (+24/‑415)	PP_1798/gmd	outer membrane efflux protein/GDP‑mannose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base G (18/20);  new base C (2/2);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

CTCAACCGTTGGGCGGGTGAGTTGAATCAGGAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGCTTCATCGACCTTAGCGCAGGCCGACCCGGGCGGTGGCTGC  >  NC_002947/2021341‑2021633                                                                                                                                                   |                                                                                                                                                      cTCAACCGTTGGGCGGGTGAGTTGAATCAGGAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAgggg                                                                                                                                                     >  6:350270/1‑149 (MQ=255) cTCAACCGTTGGGCGGGTGAGTTGAATCAGGAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAgg                                                                                                                                                       <  1:236320/147‑1 (MQ=255) cTCAACCGTTGGGCGGGTGAGTTGAATCAGGAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTACCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAgg                                                                                                                                                       <  5:131698/147‑1 (MQ=255)    aaCGGTTGGGCGGGTGAGTTGAATCAGGAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCGTGAGAGCAGCCTGCCCGGTGGGGAGGTTGCTCGCCTGGGCGCCAGTTGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCt                                                                                                                                                   <  1:224949/148‑1 (MQ=255)     aCCGTTGGGCGGGTGAGTTGAATCAGGAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGa                                                                                                                                                 <  8:217953/149‑1 (MQ=255)           gggCGGGTGAGTTGAATCAGGAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGCCCCTGAagcgc                                                                                                                                            <  1:179662/148‑1 (MQ=255)                   gAGTTGAATCAGGAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCGCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAAc                                                                                                                                   <  2:322307/149‑1 (MQ=255)                      ttGAATCAGGAAAGTCTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCATGCCCGGTGGCGAGGTTCCTCGCCTCGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAg                                                                                                                                >  6:206212/1‑149 (MQ=255)                          aTCAGGAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGaa                                                                                                                            <  2:287548/149‑1 (MQ=255)                             aGGAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAgc                                                                                                                          <  1:104473/148‑1 (MQ=255)                              ggAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAgcg                                                                                                                         >  6:219294/1‑148 (MQ=255)                              ggAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAgcg                                                                                                                         >  1:190737/1‑148 (MQ=255)                                            aCGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAAttt                                                                                                           <  2:190737/148‑1 (MQ=255)                                                            gAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCtttt                                                                                           <  2:106132/148‑1 (MQ=255)                                                            gAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCtttt                                                                                           <  3:218688/148‑1 (MQ=255)                                                            gAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGCGGGTCTGCGGTTCAGGGGCTCAAGCGCAGCGCCAACCAGTGAAGCCCACCGAATCAATTTAGCCCCGCTCGCTTTTg                                                                                          >  1:5334/1‑149 (MQ=255)                                                                ccGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTc                                                                                      <  5:219294/149‑1 (MQ=255)                                                                   gTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTg                                                                                   >  7:173218/1‑149 (MQ=255)                                                                       cTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCtttt                                                                                           <  3:23826/137‑1 (MQ=255)                                                                       cTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCtttt                                                                                           >  4:23826/1‑137 (MQ=255)                                                                        tGAGACCAGCCTGCCCGGTGGCGAGGTTCCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTcc                                                                               >  1:204055/1‑148 (MQ=255)                                                                          agaCCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGccc                                                                                                      >  7:169002/1‑123 (MQ=255)                                                                          agaCCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGccc                                                                                                      <  8:169002/123‑1 (MQ=255)                                                                            aCCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTg                                                                          <  5:350270/149‑1 (MQ=255)                                                                              ctgCCTGCCCGGTGGCGACGTTGCTGGGCTGGGCGGCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGc                                                                         <  2:5334/146‑1 (MQ=255)                                                                              cAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCg                                                                        >  3:114027/1‑149 (MQ=255)                                                                                     cccGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAgg                                                                 >  1:146140/1‑149 (MQ=255)                                                                                      ccGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGc                                                                <  4:296574/149‑1 (MQ=255)                                                                                       cGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑g                                                            <  5:206212/149‑1 (MQ=255)                                                                                                        gCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGttt                                             <  2:146140/147‑1 (MQ=255)                                                                                                         ccTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGc                                          <  2:204055/149‑1 (MQ=255)                                                                                                              gcgcCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGCTTCAt                                     >  2:111788/1‑149 (MQ=255)                                                                                                                 ccAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTc                                                                                      <  3:307274/100‑1 (MQ=255)                                                                                                                 ccAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTc                                                                                      >  4:307274/1‑100 (MQ=255)                                                                                                                        cTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGCTTCATCGACCTTAgc                           >  3:288437/1‑149 (MQ=255)                                                                                                                         tGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGCTTCATCGACCTTAgcg                          >  3:239739/1‑149 (MQ=255)                                                                                                                            aaGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGCTTCATCGACCTTAGCGCAg                       <  4:7377/149‑1 (MQ=255)                                                                                                                                 aGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGCTTCATCGACCTTAGCGCAGGCCGa                  >  5:289626/1‑149 (MQ=255)                                                                                                                                      gTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGtt                                              <  2:3884/116‑1 (MQ=255)                                                                                                                                      gTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGtt                                              >  1:3884/1‑116 (MQ=255)                                                                                                                                         tGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGCTTCATCGACCTTAGCGCAGGCCGACCCGGGCg          >  4:66450/1‑149 (MQ=255)                                                                                                                                         tGCGCTTCAGCCCCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGCGTGGCGGAGAAGATGGTTTGGCTTCATCGACCTTAGCGCAGGCCGACCCgg             >  1:119646/1‑149 (MQ=255)                                                                                                                                         tGCGCTTCAGCCCCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGCTTCATc                                    >  5:317578/1‑123 (MQ=255)                                                                                                                                         tGCGCTTCAGCCCCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGCTTCATc                                    <  6:317578/123‑1 (MQ=255)                                                                                                                                                cAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGCTTCATCGACCTTAGCGCAGGCCGACCCGGGCGGTGgctg   >  4:347248/1‑149 (MQ=255)                                                                                                                                                  ggggCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGCTTCATCGACCTTAGCGCAGGCCGACCCGGGCGGTGgctgc  >  4:163569/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CTCAACCGTTGGGCGGGTGAGTTGAATCAGGAAAGTGTGGATAACGTCAATGTCTGGTTGAGCCCGGTGCCTGAGACCAGCCTGCCCGGTGGCGAGGTTGCTCGCCTGGGCGCCAGTCGCTGAAAGCGAGGGGGTCTGCGCTTCAGGGGCTGAAGCGCAGCGCCAACCAGTGAAGCGCACAGTATCAATTTAGCCCCGGTCCCTTTTGCGTCCTGTTCCCATTGCGTGTCAGGC‑‑‑GAGGTCAAGATGGTTTGGCTTCATCGACCTTAGCGCAGGCCGACCCGGGCGGTGGCTGC  >  NC_002947/2021341‑2021633

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