Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	754,926	T→A	13.7%	Q180H (CAA→CAT)	USA300HOU_0705 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	754,926	0	T	A	13.7%	72.8 / 2.8	29	Q180H (CAA→CAT)	USA300HOU_0705	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/11);  new base A (2/2);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TCTTTATTTTGTTGAGAGGCGCTTTTGGCAACATAATCTTTATATTTTTTAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCATCTAATTTATTAGGTGTTAAGCCTGCTTTTTCTTCGTATTTTTGTGATTTTTTCATATCTT  >  CP000730/754801‑755065                                                                                                                              |                                                                                                                                            tCTTTATTTTGTTGAGAGGCGCTTTTGGCAACATAATCTTTATATTTTTTAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTAtt                                                                                                                               >  1:159982/1‑140 (MQ=255)       ttttGTTGAGAGGCGCTTTTGGCAACATAATCTTTATATTTTTTAAATGACTTGATGCGTGCTTTCTCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTATGTAACTTTTTATTCAtttt                                                                                                                         >  1:74559/1‑140 (MQ=255)       ttttGTTGAGAGGCGCTTTTGGCAACATAATCTTTATATTTTTTAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTATGTAACTTTTTATTCAtttt                                                                                                                         >  1:268353/1‑140 (MQ=255)       ttttGTTGAGAGGCGCTTTTGGCAACATAATCTTTATATTTTTTAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGGGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCAtttt                                                                                                                         >  2:51729/1‑140 (MQ=255)            ttgAGAGGCGCTTTTGGCAACATAATCTTTATATTTTTTAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTGTGTGTCGTGTTTCAATATTTTTTTGTGACTGTTTTTTCATTTTCGCGa                                                                                                                    >  1:263125/1‑140 (MQ=255)            tggagaTGCGCTTTTGGCAACATAATCTTTATATTTTTTAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGa                                                                                                                    <  1:1141/138‑1 (MQ=255)                 agGCGCTTTTGGCAACATAATCTTTATATTTTTTAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTTTTCTTTTTTTTTTCGTTTTTCAATATTTTTTTTTAACTTTTTATTCATTTTTGGGATTTCt                                                                                                               >  2:283433/1‑140 (MQ=255)                         ttGGCAACATAATCTTTATATTTTATAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGAtt                                                                                                       >  1:233391/1‑140 (MQ=255)                             cAACATAATCTTTATATTTTTTAAATGACTTGATGCGTGCTTTCTCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTCGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTATCGATTCCTTTGCGAtttttt                                                                                                   >  1:90107/1‑140 (MQ=255)                              aaCATAATCTTTATATTTTTTAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCG‑TTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTc                                                                                                 >  2:122800/1‑140 (MQ=255)                                  tAATCTTTATATTTTTTAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTCTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGc                                                                                              <  1:212216/140‑1 (MQ=255)                                          atatTTTTTAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTAt                                                                                      >  2:155504/1‑140 (MQ=255)                                            atTTTTTAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCCGCTCGTTTATCg                                                                                    >  2:75252/1‑140 (MQ=255)                                                      tGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGACttttttat                                                                          >  1:82096/1‑138 (MQ=255)                                                                 gTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCa                                                               >  1:191976/1‑140 (MQ=255)                                                                  tGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTATGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCAt                                                              <  2:74559/140‑1 (MQ=255)                                                                  tGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTATGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCAt                                                              <  2:268353/140‑1 (MQ=255)                                                                    cTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCATCt                                                            <  1:115026/140‑1 (MQ=255)                                                                     tttATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCATCTa                                                           >  2:272463/1‑140 (MQ=255)                                                                                         cgttttGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCATCTAATTTATTAGGTGTTAAGCCt                                       <  1:64051/138‑1 (MQ=255)                                                                                           ttttGTTTTTTTTTTGGTCGTTTTTTATTTTTTTTTTGTAACTTTTTATTCTTTTTTGCTTTTTCTTTGCGATTTTTTTCATCTAGTTTATCGCCTTTTTTCTCGGTTTTCTCATCTAATTTATTAGGTGTTAAGCCTGc                                     <  2:212792/140‑1 (MQ=255)                                                                                            tttGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAGCTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCATCTAGATTGTTAGGTGTTAAGCCtcgc                                    >  1:32991/1‑137 (MQ=255)                                                                                                        tgtgTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTAAAGGGATTTTATCGACTACTTTCTCATTTTTCTCATCTAATTTATTAGGTGTTAAGCCTGCTTTTTCTTCGTAt                        <  1:75252/140‑1 (MQ=255)                                                                                                         gtgtCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCATCTAATTTATTAGGTGTTAAGCCTGCTTTTTCTTCGTAtt                       >  2:41643/1‑140 (MQ=255)                                                                                                          tgtCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCATCTAATTTATTAGGTGCTTACCCTGCTTTTTCTTCGTAttt                      <  1:155504/140‑1 (MQ=255)                                                                                                           gtCGTTTTTCAATTTTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCATCTAATTTATTAGGTGTTAAGCCTGCTTTTTCTTCGTAtttt                     <  2:133111/140‑1 (MQ=255)                                                                                                            tCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCATCTAATTTATTAGGTGTTAAGCCTGCTTTTTCTTCGTAttttt                    <  2:233391/140‑1 (MQ=255)                                                                                                             cGTTTTTCAATTTTTTTTTGTAACTTTTTATTCATTTTAGCGTTTTCTTTGCGTTTTTTTTCAGCTTGTTTATCTCCTTTTTTCTCAGTTTTCTCATCTAATTTATTAGGTGTTAATCCTGCTTTTTCTTCGTATTTTtg                   <  2:81659/140‑1 (MQ=255)                                                                                                                     aaTTTTTTTTTGGAACTTTTTATTCAGTTTAGCGTTTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCATCTAATTTATTAGGTGTTAAGCCTGCTTTTTCTTCGTATTTTTGTGAttttt           <  1:312262/140‑1 (MQ=255)                                                                                                                              ttGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCATCTAATTCATTAGGTGTTAAGCCTGCTTTTTCTTCGTATTTTTGTGATTTTTTCATATCtt  <  1:214696/140‑1 (MQ=255)                                                                                                                              |                                                                                                                                            TCTTTATTTTGTTGAGAGGCGCTTTTGGCAACATAATCTTTATATTTTTTAAATGACTTGATGCGTGCTTTATCAGCTTCTTGTTGGCGTTTTTGTTCTTCTTTGTGTCGTTTTTCAATATTTTTTTGTAACTTTTTATTCATTTTAGCGATTTCTTTGCGATTTTTTTCAGCTAGTTTATCGCCTTTTTTCTCAGTTTTCTCATCTAATTTATTAGGTGTTAAGCCTGCTTTTTCTTCGTATTTTTGTGATTTTTTCATATCTT  >  CP000730/754801‑755065

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