Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,408,518	C→T	7.1%	A74T (GCC→ACC)	PP_4752 ←	Xaa‑Pro aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,408,518	0	C	T	7.1%	165.6 / 2.7	56	A74T (GCC→ACC)	PP_4752	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base C (28/24);  new base T (2/2);  total (30/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.06e-01

ACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAACGCCGTGTCGAC‑CGGCGCGATGGCATGTTCAGTGACC  >  NC_002947/5408370‑5408660                                                                                                                                                     |                                                                                                                                                 aCGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGc                                                                                                                                                  >  1:368571/1‑149 (MQ=255)       tCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑Caa                                                                                                                                              >  7:116366/1‑146 (MQ=255)               gtgCTGGCTGCTGGTGAACTCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCg                                                                                                                                    <  3:267313/148‑1 (MQ=255)                tgCTGGCTGCTGGTGAACTCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGc                                                                                                                                  >  3:323558/1‑149 (MQ=255)                  ctggctgCTGGTGAACTCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCg                                                                                                                                 >  2:329495/1‑148 (MQ=255)                       tgctgGTGAACTCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATcc                                                                                                                           >  2:11857/1‑149 (MQ=255)                         ctgGTGAACTCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCg                                                                                                                                    <  6:14562/138‑1 (MQ=255)                         ctgGTGAACTCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCg                                                                                                                                    >  5:14562/1‑138 (MQ=255)                               aaCTCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGgcgcg                                                                                                                   >  6:197970/1‑149 (MQ=255)                                 cTCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGgcgcg                                                                                                                   >  1:357526/1‑147 (MQ=255)                                  tCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGT‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAg                                                                                                                >  7:301394/1‑149 (MQ=255)                                    gAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGTC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAg                                                                                                                <  2:368571/147‑1 (MQ=255)                                        aGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTc                                                                                                          <  5:233888/149‑1 (MQ=255)                                          caccacCGGGCCGTCGGTGGCGACGAACACATAGCGGCTCGGCGAGTGCATCACCCACAAACCAAGCTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTaa                                                                                                                              <  8:317291/127‑1 (MQ=255)                                          caccacCGGGCCGTCGGTGGCGACGAACACATAGAGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTaa                                                                                                                              >  7:317291/1‑127 (MQ=255)                                             caccGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGT‑CAACAGAAt                                                                                                                                        >  7:334169/1‑114 (MQ=255)                                             caccGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGT‑CAACAGAAt                                                                                                                                        <  8:334169/114‑1 (MQ=255)                                              accGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAAt                                                                                                    <  1:262104/149‑1 (MQ=255)                                                    ccGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTc                                                                                              >  1:328070/1‑149 (MQ=255)                                                               gacgaACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAAt                                                                                   >  3:32166/1‑149 (MQ=255)                                                                acgaACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATg                                                                                  <  2:357526/149‑1 (MQ=255)                                                                      acaTAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCaa                                                                            <  1:329495/149‑1 (MQ=255)                                                                       caTAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAg                                                                           <  4:323558/149‑1 (MQ=255)                                                                        aTAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAgg                                                                          >  5:271251/1‑149 (MQ=255)                                                                                tGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGcga                                                                  >  2:326557/1‑149 (MQ=255)                                                                                       gTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGT‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATgg                                                           <  8:301394/149‑1 (MQ=255)                                                                                         gCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATggcg                                                         >  1:148398/1‑149 (MQ=255)                                                                                            tCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCa                                                      >  4:250266/1‑149 (MQ=255)                                                                                              aCCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTAGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTCCGATCGATGGCGGCATa                                                    >  3:47515/1‑149 (MQ=255)                                                                                               cccACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTTGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGTTAATGGCGCAAGGCTTTGc                                                                    <  3:102876/132‑1 (MQ=255)                                                                                               cccACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATAGGC‑CAACAGAATGCCGGCGTAATCCTGGGAGCGCAGTTGCTCGCGAATGCGCTCAAGGCGGTAATGGAGCAAGGCTTTGc                                                                    >  4:102876/1‑132 (MQ=255)                                                                                                  acaAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGcc                                                >  4:100720/1‑149 (MQ=255)                                                                                                  acaAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGcc                                                >  6:5819/1‑149 (MQ=255)                                                                                                    aaaCCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCt                                               <  2:181850/148‑1 (MQ=255)                                                                                                    aaaCCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCt                                               >  1:181850/1‑148 (MQ=255)                                                                                                      aCCAAGGTTATTGGTGTCGGTGTCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAg                                                                           <  2:273788/118‑1 (MQ=255)                                                                                                      aCCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGCACAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAg                                                                           >  1:273788/1‑119 (MQ=255)                                                                                                       ccAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAg                                           <  3:236947/149‑1 (MQ=255)                                                                                                        cAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATg                                                            <  7:98351/131‑1 (MQ=255)                                                                                                        cAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATg                                                            >  8:98351/1‑131 (MQ=255)                                                                                                        cAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAg                                           >  6:300268/1‑148 (MQ=255)                                                                                                         aaGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTa                                         <  4:129657/149‑1 (MQ=255)                                                                                                           ggTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAAc                                       <  5:300268/149‑1 (MQ=255)                                                                                                            gTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGa                                            <  8:77815/143‑1 (MQ=255)                                                                                                            gTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGa                                            >  7:77815/1‑143 (MQ=255)                                                                                                               aTTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAACGcc                                    <  5:197970/148‑1 (MQ=255)                                                                                                                  ggTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAACGCCgtg                                 <  1:11857/148‑1 (MQ=255)                                                                                                                               tAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAACGCCGTGTCGAC‑CGGCGCGAt                  >  6:79943/1‑149 (MQ=255)                                                                                                                                aGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAACGCCGTGTCGAC‑CGGCGCGATg                 <  7:353738/149‑1 (MQ=255)                                                                                                                                 gCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATg                                                                                  >  4:254802/1‑84 (MQ=255)                                                                                                                                 gCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATg                                                                                  <  3:254802/84‑1 (MQ=255)                                                                                                                                 gCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCACGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAACGCCG‑GTCGACACGgcgttact                 <  2:328070/148‑6 (MQ=255)                                                                                                                                 gCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAACGCCGTGTCGAC‑CGGCGCGATgg                <  5:5819/149‑1 (MQ=255)                                                                                                                                 gCGAATATTGATCGGATCGCC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAACGCCGTGTCGAC‑CGGCGCGATg                 <  3:250266/148‑1 (MQ=255)                                                                                                                                          gatcggatcgGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAACGCCGTGTCGAC‑CGGCGCGATGGCATGTTCAg       >  8:216581/1‑149 (MQ=255)                                                                                                                                               gatcggC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAACGCCGTGTCGAC‑CGGCGCGATGGCATGTTCAGTGAcc  >  3:554/1‑149 (MQ=255)                                                                                                                                               gatcggC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAACGCCGTGTCGAC‑CGGCGCGATGGCATGTTCAGTGAcc  >  6:113724/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 ACGTGCTCGCTGTTGTGCTGGCTGCTGGTGAACTCGAACAGCACCACCGGGCCGTCGGTGGCGACGAACACATAGCGGCTGGGCGAGTGCATCACCCACAAACCAAGGTTATTGGTGTCGGTGGCATAGCGAATATTGATCGGATCGGC‑CAACAGAATGCCGGCGTAATCCTGGGCGCGCAGTTGCTCGCGAATGCGCTCCAGGCGGTAATGGCGCAAGGCTTTGCGATCGATGGCGGCATACGCCTCGAGTAACGCCGTGTCGAC‑CGGCGCGATGGCATGTTCAGTGACC  >  NC_002947/5408370‑5408660

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