Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,631,141	C→G	7.5%	L202L (CTC→CTG)	lon‑II →	DNA‑binding ATP‑dependent protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,631,141	0	C	G	7.5%	163.2 / 2.9	53	L202L (CTC→CTG)	lon‑II	DNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base G (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CAGCATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCT  >  NC_002947/2630994‑2631284                                                                                                                                                    |                                                                                                                                                 cAGCATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGAcctcc                                                                                                                                                 >  8:169215/1‑149 (MQ=255)    cATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTg                                                                                                                                               >  6:193136/1‑148 (MQ=255)    cATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGc                                                                                                                                              <  4:201750/149‑1 (MQ=255)                 gggCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAgcg                                                                                                                                 >  6:170083/1‑149 (MQ=255)                 gggCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAgcg                                                                                                                                 >  5:108281/1‑149 (MQ=255)                          gTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCgg                                                                                                                        <  5:170083/149‑1 (MQ=255)                             gACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCg                                                                                                                     <  7:169215/149‑1 (MQ=255)                               cacCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCgg‑g                                                                                                                  >  5:332416/1‑149 (MQ=255)                                  cATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTc                                                                                                                >  1:175234/1‑148 (MQ=255)                                                        aaGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGcagc                                                                                           >  1:230096/1‑147 (MQ=255)                                                          gATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGa                                                                                                      <  8:21249/134‑1 (MQ=255)                                                          gATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGa                                                                                                      >  7:21249/1‑134 (MQ=255)                                                                   gAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGtac                                                                               >  4:250117/1‑148 (MQ=255)                                                                     aGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGtactac                                                                            <  2:230096/149‑1 (MQ=255)                                                                      gCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTAcc                                                                           <  6:108281/149‑1 (MQ=255)                                                                        aGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCt                                                                          >  5:280389/1‑148 (MQ=255)                                                                             aTTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTc                                                                                                                <  3:275925/105‑1 (MQ=255)                                                                             aTTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTc                                                                                                                >  4:275925/1‑105 (MQ=255)                                                                              ttCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAg                                                                   >  5:179909/1‑149 (MQ=255)                                                                                   gAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCa                                                                 <  6:280389/146‑1 (MQ=255)                                                                                     aaTCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTCCTGGATGCCGAAAACGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTCAATGAGCAGATg                                                             >  3:70044/1‑148 (MQ=255)                                                                                       tcgtcgACTTGACTACCCGTGTCGAGCCCGTGCGGGCGCGGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAg                                                          <  6:179909/149‑1 (MQ=255)                                                                                          tcgACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGcagcc                                                                                          <  4:76701/114‑1 (MQ=255)                                                                                          tcgACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGcagcc                                                                                          >  3:76701/1‑114 (MQ=255)                                                                                              cTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTc                                                   <  6:332416/149‑1 (MQ=255)                                                                                               ttGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAgc                                                                                                <  8:154108/103‑1 (MQ=255)                                                                                                       ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAAc                                          <  5:132122/149‑1 (MQ=255)                                                                                                       ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAAc                                          >  8:187066/1‑149 (MQ=255)                                                                                                            tCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGc                                                        >  7:39921/1‑130 (MQ=255)                                                                                                            tCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGc                                                        <  8:39921/130‑1 (MQ=255)                                                                                                            tCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                    >  4:207222/1‑134 (MQ=255)                                                                                                            tCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                    <  3:207222/134‑1 (MQ=255)                                                                                                              gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATg                                                                                                    >  7:310074/1‑84 (MQ=255)                                                                                                              gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATg                                                                                                    <  2:261083/84‑1 (MQ=255)                                                                                                              gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATg                                                                                                    <  8:310074/84‑1 (MQ=255)                                                                                                              gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATg                                                                                                    >  1:261083/1‑84 (MQ=255)                                                                                                              gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                    >  6:236991/1‑132 (MQ=255)                                                                                                              gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                    <  5:236991/132‑1 (MQ=255)                                                                                                              gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGa                                   >  2:274954/1‑149 (MQ=255)                                                                                                                gCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTACTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGAt                                  <  2:175234/148‑1 (MQ=255)                                                                                                                    gTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGa                             >  4:210307/1‑149 (MQ=255)                                                                                                                        tGGCGCTGGTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCTGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGATAGAACTTGGCGATGGTGACGaa                         <  4:70044/149‑1 (MQ=255)                                                                                                                          gcgcTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAg                        <  3:210307/148‑1 (MQ=255)                                                                                                                                  ggATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                    <  5:334556/112‑1 (MQ=255)                                                                                                                                  ggATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                    >  6:334556/1‑112 (MQ=255)                                                                                                                                      gCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGAGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCagaaagaa                                           >  1:73767/1‑118 (MQ=255)                                                                                                                                      gCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCagaaagaa                                           <  2:73767/117‑1 (MQ=255)                                                                                                                                        gcttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATg                                 <  1:52766/120‑1 (MQ=255)                                                                                                                                        gcttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATg                                 >  2:52766/6‑125 (MQ=255)                                                                                                                                        gcttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTc          >  4:83729/6‑148 (MQ=255)                                                                                                                                        cGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATg                                 >  8:21392/1‑125 (MQ=255)                                                                                                                                               gACCTGCAGGAGGTCTAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCt  <  3:83729/149‑1 (MQ=255)                                                                                                                                               gACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGc   <  1:274954/148‑1 (MQ=255)                                                                                                                                               gACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCt  >  1:142440/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CAGCATCGACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGG‑GTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCT  >  NC_002947/2630994‑2631284

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