Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,315,536	G→A	8.3%	L153L (CTC→CTT)	PP_2037 ←	aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,315,536	0	G	A	8.3%	128.2 / 3.2	48	L153L (CTC→CTT)	PP_2037	aldolase
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base A (2/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

GCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGT  >  NC_002947/2315391‑2315681                                                                                                                                                  |                                                                                                                                                  gcagAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGa                                                                                                                                               >  1:659463/1‑150 (MQ=255)   agagCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAAt                                                                                                                                             <  2:248310/150‑1 (MQ=255)      gcgagcACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGca                                                                                                                                          <  1:2850/150‑1 (MQ=255)                  gCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGa                                                                                                                                                  <  2:326635/130‑1 (MQ=255)                  gCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGa                                                                                                                                                  >  1:326635/1‑130 (MQ=255)                  gCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAAc                                                                                                                              >  1:127958/1‑150 (MQ=255)                         cgatcgaACGCCCAGCGATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGc                                                                                                                       >  2:395998/1‑150 (MQ=255)                                  gCCCAGCGATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAAc                                                                                                                              <  2:303673/134‑1 (MQ=255)                                  gCCCAGCGATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAAc                                                                                                                              >  1:303673/1‑134 (MQ=255)                                    ccAGCGATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCAtg                                                                                                            >  1:1043946/1‑150 (MQ=255)                                         gATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCga                                                                                                       <  2:1106320/150‑1 (MQ=255)                                              gcagGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAg                                                                                                                                                 >  2:591417/1‑103 (MQ=255)                                              gcagGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAg                                                                                                                                                 <  1:591417/103‑1 (MQ=255)                                              gcagGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATgg                                                                                                  <  2:1043946/150‑1 (MQ=255)                                                 ggCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCa                                                                                               >  1:153858/1‑150 (MQ=255)                                                  gCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATg                                                                                                   <  2:1251395/145‑1 (MQ=255)                                                  gCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATg                                                                                                   >  1:1251395/1‑145 (MQ=255)                                                    cATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCg                                                                                            >  1:575717/1‑150 (MQ=255)                                                      tgatgGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGa                                                                                          >  1:229301/1‑150 (MQ=255)                                                       gatgGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGaa                                                                                         >  1:156107/1‑150 (MQ=255)                                                            gAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTcc                                                                                    >  1:1072313/1‑150 (MQ=255)                                                             aaaGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTCAAGGAATg                                                                                                                                            >  2:1271317/1‑93 (MQ=255)                                                             aaaGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTCAAGGAATg                                                                                                                                            <  1:1271317/93‑1 (MQ=255)                                                                 cAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCAt                                                                               <  2:1091156/150‑1 (MQ=255)                                                                  aGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTTGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATc                                                                              >  1:110050/1‑150 (MQ=255)                                                                        gcgcgcTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAg                                                                        <  2:127958/150‑1 (MQ=255)                                                                                 tCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCtg                                                               >  2:1206668/1‑150 (MQ=255)                                                                                  cGCCAATCGCGTTCGAGATGCTTTCGCCTAGTTATGAGCCCGCCTGGACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCtgt                                                              <  2:782068/150‑1 (MQ=255)                                                                                      aaTCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAg                                                          >  2:233191/1‑150 (MQ=255)                                                                                             ttCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCgggtgg                                                   <  2:156107/150‑1 (MQ=255)                                                                                                  gatgatTTCGCCTTCTTCGTTGCCCACCGGTACCGGGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGtgtg                                              >  1:567418/1‑150 (MQ=255)                                                                                                  gatgatTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTTTGCAGCGGGTGGGtgtg                                              >  1:389907/1‑150 (MQ=255)                                                                                                   atgatTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTTGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTgg                                             <  2:110050/150‑1 (MQ=255)                                                                                                   atgatTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTgg                                             <  1:1232864/150‑1 (MQ=255)                                                                                                     gatTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGgat                                           <  2:957795/150‑1 (MQ=255)                                                                                                         tCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGatgatg                                       <  1:1135168/150‑1 (MQ=255)                                                                                                              ttcttcGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGtt                                  <  1:395998/150‑1 (MQ=255)                                                                                                                     ttGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCg                           >  1:341358/1‑150 (MQ=255)                                                                                                                         ccACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGgggc                       >  1:1064691/1‑150 (MQ=255)                                                                                                                                 tcacCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAATTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCgg               <  2:567418/149‑1 (MQ=255)                                                                                                                                      tgtgcATTCCTCAAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGgatga                                         >  1:533098/5‑119 (MQ=255)                                                                                                                                      tgtgcATTCCTCAAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGgatga                                         <  2:533098/115‑1 (MQ=255)                                                                                                                                         ccATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGCCGCGGTAcagcca       >  2:1245981/1‑150 (MQ=255)                                                                                                                                          cATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTAcagccag      <  2:389907/150‑1 (MQ=255)                                                                                                                                          cATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTAcagccag      >  1:1216716/1‑150 (MQ=255)                                                                                                                                            ttCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCt    >  1:714187/1‑150 (MQ=255)                                                                                                                                              ccTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCtgt  <  1:1206668/150‑1 (MQ=255)                                                                                                                                              ccTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCtgt  <  2:659463/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCAGCAGGCCATGATGGGAAAGCAGGATTGCGCGCTTGTCGCCAATCGCGTTCGCGATGATTTCGCCTTCTTCGTTGCCCACCGGTACACCGGGCCATTCCTTGAGGAATGCACAGTCGTCGAACAGGGGGCACAGGTCCATGTGCGAGATGGCCAGCGGAACTTCCAGCATCGACAAGGCGGCGCTGTGCAGCGGGTGGGTGTGGATGATGCAGTTGACGTCGGGGCGGGCGCGGTACAGCCAGCTGT  >  NC_002947/2315391‑2315681

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