Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,917,965	G→A	8.7%	V122I (GTC→ATC) ‡	PP_4328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,917,965	0	G	A	8.7%	160.8 / 4.1	57	V122I (GTC→ATC) ‡	PP_4328	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/23);  new base A (2/3);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCC  >  NC_002947/4917822‑4918106                                                                                                                                                |                                                                                                                                              aGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGc                                                                                                                                          <  8:226890/149‑1 (MQ=255)     cAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCCCCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCt                                                                                                                                      >  4:83659/1‑149 (MQ=255)          ccagccagCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGCCCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTtgctg                                                                                                                                 >  7:7849/1‑149 (MQ=255)          ccagccagCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgctg                                                                                                                                 <  6:223832/149‑1 (MQ=255)               atctggCTGCCCTGGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgcaagg                                                                                                                            >  3:291043/6‑149 (MQ=255)                agccGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgcaagg                                                                                                                            >  2:247759/1‑148 (MQ=255)                      tGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTc                                                                                                                     >  4:34155/1‑149 (MQ=255)                        ccccAGGGCAGCCAGATCACCGTGAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCt                                                                                                                   <  6:228566/149‑1 (MQ=255)                        ccccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCt                                                                                                                   <  3:83659/149‑1 (MQ=255)                         cccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCt                                                                                                                   >  6:24044/1‑148 (MQ=255)                         cccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCt                                                                                                                   <  5:24044/148‑1 (MQ=255)                          ccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGa                                                                                                                 >  5:242103/1‑149 (MQ=255)                                  gCCAGATCAACGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAg                                                                                                         <  6:255963/149‑1 (MQ=255)                                     aGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCt                                                                                                                   <  1:112499/136‑1 (MQ=255)                                     aGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCt                                                                                                                   >  2:112499/1‑136 (MQ=255)                                           ccGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgcaaggcaag                                                                                                                        <  3:305170/125‑1 (MQ=255)                                           ccGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgcaaggcaag                                                                                                                        >  4:305170/1‑125 (MQ=255)                                           ccGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTcc                                                                                                                    <  4:88721/129‑1 (MQ=255)                                           ccGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTcc                                                                                                                    >  3:88721/1‑129 (MQ=255)                                              gcAGCCAGACGGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAg                                                                                             <  8:263650/148‑1 (MQ=255)                                                  ccAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCAc                                                                                         >  8:301370/1‑149 (MQ=255)                                                     gACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCg                                                                                      <  1:249615/149‑1 (MQ=255)                                                      aCCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGccc                                                                                               >  4:148714/1‑139 (MQ=255)                                                      aCCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGccc                                                                                               <  3:148714/139‑1 (MQ=255)                                                       ccGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCa                                                                                              <  6:275489/139‑1 (MQ=255)                                                       ccGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCa                                                                                              >  5:275489/1‑139 (MQ=255)                                                        cGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGc                                                                                    >  3:1609/1‑148 (MQ=255)                                                        cGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgc                                                                                                                                   <  5:199560/101‑1 (MQ=255)                                                        cGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgc                                                                                                                                   >  6:199560/1‑101 (MQ=255)                                                           aaGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAAc                                                                                >  2:258609/1‑149 (MQ=255)                                                            aGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAAcc                                                                               >  5:12239/1‑149 (MQ=255)                                                               aaTCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACAgga                                                                            >  5:276603/1‑148 (MQ=255)                                                               aaTCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGCTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGc                                                                            >  4:166235/1‑149 (MQ=255)                                                                  cGCCTGGCGATCATGCTGCAACACGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCgt                                                                         <  2:320414/149‑1 (MQ=255)                                                                                 ctgcAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGt                                                          >  6:68064/1‑149 (MQ=255)                                                                                     aaCAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGc                                                      <  4:1609/149‑1 (MQ=255)                                                                                     aaCAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGc                                                      <  4:53413/149‑1 (MQ=255)                                                                                                gccagccaGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAg                                           >  1:232463/1‑149 (MQ=255)                                                                                                 ccagccagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTc                                                >  8:268124/1‑143 (MQ=255)                                                                                                 ccagccagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTc                                                <  7:268124/143‑1 (MQ=255)                                                                                                          tCGCCACCCTCACCCAGCTGGACACCAGCAAGTTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGc                                  <  1:223876/148‑1 (MQ=255)                                                                                                             ccACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGt                              <  3:34155/149‑1 (MQ=255)                                                                                                                  cTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTg                             >  3:314820/1‑145 (MQ=255)                                                                                                                  cTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTg                             <  4:314820/145‑1 (MQ=255)                                                                                                                     aCCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCaa                                              >  4:84915/1‑125 (MQ=255)                                                                                                                     aCCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCaa                                              <  3:84915/125‑1 (MQ=255)                                                                                                                        cAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAg                   <  7:301370/149‑1 (MQ=255)                                                                                                                                ttgcAGCAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCaa                                              >  5:278016/4‑114 (MQ=255)                                                                                                                                ttgcAGCAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCaa                                              <  6:278016/111‑1 (MQ=255)                                                                                                                                cacCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGAc            <  6:12239/148‑1 (MQ=255)                                                                                                                                 acCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAg          >  6:42457/1‑149 (MQ=255)                                                                                                                                 acCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCGGCTCAACAGCAACCCGGCCGGTGCCACCCAGAGCCTCGCCAg          >  6:47367/1‑149 (MQ=255)                                                                                                                                    aGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAAc                                       >  1:63940/1‑117 (MQ=255)                                                                                                                                    aGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAAc                                       <  2:63940/117‑1 (MQ=255)                                                                                                                                    aGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgccc       >  6:29144/1‑149 (MQ=255)                                                                                                                                    aGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgccc       >  6:11490/1‑149 (MQ=255)                                                                                                                                         ggTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcccgcgcc  >  7:102125/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                              AGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCGCGCC  >  NC_002947/4917822‑4918106

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