Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,089,214	T→C	11.1%	Y131Y (TAT→TAC)	potF‑IV →	putrescine‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,089,214	0	T	C	11.1%	117.3 / 5.7	45	Y131Y (TAT→TAC)	potF‑IV	putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (20/20);  new base C (3/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTGTTCCTCGCCAAGCAGATCGAAGCCGGGGTGTTCCAGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTC  >  NC_002947/6089072‑6089355                                                                                                                                               |                                                                                                                                              aCGTGTTCCTCGCCAAGCAGATCGAAGCCGGGGTGTTCCAGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTAc                                                                                                                                         >  5:166015/1‑149 (MQ=255)        ccTCGCCAAGCAGATCGAAGCCGGGGTGTTCCAGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCg                                                                                                                                  >  8:161725/1‑149 (MQ=255)            gCCAAGCAGATCGAAGCCGGGGTGTTCCAGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCt                                                                                                                              <  7:161725/149‑1 (MQ=255)                   aGATCGAAGCCGGGGTGTTCCAGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGc                                                                                                                       <  4:66124/149‑1 (MQ=255)                        gAAGCCGGGGTGTTCCAGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCaa                                                                                                                  <  6:137638/149‑1 (MQ=255)                           gCCGGGGTGTTCCAGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAAccc                                                                                                               >  6:201809/1‑149 (MQ=255)                            ccGGGGTGTTCCAGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCg                                                                                                              <  6:102767/149‑1 (MQ=255)                             cGGGGTGTTCCAGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCgg                                                                                                             <  1:35905/149‑1 (MQ=255)                                    ttCCAGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGt                                                                                                      >  5:206916/1‑149 (MQ=255)                                        aGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTACATATACGGCAcg                                                                                                                                   >  3:257795/1‑115 (MQ=255)                                        aGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTACATATACGGCAcg                                                                                                                                   <  4:257795/116‑2 (MQ=255)                                                      agcaAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGt                                                                                                                                           <  7:169889/94‑1 (MQ=255)                                                      agcaAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGt                                                                                                                                           >  8:169889/1‑94 (MQ=255)                                                         aaGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACaa                                                                                 <  6:166015/149‑1 (MQ=255)                                                               ccGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGcc                                                                           >  8:117087/1‑149 (MQ=255)                                                                  aaCTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGt                                                                        <  6:206916/149‑1 (MQ=255)                                                                       gCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGa                                                                     >  1:244470/1‑147 (MQ=255)                                                                        cAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAg                                                                  >  6:242065/1‑149 (MQ=255)                                                                         aaCACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGc                                                                 >  6:262055/1‑149 (MQ=255)                                                                           cacCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTg                                                               <  6:31601/149‑1 (MQ=255)                                                                              cTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCaa                                                                                                   <  5:97939/110‑1 (MQ=255)                                                                              cTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCaa                                                                                                   >  6:97939/1‑110 (MQ=255)                                                                              cTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAgg                                                                                                 >  7:103377/1‑112 (MQ=255)                                                                              cTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAgg                                                                                                 <  8:103377/112‑1 (MQ=255)                                                                                 gACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGAcc                                                          <  6:15096/148‑1 (MQ=255)                                                                                 gACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGAc                                                           <  5:242065/147‑1 (MQ=255)                                                                                    ccGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGAcc                                                          <  5:18038/145‑1 (MQ=255)                                                                                    ccGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGAcc                                                          >  6:18038/1‑145 (MQ=255)                                                                                      gtcACTGATGAATCTGATCGAAGCTAAGGACCCTGGCAACTAGTATGCCGTGCCGTATATGTACGGCACCGTGCTGATCTTCTTGAACCGGGCCAAGGTGAAGGCTGCACTCGGCTAGTACGCGCTGGTGGAGAGCTGGGTCGTGATCt                                                    <  8:199327/147‑1 (MQ=255)                                                                                          cTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCaa                                                <  2:244470/149‑1 (MQ=255)                                                                                           tgatgaAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAg                                               >  4:211946/1‑149 (MQ=255)                                                                                              tgaAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGaa                                            <  5:262055/149‑1 (MQ=255)                                                                                                  gCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAagag                                         >  6:25566/1‑148 (MQ=255)                                                                                                    tGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGaa                                       >  4:67549/1‑148 (MQ=255)                                                                                                                 cGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGcc                                <  3:101318/142‑1 (MQ=255)                                                                                                                 cGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGcc                                >  4:101318/1‑142 (MQ=255)                                                                                                                   aCCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGaa                                            >  8:4443/1‑128 (MQ=255)                                                                                                                   aCCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGaa                                            <  7:4443/128‑1 (MQ=255)                                                                                                                   aCCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAg                       >  3:145354/1‑149 (MQ=255)                                                                                                                    cccGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGc                      <  5:25566/149‑1 (MQ=255)                                                                                                                           aacaaGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCACTgcgg               >  1:226956/1‑149 (MQ=255)                                                                                                                             caaGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTgcggcg             >  1:136030/1‑149 (MQ=255)                                                                                                                                    acggtgccgtACATATACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCt      >  7:137867/4‑149 (MQ=255)                                                                                                                                      ggtgccgtACATATACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGc    >  6:94798/2‑149 (MQ=255)                                                                                                                                        tgccgtACATATACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTc  <  5:94798/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                              ACGTGTTCCTCGCCAAGCAGATCGAAGCCGGGGTGTTCCAGCCGCTGGATCGCAGCAAGCTGCCGAACTGGCAACACCTGGACCCGGCACTGATGAAGCTGATCGAAGCTAACGACCCGGGCAACAAGTTTGCCGTGCCGTATATGTACGGCACCGTGCTGATCGGCTTCAACCCGGCCAAGGTCAAGGCTGCACTCGGCGACAACGCGCCGGTGGACAGCTGGGACCTGATCTTCAAGGAAGAGAACATCGCCAAGCTCAAGCAGTGCGGCGTGGCCCTGCTC  >  NC_002947/6089072‑6089355

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