Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,305,233	T→C	7.1%	F194F (TTT→TTC)	ymfH →	M16C subfamily peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,305,233	0	T	C	7.1%	157.6 / 3.5	56	F194F (TTT→TTC)	ymfH	M16C subfamily peptidase
Reads supporting (aligned to +/- strand):  ref base T (21/31);  new base C (2/2);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

CAACCTGGATATAAATTAATGTTTGATACATTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAACC  >  CP000730/1305096‑1305370                                                                                                                                          |                                                                                                                                          cAACCTGGATATAAATTAATGTTTGATACATTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTat                                                                                                                                        >  1:552909/1‑141 (MQ=255)     cTGGATATAAATTAATGTTTGATACATTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTatcatc                                                                                                                                    >  1:488214/1‑141 (MQ=255)     cTGGATATAAATTAATGTTTGATACATTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTCTatcatc                                                                                                                                    >  2:103433/1‑141 (MQ=255)     cTGGATATAAATTAATGTTTGATACATTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTCTatcatc                                                                                                                                    >  2:103438/1‑141 (MQ=255)      tgaatagAAATTAATGTTTGAGACATTACGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATcatcc                                                                                                                                   <  2:671008/138‑1 (MQ=255)           ataAATTAATGTTTGATACATTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCaa                                                                                                                              <  1:106760/141‑1 (MQ=255)           ataAATTAATGTTTGATACATTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCaa                                                                                                                              <  2:47192/141‑1 (MQ=255)             aaaTTAATGTTTGATACATTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAat                                                                                                                            <  1:102584/141‑1 (MQ=255)                   aTGTTTGATACATTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGtt                                                                                                                      >  2:26606/1‑141 (MQ=255)                            aCATTGCTTGCAATGTATCAACAACATCCAATACGTGTTGATATGGCCGATAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATttgtt                                                                                                             <  1:527251/141‑1 (MQ=255)                            aCATTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATttgtt                                                                                                             <  1:517700/141‑1 (MQ=255)                            aCATTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATttgtt                                                                                                             <  2:672778/141‑1 (MQ=255)                              aTTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATttgttgt                                                                                                           >  1:497389/1‑141 (MQ=255)                               ttGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATttgttgtt                                                                                                          >  2:129509/1‑141 (MQ=255)                                  cacgCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGc                                                                                                       >  1:22742/4‑141 (MQ=255)                                  cGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGc                                                                                                       >  2:116055/1‑141 (MQ=255)                                             tCAAAAACATCCAATACGTGTAGATATTGACGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATcc                                                                                            <  2:166120/141‑1 (MQ=255)                                                       ccAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAata                                                                                  <  1:706754/141‑1 (MQ=255)                                                         aaTACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGTGATGTGGATCCTGAAGCAATATg                                                                                >  2:319781/1‑141 (MQ=255)                                                             cGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGa                                                                            >  1:143606/1‑141 (MQ=255)                                                               tgtTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAAt                                                                          >  2:257021/1‑141 (MQ=255)                                                                    atatTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTaa                                                                     <  1:295645/141‑1 (MQ=255)                                                                       ttGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAaac                                                                  <  1:193505/141‑1 (MQ=255)                                                                          ccGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAaacaac                                                               <  1:54724/141‑1 (MQ=255)                                                                             gTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACg                                                            <  2:411512/141‑1 (MQ=255)                                                                                   tAGGAAGTATACACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCAATCAAATATGGTTTTCTTTGTTGTTGCCGCTGTGGTTCCTGACGCAATATGTCGAATAGAAAAACAACACCAGGATg                                                      >  2:37202/1‑141 (MQ=255)                                                                                   tAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATg                                                      >  1:616158/1‑141 (MQ=255)                                                                                    aGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGc                                                     <  1:282912/141‑1 (MQ=255)                                                                                    aGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGc                                                     <  1:422032/141‑1 (MQ=255)                                                                                     gAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCt                                                    <  2:669970/141‑1 (MQ=255)                                                                                          tatataCGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGtaa                                               <  2:22742/141‑1 (MQ=255)                                                                                          tatataCGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTAGCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGtaa                                               <  1:183325/141‑1 (MQ=255)                                                                                           atataCGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGtaat                                              >  2:203358/1‑141 (MQ=255)                                                                                           atataCGATATTAAAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATAGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGtaat                                              >  2:559284/1‑141 (MQ=255)                                                                                            tataCGATATTACATAAGAGGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACATCACGAGGATGCTCGtaata                                             <  2:277995/141‑1 (MQ=255)                                                                                               aCGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAg                                          <  1:438808/141‑1 (MQ=255)                                                                                               aCGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAg                                          >  2:301901/1‑141 (MQ=255)                                                                                               aCGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAg                                          >  2:301786/1‑141 (MQ=255)                                                                                                        cAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTaaccaac                                 >  1:99539/1‑141 (MQ=255)                                                                                                          aaaGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAAccc                               <  1:26606/141‑1 (MQ=255)                                                                                                          aaaGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAAccc                               >  1:467965/1‑141 (MQ=255)                                                                                                              atgatTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCaaca                           <  2:190562/141‑3 (MQ=255)                                                                                                                gatTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATc                         <  1:194978/141‑1 (MQ=255)                                                                                                                      tctttATGTTATGAAACGTTTTGTGACACATGAACGATGGTATTATTTGTTGTTGGCTATGTGGATCATGACGCATTATGTCGAATAGTAAAACAAAAGGAGGATGTTCGTAATAAAGTTGACCGACCCAAAGTCGGACGa                   <  2:305014/137‑1 (MQ=255)                                                                                                                      tATCTATGTTATGAAACGTTTTATCATCCATCAAATATGCTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGa                   <  1:603564/141‑1 (MQ=255)                                                                                                                        tCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAgg                 <  1:352329/141‑1 (MQ=255)                                                                                                                        tCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAgg                 <  2:31570/141‑1 (MQ=255)                                                                                                                         cTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCGTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGa                <  2:681273/141‑1 (MQ=255)                                                                                                                             gTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACttg            <  2:463584/141‑1 (MQ=255)                                                                                                                             gTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACttg            <  1:319329/141‑1 (MQ=255)                                                                                                                                 tGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTtga        <  1:244095/141‑1 (MQ=255)                                                                                                                                 tGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTtga        >  2:396039/1‑141 (MQ=255)                                                                                                                                    aaCGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTtgatga     >  2:479094/1‑141 (MQ=255)                                                                                                                                      cGTTCTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAAc   <  1:103438/141‑1 (MQ=255)                                                                                                                                      cGTTCTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAAc   <  1:103433/141‑1 (MQ=255)                                                                                                                                       gTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAAcc  >  1:613154/1‑141 (MQ=255)                                                                                                                                          |                                                                                                                                          CAACCTGGATATAAATTAATGTTTGATACATTGCGTGCAATGTATCAACAACATCCAATACGTGTTGATATTGCCGGTAGTGTAGAAAGTATATACGATATTACAAAAGATGATTTGTATCTATGTTATGAAACGTTTTATCATCCATCAAATATGGTTTTATTTGTTGTTGGCGATGTGGATCCTGAAGCAATATGTCGAATAGTAAAACAACACGAGGATGCTCGTAATAAAGTTAACCAACCCAAAATCGAACGAGGACTTGTTGATGAACC  >  CP000730/1305096‑1305370

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