Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,383,254	A→T	10.6%	L16Q (CTG→CAG)	PP_1204 ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,383,254	0	A	T	10.6%	112.8 / 2.9	45	L16Q (CTG→CAG)	PP_1204	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (17/23);  new base T (3/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTG  >  NC_002947/1383112‑1383396                                                                                                                                               |                                                                                                                                               tCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTgcgg                                                                                                                                          <  5:294565/149‑1 (MQ=255)  cGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTgcggc                                                                                                                                         >  6:316531/1‑149 (MQ=255)    cTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTCCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACTGTGCGgcgc                                                                                                                                       <  5:273402/149‑1 (MQ=255)       cGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGcc                                                                                                                                      <  8:25763/147‑1 (MQ=255)              tCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAggg                                                                                                                              >  4:127312/1‑148 (MQ=255)                cATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCCGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTa                                                                                                                            >  8:278848/1‑148 (MQ=255)                 aTCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTaaa                                                                                                                          <  5:90297/149‑1 (MQ=255)                  tCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAg                                                                                                                         <  7:278848/149‑1 (MQ=255)                   cACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAgg                                                                                                                        >  1:88380/1‑149 (MQ=255)                      cTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGcaca                                                                                                                                                >  4:73320/1‑122 (MQ=255)                      cTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGcaca                                                                                                                                                <  3:73320/122‑1 (MQ=255)                      cTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTAc                                                                                                                     <  3:63174/149‑1 (MQ=255)                         cGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACtt                                                                                                                   <  5:155296/148‑1 (MQ=255)                            cTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCt                                                                                                               <  5:300307/149‑1 (MQ=255)                             tagCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGAAGAGCAGGGCGCCGCACAGTGCCGCGCCCCCAAGGGTAAAGGTACTTAGCtt                                                                                                              >  5:265322/3‑149 (MQ=255)                               gCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGt                                                                                                            <  7:218199/149‑1 (MQ=255)                                    cATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATg                                                                                                       <  8:179982/149‑1 (MQ=255)                                           ctcccgggTCTGCTCGCGCAGTTGGCGCTCGACGCGCGCCTCATGGTCGCGGCGCTGGCACAGCTGGGTGGCCCAGCCGCTGCCCCGCAGGCCGCCGCCCAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTc                                                                                                <  2:129734/144‑1 (MQ=255)                                             gCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGt                                                                                              >  5:317604/1‑149 (MQ=255)                                                      gCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGt                                                                                      >  5:296616/1‑148 (MQ=255)                                                       cTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCAGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTa                                                                                    >  3:118571/1‑149 (MQ=255)                                                          gAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCa                                                                                 <  2:88380/149‑1 (MQ=255)                                                            gCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATTTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAgc                                                                               <  2:108868/149‑1 (MQ=255)                                                             cAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAgcg                                                                              <  1:136942/149‑1 (MQ=255)                                                                tttGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTg                                                                                        <  4:215232/136‑1 (MQ=255)                                                                tttGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTg                                                                                        >  3:215232/1‑136 (MQ=255)                                                                     gcTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGATcgc                                                                      >  1:170571/1‑149 (MQ=255)                                                                       tCGGCGCGCTCCTCATGGTCGCTGCGCTGGTAGTTTTGGTTGGCGCAGCCGGTGAGGTGCCGGGCGCCGCTCAGTGGGGGGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTTCATTGGTTATCAGCGGCGAACGCgg                                                                    <  6:34577/149‑1 (MQ=255)                                                                       tCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCCAGGGTAAAGGTACTTAGCTTGTACAAGATAACTCGGGCATTGGTTATAACCGGCGCACGCgg                                                                    >  2:286043/1‑149 (MQ=255)                                                                         gACGCGCTCCTCATGGGCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCgg                                                                    <  1:286043/147‑1 (MQ=255)                                                                         gACGCGCTCCGCAGGGTCGCTGCGCGGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGAAAACGTACTTCGCGTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCgg                                                                    <  6:266401/147‑1 (MQ=255)                                                                                     aTGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAg                                                        <  2:166532/147‑1 (MQ=255)                                                                                     aTGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGt                                                      <  5:316531/149‑1 (MQ=255)                                                                                       ggTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGt                                                      >  8:166455/1‑147 (MQ=255)                                                                                                 ctgGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCAtt                                          <  6:296616/149‑1 (MQ=255)                                                                                                 ctgGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCAtt                                          <  7:166455/149‑1 (MQ=255)                                                                                                   gggACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCAtt                                          >  7:143303/1‑147 (MQ=255)                                                                                                      aCATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTTAAGGTACTTCGCTTGTACATGAGCCCTCGTGCATTGGTTATCAGCGTGGACCGCGGGCCGGGAGGCAGGTCAGCCCGTCAtgtcgcg                                     >  6:230003/1‑143 (MQ=255)                                                                                                        aTCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCgg                                    >  7:118110/1‑148 (MQ=255)                                                                                                                 gcgcAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTgg                          <  8:118110/149‑1 (MQ=255)                                                                                                                      gCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTc                                                     <  4:121012/117‑1 (MQ=255)                                                                                                                      gCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACGTCGCGTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTc                                                     >  3:121012/1‑117 (MQ=255)                                                                                                                            gACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTc               >  7:240056/1‑149 (MQ=255)                                                                                                                                   ggggCGCCGCACTGGGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCAAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcggtgcgg        <  5:271625/148‑1 (MQ=255)                                                                                                                                    gggCGCCGCACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGAGGCAACGCCTGGGCGTCCTGCTcggtgcggt       >  8:111352/1‑149 (MQ=255)                                                                                                                                    gggCGCCGCACTGTGCGGCCCCCCCAAGGGTAAAGGTACTTAGCATGTACATGATAACACGAGCATAGGTTATAAGCCGCCAACGCGGGCCAGCACGAAGGAGAGCACTTCAATTAGCGGCAACGCCGGGCCGTCCAGCAcggggccgt       >  2:152920/1‑149 (MQ=255)                                                                                                                                         ccgcACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTg  >  2:195364/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                               TCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTG  >  NC_002947/1383112‑1383396

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