Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	624,924	T→C	11.1%	D395G (GAT→GGT)	USA300HOU_0559 ←	glycosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	624,924	0	T	C	11.1%	80.8 / 4.1	38	D395G (GAT→GGT)	USA300HOU_0559	glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (23/10);  new base C (2/2);  total (26/12)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.45e-01

ATCTGAAGTAATCAGTTCGTTAGGTCCATATTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTATTTTTGGGACA  >  CP000730/624794‑625061                                                                                                                                   |                                                                                                                                           aTCTGAAGTAATCAGTTCGTTAGGTCCATATTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTAAATTc                                                                                                                                  >  2:224492/1‑141 (MQ=255)         ttataaGTTCGTTAGGTACATATTTAATATCATAACTAATAAATGCAAAACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGTATCGCTATATTCT‑TGGTCTa                                                                                                                         <  1:268421/136‑1 (MQ=255)             cAGTTCGTTAGGTCCATATTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGAtt                                                                                                                     >  1:94681/1‑141 (MQ=255)                 tCGTTAGGTCCATATTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCtt                                                                                                                 >  2:5268/1‑141 (MQ=255)                      aGGTCCATATTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATAGCTATATTCT‑TGGTCTAGATTGCTTaaaaa                                                                                                            >  2:241951/1‑141 (MQ=255)                             tatTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCACCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCgt                                                                                                     >  2:25637/1‑141 (MQ=255)                             tatTTAATATCATAACTAATAACTGGAAAACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCACCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCgt                                                                                                     >  1:378317/1‑141 (MQ=255)                                 tAATATCATAACTAATAACTGGCTGACCAAGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCCCTATATTGT‑TGGTCTAGATTGCTTAAAAAGACGAgtaatt                                                                                                 >  1:185638/1‑139 (MQ=255)                                 tAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCgtaagt                                                                                                 >  1:107106/1‑141 (MQ=255)                                 tAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCgtaagt                                                                                                 >  1:165386/1‑141 (MQ=255)                                        aaaaCTAATAACTCGAACACCATGTGTAAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATg                                                                                          <  2:179640/139‑1 (MQ=255)                                        aTAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATg                                                                                          <  2:327668/141‑1 (MQ=255)                                            cTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGAttt                                                                                      <  1:335714/141‑1 (MQ=255)                                                   cTGGAACACCATGTGCTAAAGATTAAATTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATc                                                                               >  2:271389/1‑141 (MQ=255)                                                        acacCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCatat                                                                          >  1:289683/1‑141 (MQ=255)                                                                  cTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCa                                                                >  1:312919/1‑141 (MQ=255)                                                                  cTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAACCGAGTACGTAAACATGATTTTCCAGTCCaatatgatgtatca                                                                >  2:262577/1‑128 (MQ=255)                                                                         ttCAAGTAGCGTTAAAGAGCAACCTACCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCAtt                                                         <  1:86199/141‑1 (MQ=255)                                                                          tCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCAttt                                                        >  2:350191/1‑141 (MQ=255)                                                                          tCAAGTAGCGATAAAGAGAAACCTGACATGCTACTTGTGATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAACTCAttt                                                        >  2:88412/1‑141 (MQ=255)                                                                           cAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAAATCAAATAGGCATCACTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCACTCCATATTTTTGTATCAATTCATTTa                                                       >  1:233113/1‑141 (MQ=255)                                                                           cAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATAAAAATCAAATAAACAACAC‑ATTTGGTGTTGTCTGGATTAGTAAAAAAACCACGTACGTAAACCTGGATTTCCAAGCCATAATTTTGTAGCACTTCAttga                                                       >  2:95754/1‑139 (MQ=255)                                                                            aaGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTaa                                                      >  2:55313/1‑141 (MQ=255)                                                                            aaGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCAGCGTAAGTAAACATGATTTTCCAATCCATATTGTTGTATCAATTCTTTTaa                                                      >  2:286535/1‑141 (MQ=255)                                                                                aGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAAtttt                                                  >  1:348339/1‑141 (MQ=255)                                                                                     tacagcgaAACCTTCCATGTTACTTGTTATTAAACTCAAACAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTCATTTTTTACt                                             <  1:56047/135‑1 (MQ=255)                                                                                      aaaGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCACTGCATTTGATTTGTTCCtt                                            >  1:33916/1‑141 (MQ=255)                                                                                       aaGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATGAATTAATTTAATTTTTTACttt                                           >  2:211782/1‑141 (MQ=255)                                                                                       aaGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGCCATCGCTATATTCT‑TGGTCTAGCTTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCACTTCATTTAATTTTTTACttt                                           >  2:45583/1‑141 (MQ=255)                                                                                                      tGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACcat                            <  1:288671/141‑1 (MQ=255)                                                                                                                   taAACTCAAATAGGCACCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGctc               <  1:25637/141‑1 (MQ=255)                                                                                                                   taAACTCAAATAGGCACCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGctc               <  2:378317/141‑1 (MQ=255)                                                                                                                      aCTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTAt            >  1:376015/1‑141 (MQ=255)                                                                                                                       cTCAAATATGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTAtt           <  1:1666/141‑1 (MQ=255)                                                                                                                       cTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTAtt           >  2:97685/1‑141 (MQ=255)                                                                                                                          aaaTAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTAttttt        <  1:154915/141‑1 (MQ=255)                                                                                                                             tAGGCATCGCCATATTCG‑TGGTCTAGATTGCTCAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGATGCTCTATTTTTggg     <  2:3522/141‑1 (MQ=255)                                                                                                                                gCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTATTTTTGGGACa  >  1:376547/1‑141 (MQ=255)                                                                                                                                   |                                                                                                                                           ATCTGAAGTAATCAGTTCGTTAGGTCCATATTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTATTTTTGGGACA  >  CP000730/624794‑625061

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