Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	850,445	C→G	9.4%	intergenic (‑25/‑780)	pcs ← / → hemA	phosphatidylcholine synthase/glutamyl‑tRNA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	850,445	0	C	G	9.4%	155.5 / 8.0	54	intergenic (‑25/‑780)	pcs/hemA	phosphatidylcholine synthase/glutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base C (27/21);  new base G (2/3);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GTATGAGCATAAATATAGTAGGAAAAATAAATATAGCTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCTCATCAA  >  NC_002947/850302‑850591                                                                                                                                                |                                                                                                                                                    gTATAAGCATAAATATAGTCAGAAAAATAAATATAGCTCAGCGTAGAGGCCGAAATAGAACACAAGTTAACACCCCAAAGACGCACTGAGTGCCTACGAGATACTGCAAAATATTACATTACTCGCGCCTATTATCTAGCTTCCTGCCAt                                                                                                                                               >  2:79470/1‑150 (MQ=255)          tAAATATAGTAGGAAAAATAAATATAGCTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGc                                                                                                                                     <  1:2050568/150‑1 (MQ=255)              tataGTAGGAAAAATAAATATAGCTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATAc                                                                                                                                 <  1:691754/150‑1 (MQ=255)                            taaataTAGCTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAAt                                                                                                                   >  2:1590011/1‑150 (MQ=255)                                 ataGCTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTTAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGttt                                                                                                              <  1:82243/150‑1 (MQ=255)                                     cTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTaa                                                                                                          >  1:1285516/1‑150 (MQ=255)                                          cGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATaa                                                                                                     <  2:2219371/150‑1 (MQ=255)                                              gagGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCtt                                                                                                 <  2:2324838/150‑1 (MQ=255)                                                 gCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGaa                                                                                                                                       >  2:1586626/1‑109 (MQ=255)                                                 gCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGaa                                                                                                                                       <  1:1586626/109‑1 (MQ=255)                                                   cGAAATAGCACACAAGTTAACGACCAAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCGTGCATTGGAAATCGCTTTTAAGATAAGCTTAACac                                                                                            >  2:862081/1‑149 (MQ=255)                                                       aTAGCACACAAGTTAACGCCCCAAATAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACt                                                                                        <  1:79470/150‑1 (MQ=255)                                                          gCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGAATACACTACAATGGAAATCGTTTTTAAGATAAGCATAACAAAACTGaa                                                                                     >  1:1352537/1‑150 (MQ=255)                                                                  gTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATaaaa                                                                             >  2:2102860/1‑150 (MQ=255)                                                                          cccAAAGAAGCACTGAGTGCCTACTAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGcttct                                                                     >  2:1236396/1‑150 (MQ=255)                                                                            cAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCAATG‑CACGAAGCATAACCTGAATTAGAAATCATTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCttg                                                                   >  1:1159027/1‑150 (MQ=255)                                                                               cctatctcTGAGTGCCTACGGGCTACTGCAGCTTGCTCCATTGCTCTCTCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCttgttg                                                                <  2:613080/143‑1 (MQ=255)                                                                               agaagCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCttgttg                                                                <  2:1285516/150‑1 (MQ=255)                                                                                    cACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATATCCTGCATTGGAAATCATTTTTAAGATAAGATTAACAATTCTGATATATAAAAATGCTTCTTGTTGgtgtg                                                           >  1:1905550/1‑150 (MQ=255)                                                                                       tGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTCATAAGACTGAAATATAAAAATGCTTCTTGTTggagatcgg                                                        >  1:613080/1‑143 (MQ=255)                                                                                           tGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATCTGGAAATGCTTCTTGTTGGTGTGTCACCaa                                                    >  1:1756721/1‑150 (MQ=255)                                                                                           tGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCaa                                                    <  1:1164160/150‑1 (MQ=255)                                                                                            gCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCaag                                                   <  1:2102860/150‑1 (MQ=255)                                                                                              cTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGAATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAAAGCTTCTTGTTGGTGTGTCACCaagta                                                 >  1:1522999/1‑150 (MQ=255)                                                                                                 cGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCt                                                                                                  >  1:916236/1‑98 (MQ=255)                                                                                                 cGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCt                                                                                                  <  2:916236/98‑1 (MQ=255)                                                                                                   agaTACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAgtgt                                            <  2:855492/150‑1 (MQ=255)                                                                                                    gaTACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGcttc                                                                      >  2:607148/1‑123 (MQ=255)                                                                                                    gaTACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGcttc                                                                      <  1:607148/123‑1 (MQ=255)                                                                                                     aTACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGc                                          <  2:29371/150‑1 (MQ=255)                                                                                                       aCTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGcat                                        <  1:704871/150‑1 (MQ=255)                                                                                                         tGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGcatca                                      <  1:1236396/150‑1 (MQ=255)                                                                                                          gCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAACGTGGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGcatcat                                     >  1:21869/1‑150 (MQ=255)                                                                                                            aaaTTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTc                                   >  1:1027079/1‑150 (MQ=255)                                                                                                            aaaTTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGAATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAAAGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTc                                   <  2:1522999/150‑1 (MQ=255)                                                                                                                    aaacgatttccaatgcagggtatGCTTCGTGC‑ATGGCAGGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTtcat                           <  1:891846‑M1/127‑1 (MQ=255)                                                                                                                     cATTGCTCGCGCCTATTATATAGATTAATGGCATG‑CAGGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTtcatc                          <  2:21869/150‑1 (MQ=255)                                                                                                                        tGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTaa                                                                                               <  1:1581685/78‑1 (MQ=255)                                                                                                                        tGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTaa                                                                                               >  2:1581685/1‑78 (MQ=255)                                                                                                                          cTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTAc                     >  1:1418175/1‑150 (MQ=255)                                                                                                                           tCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGcttc                                                                      >  2:1163682/1‑100 (MQ=255)                                                                                                                           tCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGcttc                                                                      <  1:1163682/100‑1 (MQ=255)                                                                                                                             ccaatgcagggtatGCTTCGTGC‑ATGGCAGGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTc                                   >  2:754693‑M1/15‑133 (MQ=255)                                                                                                                             ccaatgcagggtatGCTTCGTGC‑ATGGCAGGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCGTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTc                                   <  1:754693‑M1/119‑1 (MQ=255)                                                                                                                               gcCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCgctg                >  2:1363430/1‑150 (MQ=255)                                                                                                                                atgcagggtatGCTTCGTGC‑ATGGCAGGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCgctgc               >  2:891846‑M1/12‑150 (MQ=255)                                                                                                                                    ttatATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTaca           >  2:762315/1‑148 (MQ=255)                                                                                                                                      atataGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCgctgc               >  2:1158885/1‑144 (MQ=255)                                                                                                                                      atataGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCgctgc               <  1:1158885/144‑1 (MQ=255)                                                                                                                                      atataGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAAcc         >  2:1555070/1‑150 (MQ=255)                                                                                                                                       gtatGCTTCGTGC‑ATGGCAGGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCt        <  2:663758‑M1/146‑1 (MQ=255)                                                                                                                                          aGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCtcat     >  2:38586/1‑150 (MQ=255)                                                                                                                                           GCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATGAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCtcatc    >  1:152096‑M1/1‑150 (MQ=255)                                                                                                                                             TTCCTGCCATG‑AACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGAAAAGAATAACAGTACGGAGGTAGAGAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCTCATCaa  >  1:862081‑M1/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                    GTATGAGCATAAATATAGTAGGAAAAATAAATATAGCTCAGCGTAGAGGCCGAAATAGCACACAAGTTAACGCCCCAAAGAAGCACTGAGTGCCTACGAGATACTGCAAATTCTTCCATTGCTCGCGCCTATTATATAGCTTCCTGCCATG‑CACGAAGCATACCCTGCATTGGAAATCGTTTTTAAGATAAGCTTAACAAAACTGAAATATAAAAATGCTTCTTGTTGGTGTGTCACCAAGTAAGTGTGCATCATTCCAGTTCATCATTACCGCTGCTAAACCTCATCAA  >  NC_002947/850302‑850591

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