Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	6,180,730	C→G	W258S (TGG→TCG)	atpB ←	ATP synthase subunit a

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,180,730	0	C	G	96.7%	190.7 / ‑2.3	62	W258S (TGG→TCG)	atpB	ATP synthase subunit a
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base G (27/32);  minor base A (1/1);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

TCCCCCCGTACAGGCAGGGAAGGCTATCAGACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCCACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGTTGCCGAACAGACGCAGTGCCAGCGAGATCGGCTTGGCGATCAGG  >  NC_002947/6180588‑6180873                                                                                                                                               |                                                                                                                                                cccccccGTACAGGCAGGGAAGGCTATCAGACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCATGATGTGGAACACAGCCGACGCCCa                                                                                                                                          <  2:1330489/149‑1 (MQ=255)             ggcaggGAAGGCTATCAGACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGcaccac                                                                                                                              >  1:664611/1‑150 (MQ=255)              gcaggGAAGGCTATCAGACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACg                                                                                                                                              >  2:864554/1‑133 (MQ=255)              gcaggGAAGGCTATCAGACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACg                                                                                                                                              <  1:864554/133‑1 (MQ=255)                aggGAAGGCTATCAGACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACAcc                                                                                                                           <  2:664611/150‑1 (MQ=255)                  ggAAGGCTATCAGACACGCCAGGCGGTCTGAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCa                                                                                                                         <  1:1578007/150‑1 (MQ=255)                      ggCTATCAGACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCAACGCCCATTGCAGCACCACACCCAGGcc                                                                                                                     >  1:1676236/1‑150 (MQ=255)                       gCTATCAGACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCg                                                                                                                    >  2:506069/1‑150 (MQ=255)                         tATCAGACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCt                                                                                                                  >  2:436856/1‑150 (MQ=255)                         tATCAGACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCt                                                                                                                  <  1:1345792/150‑1 (MQ=255)                           tCAGACACGCCAGGCGGTCTTAATGTTTATCTGCGAAAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGa                                                                                                                <  1:964686/150‑1 (MQ=255)                           tCAGACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGa                                                                                                                <  1:1616286/150‑1 (MQ=255)                              gACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGcc                                                                                                             >  1:1363958/1‑150 (MQ=255)                               acacGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGa                                                                                                                <  1:620622/146‑1 (MQ=255)                               acacGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGa                                                                                                                >  2:620622/1‑146 (MQ=255)                                 acGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCcaca                                                                                                          <  2:963350/150‑1 (MQ=255)                                   gCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCAcagc                                                                                                        >  2:1510554/1‑150 (MQ=255)                                      aGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAgg                                                                                                     <  1:119402/150‑1 (MQ=255)                                          ggTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCgg                                                                                                 >  2:371479/1‑150 (MQ=255)                                           gTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGa                                                                                                <  2:668046/150‑1 (MQ=255)                                             cTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCAACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGc                                                                                              <  2:1676236/150‑1 (MQ=255)                                              ttAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGcc                                                                                             >  1:887355/1‑150 (MQ=255)                                              ttAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGcc                                                                                             >  2:77145/1‑150 (MQ=255)                                                 aTGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCGTTACATCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGaa                                                                                          >  1:1498476/1‑150 (MQ=255)                                                  tGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGAc                                                                                                                                               <  2:738524/96‑1 (MQ=255)                                                    gTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACAt                                                                                       >  2:1128020/1‑150 (MQ=255)                                                      tATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTa                                                                                     <  1:1510554/150‑1 (MQ=255)                                                        tctcCGGGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGACCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACa                                                                                   <  2:126960/150‑1 (MQ=255)                                                               gAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCa                                                                            >  1:175178/1‑150 (MQ=255)                                                                  ccATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGa                                                                         >  1:1305874/1‑150 (MQ=255)                                                                   cATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGAt                                                                        <  2:1135061/150‑1 (MQ=255)                                                                          aGGTAGACGATGATCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCCACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACAcc                                                                 <  2:995118/150‑1 (MQ=255)                                                                           ggTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCa                                                                <  1:436856/150‑1 (MQ=255)                                                                            gTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGa                                                                                                >  2:268524/1‑117 (MQ=255)                                                                            gTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGa                                                                                                <  1:268524/117‑1 (MQ=255)                                                                            gTAGACGATGCTCAGCATCATGACGCTCTTACCTTGCAGGGTGATGACCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGNCACAGCNGGCCGGAGCCGANCATTACAGCGATCAGGATGAACACCAg                                                               <  1:1311641/150‑1 (MQ=255)                                                                             tAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGc                                                              >  1:1715041/1‑150 (MQ=255)                                                                               gACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGc                                                                                                    <  2:281200/110‑1 (MQ=255)                                                                               gACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGc                                                                                                    >  1:281200/1‑110 (MQ=255)                                                                                  gATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGcc                                                         <  2:1363958/150‑1 (MQ=255)                                                                                       tCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCAt                                                    >  2:779370/1‑150 (MQ=255)                                                                                              catGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGtt                                             <  1:77145/150‑1 (MQ=255)                                                                                               atGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGTTg                                            >  1:1304645/1‑150 (MQ=255)                                                                                                tGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGTTGc                                           <  2:887355/150‑1 (MQ=255)                                                                                                      tGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGTTGCCGAACa                                     <  1:371479/150‑1 (MQ=255)                                                                                                         aaGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGAt                                                                        >  2:19123/1‑112 (MQ=255)                                                                                                         aaGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGAt                                                                        <  1:19123/112‑1 (MQ=255)                                                                                                             ttGCAGGGTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGTTGCCGAACAGACGCAg                              >  2:1032383/1‑150 (MQ=255)                                                                                                                   ggTGATGATCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGTTGCCGAACAGACGCAGTGCCAg                        >  2:301512/1‑150 (MQ=255)                                                                                                                        tgatCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCa                                                     <  1:472655/116‑1 (MQ=255)                                                                                                                        tgatCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCa                                                     >  2:472655/1‑116 (MQ=255)                                                                                                                          atCAGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGTTGCCGAACAGACGCAGTGCCAGCGAGATc                 >  1:1307538/1‑150 (MQ=255)                                                                                                                             aGGATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGTCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGGTCAGGATGAACACCAGCTCGCCAGCATACATGTTGCCGAACAGACGCAGTGCCAGCGAGATCGGc              <  2:1257581/150‑1 (MQ=255)                                                                                                                               gATGTGGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGTTGCCGAACAGATGCAGTGCCAGCGAGATCGGCtt            <  1:779370/150‑1 (MQ=255)                                                                                                                                 tgtgGAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGTTGCCGAACAGACGCAGTGCCAGCGAGATCGGCTTgg          <  1:1128020/150‑1 (MQ=255)                                                                                                                                     gAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAg                                                                                  >  1:1382659/1‑74 (MQ=255)                                                                                                                                     gAACACAGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAg                                                                                  <  2:1382659/74‑1 (MQ=255)                                                                                                                                     gAACACAGCCGACGCCCATTGCAGCACCACACCCACGGCGCTGAGCCACATCAGGCCGGAGCCGAACATGACAGCCACGAAAAAGAACGGGAGCGCGCCAGCATACATGTTGCCGAACAGACGCAGTGCCAGCGAGATCGGCTTGGCGAt      >  1:834127/1‑150 (MQ=255)                                                                                                                                        cacaGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAg                                                                                               >  1:838104/1‑58 (MQ=255)                                                                                                                                        cacaGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAg                                                                                               <  2:838104/58‑1 (MQ=255)                                                                                                                                         acaGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGTTGCCGAACAGACGCAGTGCCAGCGAGATCGGCTTGGCGATCAgg  <  2:1304645/150‑1 (MQ=255)                                                                                                                                         acaGCCGACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGTTGCCGAACAGACGCAGTGCCAGCGAGATCGGCTTGGCGATCAgg  <  1:1235264/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                TCCCCCCGTACAGGCAGGGAAGGCTATCAGACACGCCAGGCGGTCTTAATGGTTATCTTCGTGAGCCATCGACAGGTAGACGATGGTCAGCATCATGAAGATGAAAGCTTGCAGGGTGATGATCAGGATGTGGAACACAGCCCACGCCCATTGCAGCACCACACCCAGGCCGCTGAGCCACAGCAGGCCGGAGCCGAACATTACAGCGATCAGGATGAACACCAGCTCGCCAGCATACATGTTGCCGAACAGACGCAGTGCCAGCGAGATCGGCTTGGCGATCAGG  >  NC_002947/6180588‑6180873

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