Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,324,546	T→C	E122E (GAA→GAG)	folP ←	dihydropteroate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,324,546	0	T	C	97.3%	120.6 / ‑4.4	37	E122E (GAA→GAG)	folP	dihydropteroate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (18/18);  total (18/19)
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.87e-01

TAGCGATTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGTTTGGATGTATCGACTGAGATCCAGACTTCGAAGCGTTGAGCAATTGCCTCAACCACAGGAATAACA  >  NC_000913/3324421‑3324681                                                                                                                              |                                                                                                                                        tAGCGATTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACTGGTAAACCGGTTTCTGCAGCCGCATCCGGAACGCCAGGCTCGGAAAGGGAGCg                                                                                                                           >  1:271697/1‑140 (MQ=255)  aGCGATTCACTTCTGCATAGACCTCGTCATATTTCCGAGCTTCCTGCACGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCgg                                                                                                                          <  2:176087/140‑1 (MQ=255)  aGCGATTCACTTCTGCAAAGACTTCGTCATACTTCGGTGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCgg                                                                                                                          <  1:283392/140‑1 (MQ=255)    cGATTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGGCAGGCTCGGAAAGGGAGCGGat                                                                                                                        >  1:305879/1‑140 (MQ=255)    cGATTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGat                                                                                                                        <  1:554756/140‑1 (MQ=255)      aTTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGatat                                                                                                                      <  1:132384/140‑1 (MQ=255)      aTTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGatat                                                                                                                      >  1:429294/1‑140 (MQ=255)      aTTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGatat                                                                                                                      <  1:504233/140‑1 (MQ=255)      aTTCACTTCTGCAAACACATCTTCATACTTGGGAGCTTCTTGCGTGGTTTTGGGATTACCCTGCAGATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGCCTCGGAAAGGTAGCGGatat                                                                                                                      >  2:660891/1‑140 (MQ=255)      aTTCACTTCTGAAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCATCCGGCTCCAGAGCGCCTGGCTCGGAAAGGGAGCGGatat                                                                                                                      >  2:283392/1‑140 (MQ=255)        tCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCa                                                                                                                    <  1:22416/140‑1 (MQ=255)        tCACTTCTGCAAAGACATCGTCATACTTCCGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCa                                                                                                                    <  1:220664/140‑1 (MQ=255)                   aaGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAAtgtg                                                                                                         <  2:10546/140‑1 (MQ=255)                                aCTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTc                                                                                            >  1:524922/1‑140 (MQ=255)                                    cGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTg                                                                                        >  2:603299/1‑140 (MQ=255)                                    cGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTg                                                                                        >  2:663179/1‑140 (MQ=255)                                    cGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTg                                                                                        >  2:393434/1‑140 (MQ=255)                                      gAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGAc                                                                                      >  1:577640/1‑140 (MQ=255)                                       aGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACt                                                                                     >  1:144962/1‑140 (MQ=255)                                       aGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACt                                                                                     >  1:497510/1‑140 (MQ=255)                                       aGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACt                                                                                     <  1:570462/140‑1 (MQ=255)                                       aGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCACAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACt                                                                                     <  2:386793/140‑1 (MQ=255)                                         cTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCa                                                                                   >  2:600265/1‑140 (MQ=255)                                          ttCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCAc                                                                                  >  2:118704/1‑140 (MQ=255)                                            ccTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACgg                                                                                >  1:599163/1‑140 (MQ=255)                                                   ggTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACtt                                                                         <  1:358938/140‑1 (MQ=255)                                                          ggATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGttt                                                                  >  1:527099/1‑140 (MQ=255)                                                                  cTGCATATGCATCAGACCAACTTGTAAACCGTTTTGTGCAGCCGCCTCGAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGTTTGGATGTAt                                                          <  2:271697/140‑1 (MQ=255)                                                                  cTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGTTTGGATGTAt                                                          <  2:497510/140‑1 (MQ=255)                                                                       tatGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGTTTGGATGTATCGACt                                                     <  1:66561/140‑1 (MQ=255)                                                                          gCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCTACTTTCGCTGACACACAGATGACTTCTGGTTTGGATGTATCGACTgag                                                  >  2:207272/1‑140 (MQ=255)                                                                          gCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGTTTGGATGTATCGACTgag                                                  >  2:505392/1‑140 (MQ=255)                                                                                      ccggtaaaccggtTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGTTTGGATGTATCGACTGAGATCCAGACTTCg                                      <  2:577640/140‑1 (MQ=255)                                                                                        ggtaaaccggtTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGTTTGGATGTATCGACTGAGATCCAGACTTCGaa                                    <  2:144962/140‑1 (MQ=255)                                                                                          taaaccggtTTCTGCAGCCGCCTCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGTTTGGATGTATCGACTGAGATCCAGACTTCGAAGc                                  >  2:330116/1‑140 (MQ=255)                                                                                                                tCCAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGTTTGGATGTATCGACTGAGATCCAGACTTCGAAGCGTTGAGCAATTGCCTCAACcac            <  1:505392/140‑1 (MQ=255)                                                                                                                  cAGAGCGCCAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGTTTGGATGTATCGACTGAGATCCAGACTTCGAAGCGTTGAGCAATTGCCTCAACCACAg          <  1:393434/140‑1 (MQ=255)                                                                                                                          cAGGCTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGTTTGGATGTATCGACTGAGATCCAGACTTCGAAGCGTTGAGCAATTGCCTCAACCACAGGAATAaca  <  1:660891/140‑1 (MQ=255)                                                                                                                              |                                                                                                                                        TAGCGATTCACTTCTGCAAAGACATCGTCATACTTCGGAGCTTCCTGCATGGTTTTTGGATTTCCCTGCATATGCATCAGACAAACCGGTAAACCGGTTTCTGCAGCCGCCTCCAGAGCGCCAGGTTCGGAAAGGGAGCGGATATCATTAATAATGTGAGCGCCAACTTTCGCTGACTCACGGATGACTTCTGGTTTGGATGTATCGACTGAGATCCAGACTTCGAAGCGTTGAGCAATTGCCTCAACCACAGGAATAACA  >  NC_000913/3324421‑3324681

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑