Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,385,150	A→C	13.7%	intergenic (‑120/+41)	proS ← / ← oprD	proline‑‑tRNA ligase/basic amino acid specific porin OprD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,385,150	0	A	C	13.7%	129.8 / 15.2	51	intergenic (‑120/+41)	proS/oprD	proline‑‑tRNA ligase/basic amino acid specific porin OprD
Reads supporting (aligned to +/- strand):  ref base A (21/23);  new base C (4/3);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

AAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGA  >  NC_002947/1385011‑1385287                                                                                                                                             |                                                                                                                                          aaaTATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCggg                                                                                                                                   <  7:81232/148‑1 (MQ=255)      ttGACTGGTGCGCAGGGTTAGCCGTGTCGATTGCCTCAGACGTAGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACGCCGGGCttttt                                                                                                                             <  4:43662/149‑1 (MQ=255)                cgcATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATCGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCgg                                                                                                                    <  2:41300/148‑1 (MQ=255)                   aTGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGc                                                                                                                                  <  5:215072/131‑1 (MQ=255)                   aTGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGc                                                                                                                                  >  6:215072/1‑131 (MQ=255)                      gTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGTGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACGCCGGGCTTTTTc                                                                                                                            <  5:106642/134‑1 (MQ=38)                      gTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACGCCGGGCTTTTTc                                                                                                                            >  6:106642/1‑134 (MQ=39)                         aGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGAAGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGc                                                                                                            >  1:88991/1‑148 (MQ=255)                            ctagtagATTTCTTAAGACGTTTAAATATCCCTGAATTTTACGGGGT‑GTAAACGAAGGCGTACAGGATGCCCCTTGGCGTGATACTTGTCAGCCAAGAAAAAGCTCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGcc                                                                                                       <  4:43055/143‑1 (MQ=255)                             gtgtCGATTGCCTAAGACGTTGAATTAGGCGCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGGCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCa                                                                                                      <  2:285317/149‑1 (MQ=255)                                   aTTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTa                                                                                                >  3:32694/1‑149 (MQ=255)                                      gCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACa                                                                                              <  4:222778/148‑1 (MQ=255)                                        cTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGa                                                                                           >  6:187633/1‑149 (MQ=255)                                               gTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAg                                                                                    >  5:236500/1‑149 (MQ=255)                                                ttGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGc                                                                                   <  6:236500/149‑1 (MQ=255)                                                   aaaTAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCgg                                                                                 <  5:199762/148‑1 (MQ=255)                                                           ccGCATTGTACGGGGT‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCtttt                                                                                                                              <  6:333440/95‑1 (MQ=255)                                                           ccGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTTACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCAc                                                                        <  6:50976/149‑1 (MQ=255)                                                           ccGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCtttt                                                                                                                              >  5:333440/1‑95 (MQ=255)                                                               aTTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATc                                                                    <  6:111946/149‑1 (MQ=255)                                                                            ggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGc                                                                                   >  2:354454/1‑122 (MQ=255)                                                                            ggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGc                                                                                   <  1:354454/122‑1 (MQ=255)                                                                            ggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTtcg                                                        >  3:350306/1‑149 (MQ=255)                                                                                 cGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcgat                                                   <  4:350306/149‑1 (MQ=255)                                                                                   aaGGCGTACAGGATGCACCAGCACGTGACACTTGTCAGCCAAGAAAAAGCCCAGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCAATGCAAGCCAGGCTTACAGGAAGCTCAGCGGGAACTCCACCATCAGGCGCAATTCGtagatc                                                  >  2:82609/1‑148 (MQ=255)                                                                                     ggCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCGGGCTTACAGGATGCTCAGCGGGTACTCCAc                                                                        <  6:48713/124‑1 (MQ=255)                                                                                     ggCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCGGGCTTACAGGATGCTCAGCGGGTACTCCAc                                                                        >  5:48713/1‑124 (MQ=255)                                                                                         tACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATggg                                            >  4:302010/1‑148 (MQ=255)                                                                                              gATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCg                                      >  6:353806/1‑149 (MQ=255)                                                                                                  ccccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtagt                                  <  1:82609/149‑1 (MQ=255)                                                                                                   cccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGCCGATCGATGGGGAATCgtagta                                 <  5:187633/149‑1 (MQ=255)                                                                                                    ccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtag                                   >  5:324288/1‑146 (MQ=255)                                                                                                     cAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGa                               <  5:353806/149‑1 (MQ=255)                                                                                                           gTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATc                                                                    >  7:194481/1‑106 (MQ=255)                                                                                                           gTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATc                                                                    <  8:194481/106‑1 (MQ=255)                                                                                                                acTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCTGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg                      >  1:253051/1‑147 (MQ=255)                                                                                                                acTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCTGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg                      >  1:118786/1‑147 (MQ=255)                                                                                                                acTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg                      >  5:33770/1‑147 (MQ=255)                                                                                                                 cTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg                      >  5:276163/1‑146 (MQ=255)                                                                                                                 cTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCgg                      >  7:245507/1‑146 (MQ=255)                                                                                                                  ttGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGgaag                   >  5:157748/1‑148 (MQ=255)                                                                                                                           aaGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCgg                                                                                 >  2:189225/1‑77 (MQ=255)                                                                                                                           aaGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCgg                                                                                 <  1:189225/77‑1 (MQ=255)                                                                                                                             gAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTgg       <  6:276163/149‑1 (MQ=255)                                                                                                                             gAAAAAGCCCGGCAATCACGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAgg                                                                 <  2:299453/91‑1 (MQ=37)                                                                                                                             gAAAAAGCCCGGCAATCACGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAgg                                                                 >  1:299453/1‑91 (MQ=37)                                                                                                                             gAAAAAGCCCGGCAATCACGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTg        >  7:136036/1‑148 (MQ=255)                                                                                                                             gAAAAAGCCCGGCAATCACGCCGGGCGTTTTCAGGTGCGGGGCCCTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGGTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTgg       >  7:140750/1‑149 (MQ=255)                                                                                                                              taacaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGgcg                                                               <  6:338109/88‑1 (MQ=255)                                                                                                                                aaaGCCCGGCGTGATTGCCGGGCTTTTTCTGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACACCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGCGGCCt    <  2:118786/149‑1 (MQ=255)                                                                                                                                  aGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGa  <  2:88991/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                          AAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGG‑GGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGA  >  NC_002947/1385011‑1385287

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