Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	801,869	G→A	7.8%	A118A (GCG→GCA)	obg →	GTPase Obg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	801,869	0	G	A	7.8%	144.6 / 2.2	51	A118A (GCG→GCA)	obg	GTPase Obg
Reads supporting (aligned to +/- strand):  ref base G (21/26);  new base A (2/2);  total (23/28)
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.81e-01

CTCCAACGGCGGTAGCACCGACTGCACCGGCAAGAAGGGTGAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGAAGGTGTTGGCCGACGTCGGCCTGCTGGGCTTG  >  NC_002947/801725‑802016                                                                                                                                                 |                                                                                                                                                    cTCCAACGGCGGTAGCACCGACTGCACCGGCAAGAAGGGTGAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAg                                                                                                                                                 <  6:63721/149‑1 (MQ=255)      aCGGCGGTAGCACCGACTGCACCGGCAAGAAGGGTGAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGtggc                                                                                                                                            <  1:262839/149‑1 (MQ=255)        ggcggTAGGACCGACTGCTCCGGCAAGAAGGGTGATGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGtggctg                                                                                                                                          <  8:301927/149‑1 (MQ=255)         gcggTAGCACCGACTGCACCGGCAAGAAGGGTGAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGtggctgg                                                                                                                                         <  6:296036/149‑1 (MQ=255)           ggTAGCACCGACTGCACCGGCAAGAAGGGTGAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCACGGTGGCTGGCa                                                                                                                                       >  2:72708/1‑149 (MQ=255)                  ccGACTGCACCGTCAAGAAGGGTGAAGACCTGTTCCTGCGCGTGCCGGTCGTCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGTCCAGATGATGATGGTCGAGCAATGTGGCTGGCACGGGCtt                                                                                                                                <  3:203557/149‑1 (MQ=255)                  ccGACTGCACCGGCAAGAAGGGTGAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACCGGCtt                                                                                                                                >  5:90949/1‑149 (MQ=255)                         cACCGGCAAGAAGGGTGAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCaa                                                                                                                           <  7:25785/147‑1 (MQ=255)                          aCCGGCAAGAAGGGTGAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAacac                                                                                                                        >  3:258733/1‑149 (MQ=255)                                   aagGGTGAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCaa                                                                                                               <  2:57196/149‑1 (MQ=255)                                      ggTGAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTc                                                                                                            >  6:179284/1‑149 (MQ=255)                                         gAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGc                                                                                                                             <  8:30845/129‑1 (MQ=255)                                         gAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGc                                                                                                                             >  7:30845/1‑129 (MQ=255)                                           aGACCTGTTCCTGCGCGTGCAGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCa                                                                                                       <  1:247642/149‑1 (MQ=255)                                               cTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCaa                                                                                                   >  6:124290/1‑149 (MQ=255)                                               cTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCaa                                                                                                   >  4:189873/1‑149 (MQ=255)                                                 gTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAAcc                                                                                                 >  3:33409/1‑149 (MQ=255)                                                         gcgTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTgcgccg                                                                                         >  2:229743/1‑149 (MQ=255)                                                               cGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCAAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACAAAACGTgagcagagc                                                                                      >  4:40012/1‑141 (MQ=255)                                                                ggTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGa                                                                                  <  4:33409/149‑1 (MQ=255)                                                                  tCGGCTCCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGac                                                                                 >  3:47132/1‑148 (MQ=255)                                                                  tCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCa                                                                                    >  7:236549/1‑145 (MQ=255)                                                                  tCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCa                                                                                    <  8:236549/145‑1 (MQ=255)                                                                           ccGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCg                                                                                                                     >  7:291359/1‑103 (MQ=255)                                                                           ccGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCg                                                                                                                     <  8:291359/103‑1 (MQ=255)                                                                           ccGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGc                                                                       <  2:141673/149‑1 (MQ=255)                                                                                aTCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAAcc                                                                  >  6:243600/1‑149 (MQ=255)                                                                                   gATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCgg                                                                <  5:129133/148‑1 (MQ=255)                                                                                   gATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCgg                                                                <  5:179284/148‑1 (MQ=255)                                                                                    aTGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGt                                                              <  4:258733/149‑1 (MQ=255)                                                                                        ttCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGagc                                                          <  4:47132/149‑1 (MQ=255)                                                                                        ttCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGag                                                           >  3:146298/1‑148 (MQ=255)                                                                                        ttCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCACAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGag                                                           >  5:217867/1‑148 (MQ=255)                                                                                        ttCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCACAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGag                                                           >  5:213101/1‑148 (MQ=255)                                                                                          ccACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGagca                                                         <  3:189873/148‑1 (MQ=255)                                                                                           cACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGTGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGagcag                                                        <  5:324218/148‑1 (MQ=255)                                                                                            aCGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCg                                                      >  4:130447/1‑149 (MQ=255)                                                                                                            gTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGaa                                      <  5:294856/149‑1 (MQ=255)                                                                                                                    gTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGAAGgt                              >  3:148402/1‑149 (MQ=255)                                                                                                                    gTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGAAGgt                              >  7:29475/1‑149 (MQ=255)                                                                                                                    gTCACCCCCGGACAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCTCCGAGCCAGACCACACCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGaat                                <  3:40012/147‑2 (MQ=255)                                                                                                                       aCCCCCGGCCAGAAGCTGATTGTCGCGCAAGGTGTCTGGCACGGGCTTGGCAACACCCTTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATTAATGTGtt                           <  6:90949/149‑1 (MQ=255)                                                                                                                        cccccGGCCAGAAGCTGATGGTCGCACAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGAAGGTGTTg                          <  6:213101/149‑1 (MQ=255)                                                                                                                        cccccGGCCAGAAGCTGATGGTCGCACAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGAAGGTGTTg                          <  6:217867/149‑1 (MQ=255)                                                                                                                         ccccGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGa                                  <  5:288448/140‑1 (MQ=255)                                                                                                                         ccccGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGa                                  >  6:288448/1‑140 (MQ=255)                                                                                                                           ccGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGAAGGTGTTGGcc                       >  5:153034/1‑149 (MQ=255)                                                                                                                           ccGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGAAGGTGTTGGc                        >  1:227337/1‑148 (MQ=255)                                                                                                                           ccGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGAAGGTGTTGGc                        <  2:227337/148‑1 (MQ=255)                                                                                                                                    aGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGAAGGTGTTGGCCGACGTCGGc              <  1:229743/149‑1 (MQ=255)                                                                                                                                         aTGGTCGCGCAAGGGGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGAAGGTGTTGGCCGACGTCGGCctg           <  1:72708/147‑1 (MQ=255)                                                                                                                                                cgcAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGAAGGTGTTGGCCGACGTCGGCCTGCTGGGCTTg  <  5:124290/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CTCCAACGGCGGTAGCACCGACTGCACCGGCAAGAAGGGTGAAGACCTGTTCCTGCGCGTGCCGGTCGGCACCACCGTGATCGATGCTTCCACGCAGGAAGTGATCGGTGACCTGGTCACCCCCGGCCAGAAGCTGATGGTCGCGCAAGGTGGCTGGCACGGGCTTGGCAACACCCGTTTCAAGTCGAGCACCAACCGTGCGCCGCGCCAGACCACCCCAGGCAAACCGGGTGAGCAGCGTGACCTGAAGATGGAAATGAAGGTGTTGGCCGACGTCGGCCTGCTGGGCTTG  >  NC_002947/801725‑802016

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