Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,572,548	A→T	11.5%	L351H (CTC→CAC)	yddA ←	putative multidrug transporter subunit of ABC superfamily, membrane component/ATP‑binding component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,572,548	0	A	T	11.5%	115.8 / 2.1	44	L351H (CTC→CAC)	yddA	putative multidrug transporter subunit of ABC superfamily, membrane component/ATP‑binding component
Reads supporting (aligned to +/- strand):  ref base A (18/19);  new base T (3/2);  total (23/21)
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.94e-01

TTTCAGTAATAGCCATTTGCCTGGCGAAACATGAAAGTTCAGGTTCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTGGCGCGATTTCATCAGTCCGCCCAGATTGATCTGCCC  >  CP009273/1572403‑1572690                                                                                                                                                  |                                                                                                                                               tttCAGTAATAGCCATTTGCCTGGCGAAACATGCAAGTTCAGGTTGTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTCCCGCATGTTGGCAATTTTTAGGCTTATTGGTAGGGCGCGGTTCAGTGTGtt                                                                                                                                             >  1:41602/1‑149 (MQ=255) tttCAGTAATAGCCATTTGCCTGGCGAAACATGAAAGTTCAGGTTCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCCGTGAGtt                                                                                                                                             >  2:325804/1‑149 (MQ=255) tttCAGTAATAGCCATTTGCCTGGCGAAACATGAAAGTTCAGGTTCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGtt                                                                                                                                             >  2:563164/1‑149 (MQ=255)  ttCAGTAATAGCCATTTGCCTGGCGAAACATGAAAGTTCAGGTTCTCTAATATAATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGAAGGGCGCTGTCCAGTGCGTTg                                                                                                                                            >  2:228356/1‑149 (MQ=255)    ccgcaaTAGCCATTTGCCTGGCGAAACATGAAAGTTCAGGTTCTCCAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGAt                                                                                                                                          <  2:241138/145‑1 (MQ=255)      gcaaTAGCCATTTGCCTGGCGAAACATGAAAGTTCAGGTTCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATgg                                                                                                                                        >  1:766600/3‑149 (MQ=255)            gCCATTTGCCTGGCGAAAGATGAAAGTTCAGGTTCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGGACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCa                                                                                                                                  <  1:70927/149‑1 (MQ=255)            gCCATTTGCCTGGCGAAACATGAAAGTTCAGGTTCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCa                                                                                                                                  <  2:783786/149‑1 (MQ=255)                 ttGCCTGGCGAAACATGAAAGTTCAGGTTCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACaa                                                                                                                             <  2:713442/149‑1 (MQ=255)                   gCCTGGCGAAACATGAAAGTTCAGGTTCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCAAATTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGAGTGTTGATGGAACGCaatcaatc                                                                                                                           >  1:24985/1‑142 (MQ=255)                   gCCTGGCGAAACATGAAAGTTCAGGTTCTCTAATATGATCTTATTATCAGGCGGACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGTGTTGATGGAACTCATACAAGc                                                                                                                           >  1:614982/1‑149 (MQ=255)                     cTGGCGAAACATGAAAGTTCAGGTTCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGcga                                                                                                                         >  1:519169/1‑149 (MQ=255)                     cTGGCGAAACATGAAAGTTCAGGTTCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGcga                                                                                                                         >  1:519052/1‑149 (MQ=255)                                           gTTCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAg                                                                                                   <  2:833932/149‑1 (MQ=255)                                            ttCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAGGCGATCGATAACCGCAGCCAGTTCAGc                                                                                                  <  2:801855/149‑1 (MQ=255)                                              ctctcgcatGAGCTTAGTAGCGGGCGCACGAATACTCGCATCAGCCACTCGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCaa                                                                                                <  1:426990/142‑1 (MQ=255)                                                 tAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGtt                                                                                             >  2:711982/1‑149 (MQ=255)                                                        aTCTTATTATCAGGCGTACGAATACTCGCATCAGCCAATTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATATCCGCAGCCAGTTCAGCAAGTTCGTCata                                                                                      <  2:568262/149‑1 (MQ=255)                                                          cTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTCAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATAtt                                                                                    >  2:578533/1‑149 (MQ=255)                                                                atCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTataaa                                                                              >  1:790720/1‑149 (MQ=255)                                                                atCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTataaa                                                                              >  1:818257/1‑149 (MQ=255)                                                                atCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTAAGCAAATTCGTCATATGTataaa                                                                              >  1:169777/1‑149 (MQ=255)                                                                 tCAGGCGTACGAATACTAGCATCAGCCACTTGTACCACATGTTGGCAATTTGTAGGCTTGATCGTAGGGCGCTGTTCAGTGGGTTTATGGACATCATAAAAGCCAGCGATGACCGCAGCCAGTTCAGCAAGTTCGCCATATTTataaat                                                                             >  1:2185/1‑149 (MQ=255)                                                                                tccGCATCAGCCACTGGTACCGCATGTTGGGAATTTTTAGGCTTATTCGTAGGGCGCTGTGCAGCGAGTTGATGGAACTAATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCTTATTTATAAATAAACCAGCTTAAAtt                                                              <  2:236553/147‑1 (MQ=255)                                                                                   gCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGtt                                                           >  2:430747/1‑149 (MQ=255)                                                                                         gCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGGTCGGCATATTTATAAATAAACAAGATGAAATCGTTCGCTAc                                                     >  1:461106/1‑149 (MQ=255)                                                                                         gCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGAGGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTAATATTTATAAATAAACCAGCTTAAATTGTTCGATAc                                                     >  2:797382/1‑149 (MQ=255)                                                                                          ccACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATAca                                                    >  1:513599/1‑148 (MQ=255)                                                                                                        tGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTg                                      <  2:634858/149‑1 (MQ=255)                                                                                                        tGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTg                                      <  2:22847/149‑1 (MQ=255)                                                                                                          ttGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGGTGGAACTCATACAAGCGGTCGATAACCGCAGCCAGTTCAGCAAGTTCGTCGTATTTTTAAATAAGGCAGCTTGAATTGTTGGATAGCAgtgtaatgggttggc                                    >  1:125004/1‑135 (MQ=255)                                                                                                               aagttttAGGCTTATTCGTAGGGCGCTGTTCAGTGTGTTGATGGCAAACCGACAATCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTGGCGCGAt                               <  2:24985/146‑1 (MQ=255)                                                                                                               aaTTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTGGCGCGAt                               <  2:802551/149‑1 (MQ=255)                                                                                                               aaTTTTTAGGCTTACTCGTAGGGCGCTGTTCAGTGTATTCATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAACTGCCTGGCGCGAt                               <  2:614982/149‑1 (MQ=255)                                                                                                                    ttAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTGGCGCGATTtcat                          <  2:655117/149‑1 (MQ=255)                                                                                                                             ttCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATTCCTGGCGCGATTTCATCAGTCCGcc                 <  2:519052/149‑1 (MQ=255)                                                                                                                             ttCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATCCCTGGCGCGATTTCATCAGTCCGcc                 <  2:519169/149‑1 (MQ=255)                                                                                                                              tCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTGGCGCGATTGCATCAGTCCGccc                >  2:412406/1‑149 (MQ=255)                                                                                                                                     gcgcCGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTGGCGCGATTTCATCAGTCCGCCCAGATTGa         <  1:409575/149‑1 (MQ=255)                                                                                                                                      cgcTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTGGCGCGATTTCATCAGTCCGCCCAGATTGAt        >  1:656266/1‑149 (MQ=255)                                                                                                                                         tGTTCAGTGAGTTGATAGAACTCATACAAACGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTGGCGCGATTTCATCAGTCCGCCCAGATTGATCTg     >  2:546255/1‑149 (MQ=255)                                                                                                                                          gTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTGGCGCGATTTCATCAGTCCGCCCAGATTGATCTGc    <  2:790720/149‑1 (MQ=255)                                                                                                                                          gTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTGGCGCGATTTCATCAGTCCGCCCAGATTGATCTGc    <  2:818257/149‑1 (MQ=255)                                                                                                                                            tCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTGGCGCGATTTCATCAGTCCGCCCAGATTGATCTGccc  <  1:65146/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                               TTTCAGTAATAGCCATTTGCCTGGCGAAACATGAAAGTTCAGGTTCTCTAATATGATCTTATTATCAGGCGTACGAATACTCGCATCAGCCACTTGTACCGCATGTTGGCAATTTTTAGGCTTATTCGTAGGGCGCTGTTCAGTGAGTTGATGGAACTCATACAAGCGATCGATAACCGCAGCCAGTTCAGCAAGTTCGTCATATTTATAAATAAACCAGCTTAAATTGTTCGATACCAGCATAAATGCCTGGCGCGATTTCATCAGTCCGCCCAGATTGATCTGCCC  >  CP009273/1572403‑1572690

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