Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,624,668	G→A	7.7%	Q365* (CAA→TAA)	PP_4943 ←	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,624,668	0	G	A	7.7%	145.6 / 2.6	53	Q365* (CAA→TAA)	PP_4943	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (24/24);  new base A (2/2);  total (26/27)
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.96e-01

CAGCGGATGGGTGGCAACACTGCCGCCAAAGCGCGCCCACGCCTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTGGGCCAGCTGCGCGGCGTCGCCCAGCGGGAA  >  NC_002947/5624520‑5624815                                                                                                                                                     |                                                                                                                                                    caGCGGATGGGTGGCAACACTGCCGCCAAAGCGCGCCCACGCCTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCGCCACCAACGCGCGCACCTGTGGGAGTTg                                                                                                                                                     >  6:308023/1‑149 (MQ=255)   gCGGATGGGTGGCAACACTGCCGCCAAAGCGCGCCCACGCCTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGc                                                                                                                                                    >  6:269131/1‑148 (MQ=255)              gCAACACTGCCGCCAAAGCGCGCCCACGCCTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTc                                                                                                                                         <  2:105738/148‑1 (MQ=255)                  cacTGCCGCCAAAGCGCGCCCACGCCTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCa                                                                                                                                    <  5:269131/149‑1 (MQ=255)                           cAAAGCGCGCCCACGCCTGGTCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGGCTGGATACGATTGAGCATCAGGCCTGCCCCACCATCGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGAt                                                                                                                           <  3:75247/149‑1 (MQ=255)                           cAAAGCGCGCCCACGCCTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGAt                                                                                                                           <  1:311431/149‑1 (MQ=255)                           cAAAGCGCGCCCACGCCTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGAt                                                                                                                           >  7:279648/1‑149 (MQ=255)                                   gcCCACGCCTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTCTGCGAGTTGCCAGAAGGCCTCGCCCACCGCCTCCTCAGAGaac                                                                                                                   >  5:254756/1‑148 (MQ=255)                                     ccACGCCTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGTGCACCGCCTGATCAGAGAAACg                                                                                                                 >  3:81003/1‑149 (MQ=255)                                      cACGCCTGGGCATAAGCCTTGGCATCGATAGCTTTCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCTTCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGc                                                                                                                <  7:194415/149‑1 (MQ=255)                                        cGCCTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTcgc                                                                                                                                       >  1:29140/1‑124 (MQ=255)                                        cGCCTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTcgc                                                                                                                                       <  2:29140/124‑1 (MQ=255)                                          ccTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCAc                                                                                                            >  1:266414/1‑149 (MQ=255)                                            tGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTACCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCAcc                                                                                                           >  6:325879/1‑148 (MQ=255)                                                cATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGgc                                                                                                      <  5:29075/149‑1 (MQ=255)                                                cATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGgc                                                                                                      >  4:197288/1‑149 (MQ=255)                                                cATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGgc                                                                                                      >  8:28347/1‑149 (MQ=255)                                                cATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGgc                                                                                                      <  2:178362/149‑1 (MQ=255)                                                   aaGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGgcgcg                                                                                                   >  2:207106/1‑148 (MQ=255)                                                    aGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGCCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTACCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCt                                                                                                   <  5:325879/148‑1 (MQ=255)                                                     gCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTg                                                                                                  >  6:30986/1‑148 (MQ=255)                                                           gCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTTGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCAtg                                                                                           >  1:219408/1‑149 (MQ=255)                                                             aTCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCAtgtg                                                                                         >  5:115483/1‑149 (MQ=255)                                                              tCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGaa                                                                                                                                               <  8:169406/94‑1 (MQ=255)                                                              tCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGaa                                                                                                                                               >  7:169406/1‑94 (MQ=255)                                                              tCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGc                                                                                        <  6:203421/149‑1 (MQ=255)                                                                gATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCt                                                                                      <  5:30986/149‑1 (MQ=255)                                                                    gCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGca                                                                                  <  2:375532/149‑1 (MQ=255)                                                                      ttGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGca                                                                                  >  6:267587/1‑147 (MQ=255)                                                                        gCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCACCCCTGGGGGAGTTGCCAGAAGGCCTCGCGCCCCGCCGGGTCCGAGAAACGCGCCCCCAGGCGCTGCCGCATGTGCCTGGccccgg                                                                              >  1:111012/1‑144 (MQ=255)                                                                              cgcgTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGgc                                                                        <  2:266414/149‑1 (MQ=255)                                                                                      cgGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCa                                                                <  3:197288/149‑1 (MQ=255)                                                                                           ggCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTTGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTc                                                           <  2:219408/149‑1 (MQ=255)                                                                                            gCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTc                                                           >  2:18261/1‑148 (MQ=255)                                                                                             ccTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCa                                                         <  1:203408/149‑1 (MQ=255)                                                                                                 gATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTAc                                                     >  7:355906/1‑149 (MQ=255)                                                                                                   tACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACgg                                                   <  8:279648/149‑1 (MQ=255)                                                                                                      gATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTgg                                                                         <  4:200495/124‑1 (MQ=255)                                                                                                      gATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTgg                                                                         >  3:200495/1‑124 (MQ=255)                                                                                                      gATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCa                                                >  5:131447/1‑149 (MQ=255)                                                                                                      gATTGAGCATCAGGCCTGCCACACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCa                                                >  4:243723/1‑149 (MQ=255)                                                                                                       aTTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCAt                                               <  5:139218/149‑1 (MQ=255)                                                                                                              aTCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGtgctgc                                        <  5:267587/149‑1 (MQ=255)                                                                                                                         cccACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTGGGc                             >  6:161118/1‑149 (MQ=255)                                                                                                                          ccaccaACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTGGGc                             >  1:148619/1‑148 (MQ=255)                                                                                                                           caccaACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTGGGc                             >  6:346527/1‑147 (MQ=255)                                                                                                                           caccaACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTGGGc                             >  3:171600/1‑147 (MQ=255)                                                                                                                            accaACGCGCGCTCCTGTGGGAGTTGCCAGCAGGCCTCGCGGACCGCGTGATCAGAGAACCGCTCCACCAGGCGCTGCAGGATGTGCCTGGCACAGGCCTTGCGCTGCCCAGCGTCCAGTACGGCCTTGTGCTGCAGGGCCTGGGCCAg                          <  6:254756/149‑1 (MQ=255)                                                                                                                                  gcgcgcACCTGTGGGAGTTTCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTGGGCCAGCTgcgc                    <  1:18261/149‑1 (MQ=255)                                                                                                                                        aCCTGTGGGAGTTACCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTgg                               >  6:241519/1‑132 (MQ=255)                                                                                                                                        aCCTGTGGGAGTTACCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTgg                               <  5:241519/132‑1 (MQ=255)                                                                                                                                                gAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTGGGCCAGCTGCGCGGCGTCGCCCAGCg      <  6:110119/149‑1 (MQ=255)                                                                                                                                                  gTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTGGGCCAGCTGCGCGGCGTCGCCCAGCggg    <  5:161118/149‑1 (MQ=255)                                                                                                                                                   ttGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTGGGCCAGCTGCGCGGCGTCGCCCAGCGGGa   <  7:28347/149‑1 (MQ=255)                                                                                                                                                     gCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTGGGCCAGCTGCGCGGCGTCGCCCAGCGGGaa  <  1:207106/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CAGCGGATGGGTGGCAACACTGCCGCCAAAGCGCGCCCACGCCTGGGCATAAGCCTTGGCATCGATAGCTTGCCAACCGCGTTCGCGGAAGGCCTGGATACGATTGAGCATCAGGCCTGCCCCACCAACGCGCGCACCTGTGGGAGTTGCCAGAAGGCCTCGCGCACCGCCTGATCAGAGAAACGCTCACCCAGGCGCTGCAGCATGTGCCTGGCACAGGCCTGGCGCTGCCCAGCGTCCAGTACGGCCATGTGCTGCAGGCCCTGGGCCAGCTGCGCGGCGTCGCCCAGCGGGAA  >  NC_002947/5624520‑5624815

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