Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	399,528	A→G	7.1%	V19A (GTG→GCG) ‡	PP_0332 ←	acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	399,528	0	A	G	7.1%	157.4 / 2.6	56	V19A (GTG→GCG) ‡	PP_0332	acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base G (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

GGATGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGACGG  >  NC_002947/399385‑399666                                                                                                                                                |                                                                                                                                           ggATGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTa                                                                                                                                       <  8:297151/149‑1 (MQ=255)   aTGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTAtt                                                                                                                                     <  4:88123/149‑1 (MQ=255)    tGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTAtt                                                                                                                                     >  4:209498/1‑148 (MQ=255)      cacaGGGCGAAGATGGCGTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGag                                                                                                                                  <  8:277792/149‑1 (MQ=255)       acaGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCGTGGTATTGatc                                                                                                                                 <  2:307899/149‑3 (MQ=255)        caGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCGTGGTATTGatct                                                                                                                                <  5:235871/149‑4 (MQ=255)              gAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGagcagcag                                                                                                                            >  5:276512/1‑147 (MQ=255)              gAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGagcagcag                                                                                                                            >  5:176124/1‑147 (MQ=255)                aGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCa                                                                                                                                            >  1:36169/1‑129 (MQ=255)                aGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCa                                                                                                                                            <  2:36169/129‑1 (MQ=255)                           tCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTga                                                                                                              <  3:274753/148‑1 (MQ=255)                             gATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGaa                                                                                                           >  2:147166/1‑149 (MQ=255)                                   cGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCgg                                                                                                     >  6:256247/1‑149 (MQ=255)                                    gACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGt                                                                                                    >  2:295164/1‑149 (MQ=255)                                      ccAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTca                                                                                                  <  5:29838/149‑1 (MQ=255)                                        aGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTAtt                                                                                                                                     >  3:156950/1‑112 (MQ=255)                                        aGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTAtt                                                                                                                                     <  4:156950/112‑1 (MQ=255)                                          gTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCGACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTcagcag                                                                                              >  5:2894/1‑149 (MQ=255)                                          gTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTcagcag                                                                                              <  1:295164/149‑1 (MQ=255)                                             tCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTcagcag                                                                                              >  1:173523/1‑146 (MQ=255)                                             tCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTcagcag                                                                                              <  2:173523/146‑1 (MQ=255)                                               ggCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGa                                                                                          <  6:276512/148‑1 (MQ=255)                                                    cccacATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCAt                                                                                     >  8:161167/1‑148 (MQ=255)                                                      cacaTCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATaaa                                                                                  <  3:209498/149‑1 (MQ=255)                                                      cacaTCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATaaa                                                                                  <  2:101387/149‑1 (MQ=255)                                                           cACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTAtt                                                                                                                                     <  4:148161/93‑1 (MQ=255)                                                           cACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTAtt                                                                                                                                     >  3:148161/1‑93 (MQ=255)                                                             ccGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGActc                                                                           >  4:15947/1‑149 (MQ=255)                                                             ccGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGActc                                                                           <  8:223996/149‑1 (MQ=255)                                                             ccGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGAct                                                                            >  7:223996/1‑148 (MQ=255)                                                                cccACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCttg                                                                        <  7:90070/149‑1 (MQ=255)                                                                  cACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCttgtt                                                                      <  5:256247/149‑1 (MQ=255)                                                                     gacgCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTtgtg                                                                   >  2:12459/1‑149 (MQ=255)                                                                      acgCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTgg                                                                  <  7:79560/149‑1 (MQ=255)                                                                          gTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTgg                                                                  >  3:225666/1‑145 (MQ=255)                                                                           tAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGAt                                                             <  3:15947/149‑1 (MQ=255)                                                                                   gCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGtt                                                     <  1:8467/149‑1 (MQ=255)                                                                                         gCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAAt                                               <  1:12459/149‑1 (MQ=255)                                                                                          ccTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAAt                                               >  1:189466/1‑148 (MQ=255)                                                                                              gCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGAt                                            >  4:268567/1‑147 (MQ=255)                                                                                              gCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGAt                                            <  6:176124/147‑1 (MQ=255)                                                                                              gCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATaa                                          <  4:225666/149‑1 (MQ=255)                                                                                                 acacCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCa                                       <  6:2894/149‑1 (MQ=255)                                                                                                  cacCAGCTTGAGCAGGTCGAAGAGCAGCAGCTGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGACGGTCAGCAGTCGACGCATAAAGATACTCTTGTTGTgg                                                                  <  5:152971/121‑1 (MQ=255)                                                                                                   acCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGActc                                                                           >  6:178709/1‑111 (MQ=255)                                                                                                   acCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGActc                                                                           <  5:178709/111‑1 (MQ=255)                                                                                                       gCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTgg                                                                  <  6:49648/116‑1 (MQ=255)                                                                                                       gCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTgg                                                                  >  5:49648/1‑116 (MQ=255)                                                                                                         ttGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTg                               >  1:44954/1‑149 (MQ=255)                                                                                                         ttGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTg                               >  1:49537/1‑149 (MQ=255)                                                                                                                gggCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAg                                         >  7:73022/1‑132 (MQ=255)                                                                                                                gggCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAg                                         <  8:73022/132‑1 (MQ=255)                                                                                                                           agcagcGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACTCCAAATGTca              >  6:70057/1‑147 (MQ=255)                                                                                                                            gcagcGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCaa            >  7:132510/1‑149 (MQ=255)                                                                                                                                    ctgctCAATACCGTGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGAc    >  7:288432/5‑149 (MQ=255)                                                                                                                                      gctCAATACCGTGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGAcgg  >  7:4753/3‑149 (MQ=255)                                                                                                                                               cACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGtt                                                     <  6:46731/89‑1 (MQ=255)                                                                                                                                               cACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGtt                                                     >  5:46731/1‑89 (MQ=255)                                                                                                                                                |                                                                                                                                           GGATGCACAGGGCGAAGATGGCTTTGTCGATTTCCGACCAGGTTTCGGCCACCCACATCACCGCCCACGACGCGTAGTCGCGGCCACGGCCTGGCAACACCAGCTTGAGCAGGGCGAAGACCAGCAGCGGGCAGATCAATACCACGGTATTGAGCAGCAGCAGGGTGGTAGTGAGAATGCCGGTCAGCAGGCGACGCATAAAGAGACTCTTGTTGTGGGAGATAAACGGTTGGCAATGATAAGCAGGGCGGTGGCCTACGCCAAATGTCCAATGCCCGACGG  >  NC_002947/399385‑399666

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