Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	950,367	G→A	7.5%	G48D (GGC→GAC)	cyoA →	cytochrome bo terminal oxidase subunit II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	950,367	0	G	A	7.5%	151.9 / 2.3	53	G48D (GGC→GAC)	cyoA	cytochrome bo terminal oxidase subunit II
Reads supporting (aligned to +/- strand):  ref base G (23/26);  new base A (2/2);  total (25/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 = 7.00e-01

GAGTAAAAAGCGTTACCCCAGACTGTTTGGCATATTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCAAGATCGAGGCCGCGGTATGGATCATC  >  NC_002947/950227‑950506                                                                                                                                              |                                                                                                                                            gAGTAAAAAGCGTTACCCCAGACTGTTTGGCATATTGCCCTTTTTAGGCATGCTATTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCAAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGACTGAtgc                                                                                                                                     >  4:84753/1‑149 (MQ=255)  aGTAAAAAGCGTTACCCCAGACTGTTTGGCATATTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGAtgc                                                                                                                                     <  8:140875/148‑1 (MQ=255)  aGTAAAAAGCGTTACCCCAGACTGTTTGGCATATAGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGAtgct                                                                                                                                    >  1:49377/1‑149 (MQ=255)   ttAAAAACCGTTACCCCAGACTTTTTGGCATATTGCCCTTTTTAGGGATTCTTTTACTCATTGGGTGGAACTGGACCCTTCTCGACCCGAAGGGCCATGTCGGCATTGAGCAAAAGAACTCTGATCCTGATCGCTACCGGCCTGAtgc                                                                                                                                     <  1:22383/147‑1 (MQ=255)           cGTTACCCCAGACTGTTTGGCATATTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTggtggt                                                                                                                           <  3:268089/149‑1 (MQ=255)                cccgAGACTGTTTTGCATATTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATcc                                                                                                                      <  2:325536/149‑1 (MQ=255)                 cccAGACTGTTTGGCATATTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGc                                                                                                                                            >  8:238740/1‑126 (MQ=255)                 cccAGACTGTTTGGCATATTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGc                                                                                                                                            <  7:238740/126‑1 (MQ=255)                    aGACTGTTTGGCATATTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCtgt                                                                                                                  <  8:131898/149‑1 (MQ=255)                    aGACTGTTTGGCATATTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCtgt                                                                                                                  >  6:137197/1‑149 (MQ=255)                      aCTGTTTGGCATATTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑ATGATCCTGATCGCTACCGGCCTGATGCTGCTGTTGGTAATCCCTGTGa                                                                                                                >  5:266182/1‑149 (MQ=255)                            tGGCATATTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGatcatc                                                                                                           <  1:327890/148‑1 (MQ=255)                                atatTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGACCTGATGCTGCTGGTGGTAATcc                                                                                                                      <  6:10541/133‑1 (MQ=255)                                atatTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGACCTGATGCTGCTGGTGGTAATcc                                                                                                                      >  5:10541/1‑133 (MQ=255)                                  atTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGAccc                                                                                                    >  1:19232/1‑149 (MQ=255)                                     gCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTgg                                                                                                 >  6:178499/1‑149 (MQ=255)                                      cccTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGc                                                                                                <  5:178499/149‑1 (MQ=255)                                       ccTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCg                                                                                               >  5:340774/1‑149 (MQ=255)                                                    gCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGc                                                                                                >  8:246510/1‑135 (MQ=255)                                                    gCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGc                                                                                                <  7:246510/135‑1 (MQ=255)                                                      ttttACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTAc                                                                                >  5:170458/1‑149 (MQ=255)                                                              aGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTc                                                                        >  7:90398/1‑149 (MQ=255)                                                               gTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTc                                                                        >  3:204301/1‑148 (MQ=255)                                                                  ggTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCcaa                                                                     <  2:19232/148‑1 (MQ=255)                                                                            gACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCa                                                          >  1:52602/1‑149 (MQ=255)                                                                            gACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCa                                                          >  7:312443/1‑149 (MQ=255)                                                                            gACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCa                                                          >  6:144335/1‑149 (MQ=255)                                                                            gACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCGGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCa                                                          >  5:270782/1‑149 (MQ=255)                                                                              cccTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCAc                                                         >  8:119109/1‑148 (MQ=255)                                                                              cccTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCAc                                                         <  7:119109/148‑1 (MQ=255)                                                                               ccTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGCTGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCATCACGGCAGCCACCt                                                       >  8:90848/1‑149 (MQ=255)                                                                               ccTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCAcc                                                        >  4:180433/1‑148 (MQ=255)                                                                               ccTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCAcc                                                        >  7:331846/1‑148 (MQ=255)                                                                                   cTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTacac                                                   <  4:204301/149‑1 (MQ=255)                                                                                       aCCCGAATGGCCAGGGCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCt                                               <  6:266182/149‑1 (MQ=255)                                                                                       aCCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCt                                               <  2:49377/149‑1 (MQ=255)                                                                                       aCCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCt                                               <  4:172834/149‑1 (MQ=255)                                                                                              gggCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGt                                        <  6:340774/149‑1 (MQ=255)                                                                                               ggCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTc                                       <  6:170458/149‑1 (MQ=255)                                                                                                    ggTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCcaa                                                                     <  3:147905/114‑1 (MQ=255)                                                                                                    ggTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCcaa                                                                     >  4:147905/1‑114 (MQ=255)                                                                                                         gCATTGATCAAAAGAAC‑CTGATCGTGATCGCTACGGGCCTGATGCTTCTTGTGGTAATCCCTGTGATCATCATTACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACCCCCCTGACTGGTCGCACTCGAcc                             <  8:90398/149‑1 (MQ=255)                                                                                                           aTTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCaa                           <  2:52602/149‑1 (MQ=255)                                                                                                                 cAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCAAGATCGa                     <  3:180433/149‑1 (MQ=255)                                                                                                                 cAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCAAGATCGa                     <  8:312443/149‑1 (MQ=255)                                                                                                                         c‑ctgatcctgatCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCAAGATCGAGGCCGCgg             <  5:144335/149‑1 (MQ=255)                                                                                                                            tgatcctgatcGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTg                                                                                                  >  2:183541/1‑62 (MQ=255)                                                                                                                            tgatcctgatcGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTg                                                                                                  <  1:183541/62‑1 (MQ=255)                                                                                                                                 ctgatcGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCAAGATCGAGGCCGCGGTATGGa       <  5:137197/148‑1 (MQ=255)                                                                                                                                 ctgatcGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCAAGATCGAGGCCGCGGTATGGa       >  3:270368/1‑148 (MQ=255)                                                                                                                                   gatcGCTACCGACCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCAAGATc                       >  5:151427/1‑130 (MQ=255)                                                                                                                                   gatcGCTACCGACCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCAAGATc                       <  6:151427/130‑1 (MQ=255)                                                                                                                                     tcGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCAAGATCGAGGCCGCGGTATGGatcatc  <  8:331846/149‑1 (MQ=255)                                                                                                                                              gCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCAAGATCGAGGCCGCGGTATGGat      <  7:248335/136‑1 (MQ=255)                                                                                                                                              gCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCAAGATCGAGGCCGCGGTATGGat      >  8:248335/1‑136 (MQ=255)                                                                                                                                              |                                                                                                                                            GAGTAAAAAGCGTTACCCCAGACTGTTTGGCATATTGCCCTTTTTAGGCATGCTTTTACTCAGTGGGTGCAACTGGACCCTGCTCGACCCGAAGGGCCAGGTCGGCATTGAGCAAAAGAAC‑CTGATCCTGATCGCTACCGGCCTGATGCTGCTGGTGGTAATCCCTGTGATCATCATGACCCTGGCGTTCGCCTGGAAGTACCGTGCTTCCAACAAGGCAGCCACCTACACCCCTGACTGGTCGCACTCGACCAAGATCGAGGCCGCGGTATGGATCATC  >  NC_002947/950227‑950506

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