Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,151,090	A→G	9.9%	K66K (AAA→AAG)	USA300HOU_1076 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,151,090	0	A	G	9.9%	104.1 / 3.9	41	K66K (AAA→AAG)	USA300HOU_1076	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/10);  new base G (2/2);  total (29/12)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

CAGAGCACATACGATATGTAGCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCAGACGAGGTTTTTG  >  CP000730/1150966‑1151230                                                                                                                              |                                                                                                                                             cAGAGCACATACGATATGTAGCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATCAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAATAATAAAACTagaagaaa                                                                                                                              >  2:445675/1‑140 (MQ=255)        caTACGATATGTAGCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAGGTACTACTAGAAGAAGacaacca                                                                                                                       >  2:380413/1‑138 (MQ=255)        caTACGATATGTAGCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAGGTACTACTAGAAGAAGAAAATCg                                                                                                                       >  2:444432/1‑141 (MQ=255)             aatatGTAGCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGtt                                                                                                                  <  1:552659/140‑1 (MQ=255)                     gCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAaacaa                                                                                                          >  2:264782/1‑141 (MQ=255)                       acaTTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACGACTAGAAGAAGAAAATCGACGTTGGAAaaaaca                                                                                                        >  2:162959/1‑141 (MQ=255)                             aGTTGATGATAAAATAAAAGAATTAGGGCATAAAGCAGCAGGTTTAGATAATTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAAACGACGTTTAAAACAACAACTTca                                                                                                  >  2:149926/1‑141 (MQ=255)                                      aaaaagaaaagAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTgca                                                                                         <  1:203261/140‑1 (MQ=255)                                       aaaataaaaGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTgcag                                                                                        <  1:614097/141‑1 (MQ=255)                                             agaggatgagagTATAAAGCAGCAGATTTAGATGCTTGAAGAGAGGCAATACTAAC‑TACTGTGAATATTATGCATAAAAAAGTACTACTAAAAGAAGATAATCGAGGTTTGAAACAACAAATTCACAAATTGCAGCAGCGt                                                                                  >  2:341806/12‑141 (MQ=255)                                                gAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTACGCATGAAAAAGTAATAGTAGAAGAAGAAATACGACGTTTGAAACAACAAATTCACCAATTGCAGCACGGTGAg                                                                               >  2:658034/1‑141 (MQ=255)                                                       ggTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATaaa                                                                        >  2:332147/1‑141 (MQ=255)                                                         tataAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATg                                                                      >  2:391695/1‑141 (MQ=255)                                                         tataAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATAAAAACACGCT‑CGGGGATGATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATg                                                                      <  1:379919/141‑1 (MQ=255)                                                               gcagcaGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCAtt                                                                <  1:367631/141‑1 (MQ=255)                                                               gcagcaGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCAtt                                                                >  2:415062/1‑141 (MQ=255)                                                                  gcagGTTTAGATACTTCAAGAACAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGAt                                                             <  1:424246/141‑1 (MQ=255)                                                                    agGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGAttt                                                           >  1:436986/1‑141 (MQ=255)                                                                    agGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATAGTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAGATTCACAAATTGGAGCAGCGTGACGAATAACTGGTCAttatatt                                                           >  2:514225/1‑136 (MQ=255)                                                                            aTACTTCAAGAAAAGCACTACTAAC‑AACGGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATaa                                                   >  2:472903/1‑141 (MQ=255)                                                                            aTACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATaa                                                   >  1:412471/1‑141 (MQ=255)                                                                            aTACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATaa                                                   >  1:107134/1‑141 (MQ=255)                                                                            aTACTTCAAGAAAAGCAATACTAAC‑TACTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATaa                                                   <  1:257553/141‑1 (MQ=255)                                                                            aTACTGCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATaa                                                   <  1:546299/141‑1 (MQ=255)                                                                             tACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTGTCATAAt                                                  >  1:598472/1‑141 (MQ=255)                                                                             tACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAAt                                                  >  1:365474/1‑141 (MQ=255)                                                                               cTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCa                                                >  2:243268/1‑141 (MQ=255)                                                                                 tCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAGGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCAtt                                              <  1:444432/141‑1 (MQ=255)                                                                                 tCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAGGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCAtt                                              <  1:380413/141‑1 (MQ=255)                                                                                    aGACAAGCAATACTAAC‑TGCTGTCAATATTATGCATGAAAACGTACTAATAGAAGCAGAACATAGACGTTTGAAACACCAAATTCTCCAATTGCGACAGGGTGAGCAATAAATGGCCATTGATTTTATCATAATCAtcttt                                           >  2:187680/1‑140 (MQ=255)                                                                                      aaaaGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGACGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGGCATTGATTTTATCATAATCATTTTCtt                                         >  1:26528/1‑141 (MQ=255)                                                                                       aaaGCAATACTAAC‑TGCTGGGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGTTCATTGATTTTATCATAATCATTTTCttt                                        >  2:376588/1‑141 (MQ=255)                                                                                        aaGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTtg                                       >  2:555746/1‑141 (MQ=255)                                                                                         aGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTtgt                                      >  1:308918/1‑141 (MQ=255)                                                                                           cAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTgtgt                                    >  1:343775/1‑141 (MQ=255)                                                                                           cAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTgtgt                                    >  2:574434/1‑141 (MQ=255)                                                                                              tACTAAC‑TGCTGTGAATATTGTGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTAtt                                 >  1:468234/1‑141 (MQ=255)                                                                                                              atatTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGAt                  <  1:243268/141‑1 (MQ=255)                                                                                                               tattatGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGAtt                 >  2:254061/1‑141 (MQ=255)                                                                                                                   atGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCAGa             >  2:526100/1‑141 (MQ=255)                                                                                                                    tGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCAGAc            <  1:391695/141‑1 (MQ=255)                                                                                                                            aaaGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCAGACGAGGtttt    <  1:503388/141‑1 (MQ=255)                                                                                                                             aaGTACTACTAGAAGAAGAAAATCGACGTTTGCAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTACCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCAGACGAGGttttt   <  1:445675/141‑1 (MQ=255)                                                                                                                              aGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCGTACGCGGTTTTTg  >  2:398802/1‑141 (MQ=255)                                                                                                                              |                                                                                                                                             CAGAGCACATACGATATGTAGCACATTTAGTTGATGATAAAATAAAAGAATTAGGGTATAAAGCAGCAGGTTTAGATACTTCAAGAAAAGCAATACTAAC‑TGCTGTGAATATTATGCATGAAAAAGTACTACTAGAAGAAGAAAATCGACGTTTGAAACAACAAATTCACAAATTGCAGCAGCGTGAGCAATAAATGGTCATTGATTTTATCATAATCATTTTCTTTGTGTATTTTGTCATCGTTGGATTCAGACGAGGTTTTTG  >  CP000730/1150966‑1151230

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