Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	369,640	T→C	8.6%	G240G (GGT→GGC)	PP_0308 →	dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	369,640	0	T	C	8.6%	124.9 / 2.0	46	G240G (GGT→GGC)	PP_0308	dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (19/23);  new base C (2/2);  total (21/25)
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.20e-01

ACAAGTCGGACGAAGAGCTGAAGTTCATCGCCGATCACGGCGGTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCTGACCAACTTCGGCAAGATCATCAACCCGCT  >  NC_002947/369495‑369765                                                                                                                                                  |                                                                                                                              acaaGTCGGACGAAGAGCTGAAGTTCATCGCCGATCACGGCGGTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTg                                                                                                                              >  8:102642/1‑147 (MQ=255)     gTCGGACGAAGAGCTGAAGTTCATCGCCGATCACGGCGGTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGaa                                                                                                                            <  6:282245/145‑1 (MQ=255)     gTCGGACGAAGAGCTGAAGTTCATCGCCGATCACGGCGGTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGaa                                                                                                                            >  5:282245/1‑145 (MQ=255)        ggACGAAGAGCTGAAGTTCATCGCCGATCACGGCGGTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTTCACCATGAACATCGTCGGTGAAGACGcc                                                                                                                      <  1:203104/148‑1 (MQ=255)          aCGAAGAGCTGAAGTTCATCGCCGATCACGGCGGTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCatc                                                                                                                   <  8:258538/149‑1 (MQ=255)               gagCTGAAGTTCATCGCCGATCACGGCGGTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCatcggtat                                                                                                              <  7:102642/149‑1 (MQ=255)                     aaGTTCATCGCCGATCACGGCGGTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCatcggtatcggt                                                                                                          <  5:141299/147‑1 (MQ=255)                      aGTTCATCGCCGATCACGGCGGTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTAcc                                                                                                       >  6:147600/1‑149 (MQ=255)                        ttCATCGCCGATCACGGCGGTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGa                                                                                                     >  2:248917/1‑149 (MQ=255)                                       ggcggTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCa                                                                                      >  8:73005/1‑149 (MQ=255)                                        gcggTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCAc                                                                                     <  2:9582/149‑1 (MQ=255)                                            tttCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGc                                                                                  >  6:251770/1‑148 (MQ=255)                                                    gCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGActtc                                                                         >  8:333211/1‑149 (MQ=255)                                                      gTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGActtctt                                                                       <  7:289348/149‑1 (MQ=255)                                                       tCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACttcttc                                                                      <  7:73005/149‑1 (MQ=255)                                                         aCCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACTATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCg                                                                     >  5:47518/1‑148 (MQ=255)                                                          ccATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGaa                                                                   <  2:99314/149‑1 (MQ=255)                                                           cATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGaa                                                                   >  1:153716/1‑148 (MQ=255)                                                                  gcgcCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGAc                                                           <  1:248917/149‑1 (MQ=255)                                                                  gcgcCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGAc                                                           >  7:41007/1‑149 (MQ=255)                                                                           cTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACaa                                                  <  6:77798/149‑1 (MQ=255)                                                                            tGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGGACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAg                                                 >  8:333957/1‑149 (MQ=255)                                                                              gCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGCGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGActtctt                                                                       >  6:52580/1‑125 (MQ=255)                                                                              gCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGCGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGActtctt                                                                       <  5:52580/125‑1 (MQ=255)                                                                              gCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGCGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGActtctt                                                                       <  5:47231/125‑1 (MQ=255)                                                                              gCCAAAGGCATTGACTCGAACATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGCGAAGACGCCATCGGTATCGGTAACGACTTCACCCAGGGTACCGGCCAGGActtctt                                                                       >  6:47231/1‑125 (MQ=255)                                                                                 aaaGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCt                                             >  1:206286/1‑148 (MQ=255)                                                                                 aaaGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCt                                             <  2:206286/148‑1 (MQ=255)                                                                                           aCTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACTATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGt                                  <  6:47518/149‑1 (MQ=255)                                                                                           aCTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTAcg                                          <  8:188021/141‑1 (MQ=255)                                                                                           aCTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTAcg                                          >  7:188021/1‑141 (MQ=255)                                                                                           aCTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCccgccg                                   <  5:147600/148‑1 (MQ=255)                                                                                            cTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTc                                 <  5:17010/149‑1 (MQ=255)                                                                                             tCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCt                                >  5:371351/1‑149 (MQ=255)                                                                                               gACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGGACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCt                                <  7:333957/147‑1 (MQ=255)                                                                                               gACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCt                                <  5:342189/147‑1 (MQ=255)                                                                                                aCCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCTGAc                             <  6:371351/149‑1 (MQ=255)                                                                                                 ccATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCTGAcc                            <  5:117053/149‑1 (MQ=255)                                                                                                  cATCGATGACTACGCCGAGGACATCGAGTACACCATGAACATCGTCGGTGAAGAAGACATCGGTATCGGTACCGACTTCACCCAGGGTCACGGACATGACTTCTTCGAATACCTGAACCCTTACAAGGTCTACGCCCGACGTCTgcac                            >  4:229922/1‑145 (MQ=255)                                                                                                            tACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCTGACCAACTTCGGc                   >  4:271392/1‑147 (MQ=255)                                                                                                                 cGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCTGACCAACTTCGGCAAGatca            <  5:251770/149‑1 (MQ=255)                                                                                                                        aTCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCTGACCAACTTCGGCAAGATCATCAcccc     >  5:303173/1‑149 (MQ=255)                                                                                                                         tCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCTGACCAACTTCGGCAAGATCATCAACCCg    >  5:28980/1‑149 (MQ=255)                                                                                                                           gAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCTGACCAACTTCGGCAAGATCATCAACCCGCt  <  7:333211/149‑1 (MQ=255)                                                                                                                                             gtcgGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCTGACCAACTTCGGCAAg                >  2:137564/1‑117 (MQ=255)                                                                                                                                             gtcgGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCTGACCAACTTCGGCAAg                <  1:137564/117‑1 (MQ=255)                                                                                                                                                  |                                                                                                                              ACAAGTCGGACGAAGAGCTGAAGTTCATCGCCGATCACGGCGGTTTCGTCGGCGTCACCATGTTCGCGCCGTTCCTGGCCAAAGGCATTGACTCGACCATCGATGACTACGCCGAGGCCATCGAGTACACCATGAACATCGTCGGTGAAGACGCCATCGGTATCGGTACCGACTTCACCCAGGGTCACGGCCAGGACTTCTTCGAATACCTGACCCATGACAAGGGCTACGCCCGCCGTCTGACCAACTTCGGCAAGATCATCAACCCGCT  >  NC_002947/369495‑369765

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