Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,403,678	A→T	16.0%	V185D (GTT→GAT)	USA300HOU_2267 ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,403,678	0	A	T	16.0%	70.0 / 5.2	25	V185D (GTT→GAT)	USA300HOU_2267	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (10/11);  new base T (2/2);  total (12/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 = 1.00e+00

AATGTTTTAATGGTACTTTTACATCATTTAATATTAAATCATGACTTTCTGTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGTATATGTCTTCACACCATTTAAAATATACCCATCATCAGCTTTAACAGCATGTGTACTTGGTCTTCCCCCTCT  >  CP000730/2403550‑2403814                                                                                                                                 |                                                                                                                                         aaTGTTTTAATGGTACTTTTACATCATTTAATATTAAATCATGACTTTCTGTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAAtt                                                                                                                              <  2:176427/141‑1 (MQ=255)      tttAATGGTACTTTTACATCATTTAATATTAAATCATGACTTTCTGTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACTAACACTTTCTAATtcttct                                                                                                                         >  1:34071/1‑141 (MQ=255)                 ttttACATCATTTAATATTAAATCATGACTTTCTGTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAagcagc                                                                                                              >  2:187108/1‑141 (MQ=255)                    tACATCATTTAATATTAAATCATGACTTTCTGTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCATCACTTTCTAATTCTTCTATATAAGCAGcaac                                                                                                           >  2:146472/1‑141 (MQ=255)                    tACATCATTTAATATTAAATCATGACTTTCTGTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCATCACTTTCTAATTCTTCTATATAAGCAGcaac                                                                                                           >  2:162333/1‑141 (MQ=255)                           tttAATATTAAATCATGACTTTCTGTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAataata                                                                                                    <  1:246929/141‑1 (MQ=255)                            ttAATATTAAATCATGACTTTCTGTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAatat                                                                                                   >  2:151918/1‑141 (MQ=255)                            ttAATATTAAATCATGACTTTCTGTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAatat                                                                                                   >  2:139675/1‑141 (MQ=255)                                        cATGACTTTCTGTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCtt                                                                                       <  1:239312/141‑1 (MQ=255)                                                   gTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGAt                                                                            >  1:43935/1‑141 (MQ=255)                                                   gTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGAt                                                                            >  1:124250/1‑141 (MQ=255)                                                   gTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGAt                                                                            >  1:235293/1‑141 (MQ=255)                                                   gTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGAt                                                                            >  1:23485/1‑141 (MQ=255)                                                      gctcgcATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGta                                                                         >  2:196178/1‑141 (MQ=255)                                                              aCCTAATACTTCCCAATTATCAGCAATTTCAACACCAGCCAAATTCCTGTCTACTAAAAAAAAACCAACACTTTCTAATTCTTCTATATAATCAGCAACAATAATATGTGTTAGTTCTTTACTCATTGATGTATATGTCtt                                                                 <  2:75076/141‑1 (MQ=255)                                                              aCCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGTATATGTCtt                                                                 <  2:133266/141‑1 (MQ=255)                                                                                aTCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGTATATGTCTTCACACCATTTAAAAtata                                               <  1:104549/141‑1 (MQ=255)                                                                                  cGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGTATATGTCTTCACACCATTTAAAATATAcc                                             <  1:219402/141‑1 (MQ=255)                                                                                         ttCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCATCACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGTATATGTCTTCACACCATTTAAAATATACCCatcatc                                      <  1:146472/141‑1 (MQ=255)                                                                                         ttCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCATCACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGTATATGTCTTCACACCATTTAAAATATACCCatcatc                                      <  1:162333/141‑1 (MQ=255)                                                                                              cacCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGTATATGTCTTCACACCATTTAAAATATACCCATCATCAGCtt                                 <  1:151918/141‑1 (MQ=255)                                                                                              cacCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGTATATGTCTTCACACCATTTAAAATATACCCATCATCAGCtt                                 <  1:175330/141‑1 (MQ=255)                                                                                              cacCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGTATATGTCTTCACACCATTTAAAATATACCCATCATCAGCtt                                 <  1:139675/141‑1 (MQ=255)                                                                                              cacCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGTATATGTCTTCACACCATTTAAAATATACCCATCATCAGCtt                                 <  1:106208/141‑1 (MQ=255)                                                                                                                             aCCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGTATATGTCTTCACACCATTTAAAATATACCCATCATCAGCTTTAACAGCATGTGTACCTGGTCTTCCCCctct  >  1:220293/1‑141 (MQ=255)                                                                                                                                 |                                                                                                                                         AATGTTTTAATGGTACTTTTACATCATTTAATATTAAATCATGACTTTCTGTCGCTCGCATACCTAATACGTCCCAATTATCGGCAATTTCAACACCAGGCAAATTCCTGTCTACTAAGAAAAAACCAACACTTTCTAATTCTTCTATATAAGCAGCAACAATAATATGTGTTAGTGCTTTACTCATTGATGTATATGTCTTCACACCATTTAAAATATACCCATCATCAGCTTTAACAGCATGTGTACTTGGTCTTCCCCCTCT  >  CP000730/2403550‑2403814

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