Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,880,549	T→C	T71A (ACC→GCC)	leuE ←	leucine efflux protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,880,549	0	T	C	98.4%	224.7 / NA	64	T71A (ACC→GCC)	leuE	leucine efflux protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (23/40);  total (23/41)
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

TACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATAC  >  NC_000913/1880456‑1880644                                                                                              |                                                                                                tACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCttaat                                                                                           <  2:1077210/100‑1 (MQ=255) tACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaat                                                                                           >  2:1545468/1‑100 (MQ=255) tACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaat                                                                                           <  1:602073/100‑1 (MQ=255)   cccTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaatta                                                                                         <  1:42298/100‑1 (MQ=255)    ccTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaattaa                                                                                        <  2:1544915/100‑1 (MQ=255)    ccTTCAGGGTCGCGTAAAGAATTTTACGCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCttaattaa                                                                                        <  2:936931/100‑1 (MQ=255)      ttCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATg                                                                                      <  1:381520/100‑1 (MQ=255)       tCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGt                                                                                     >  2:237115/1‑100 (MQ=255)       tCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGt                                                                                     <  1:547982/100‑1 (MQ=255)        cAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATg                                                                                      <  1:1556223/98‑1 (MQ=255)        cAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTc                                                                                    <  2:1755091/100‑1 (MQ=255)         aGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCg                                                                                   <  2:957854/100‑1 (MQ=255)           ggTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGcc                                                                                 <  2:827666/100‑1 (MQ=255)           ggTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGcc                                                                                 <  2:429136/100‑1 (MQ=255)           ggTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGc                                                                                  <  1:1510166/99‑1 (MQ=255)           ggTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGc                                                                                  <  1:375099/99‑1 (MQ=255)            gTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCa                                                                                >  2:680249/1‑100 (MQ=255)            gTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCa                                                                                <  2:819873/100‑1 (MQ=255)            gTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCAACAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGcc                                                                                 <  2:1245819/99‑1 (MQ=255)              cgcgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACt                                                                              >  2:1918017/1‑100 (MQ=255)               gcgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTc                                                                             >  2:699419/1‑100 (MQ=255)               gcgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTc                                                                             <  1:978018/100‑1 (MQ=255)                cgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTcc                                                                            >  2:608033/1‑100 (MQ=255)                cgTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTcc                                                                            >  2:1397560/1‑100 (MQ=255)                    aaGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGcc                                                                        <  1:80655/100‑1 (MQ=255)                    aaGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGcc                                                                        <  2:754599/100‑1 (MQ=255)                        nttttACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATg                                                                    >  1:89088/2‑100 (MQ=255)                            tACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAg                                                                <  2:716405/100‑1 (MQ=255)                                   cAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATc                                                         <  1:680248/100‑1 (MQ=255)                                        agagCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATa                                                     <  2:75747/99‑1 (MQ=255)                                         gagCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATAcc                                                   >  1:819006/1‑100 (MQ=255)                                          agCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCg                                                  <  1:528087/100‑1 (MQ=255)                                             aaaTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCAt                                               <  2:353457/100‑1 (MQ=255)                                             aaaTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCAt                                               <  1:791665/100‑1 (MQ=255)                                             aaaTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCAt                                               >  1:785249/1‑100 (MQ=255)                                              aaTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATc                                              <  1:566049/100‑1 (MQ=255)                                              aaTAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATc                                              <  1:699418/100‑1 (MQ=255)                                                 aaaaCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGcc                                           <  2:734342/100‑1 (MQ=255)                                                 aaaaCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGcc                                           <  2:1773903/100‑1 (MQ=255)                                                  aaaCGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATAACGCATCGCCa                                          >  1:909387/1‑100 (MQ=255)                                                     cGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaatn                                       <  2:1571739/100‑2 (MQ=255)                                                      gCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaataa                                      >  1:1561868/1‑100 (MQ=255)                                                       cACCAAGAT‑ACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaataaat                                    <  2:1379159/100‑1 (MQ=255)                                                        aCCAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaataaat                                    <  2:45522/100‑1 (MQ=255)                                                         ccAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCaataaa                                     >  1:1444688/1‑98 (MQ=255)                                                          cAAGATAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATac                                  <  2:988708/100‑1 (MQ=255)                                                               tAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGc                             >  2:204533/1‑100 (MQ=255)                                                               tAACGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGc                             >  2:1852346/1‑100 (MQ=255)                                                                 aCGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAg                           <  2:243991/100‑1 (MQ=255)                                                                  cGTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAgg                          >  2:1612963/1‑100 (MQ=255)                                                                   gTACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGc                         >  2:1757411/1‑100 (MQ=255)                                                                    tACAATGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGcc                        <  1:264028/100‑1 (MQ=255)                                                                        aTGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCaa                    <  2:1424701/100‑1 (MQ=255)                                                                         tGTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAg                   <  2:1561879/100‑1 (MQ=255)                                                                          gTTGAATAATATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGa                  >  1:1799278/1‑100 (MQ=255)                                                                              aataatATCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAAc              <  1:620919/100‑1 (MQ=255)                                                                                taataTCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAAcc             <  1:706403/99‑1 (MQ=255)                                                                                taataTCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCg            >  2:97943/1‑100 (MQ=255)                                                                                  atatCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGcc          >  2:332652/1‑100 (MQ=255)                                                                                  atatCGGGGTGGCCTTAATTAATGCCGCCACTCCAGCCCATGCCAGAAACATCACTACCACCTCGCCAATAAATACACCGCCGGCCGCACGCTAaccacc          >  1:865603/1‑100 (MQ=255)                                                                                    atCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCtt        >  1:1623474/1‑100 (MQ=255)                                                                                     tCGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCttt       <  2:902081/100‑1 (MQ=255)                                                                                      cGGGGTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTc      <  1:1918004/100‑1 (MQ=255)                                                                                       ggggTGGCCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCa     >  1:1735167/1‑100 (MQ=255)                                                                                         ggtggcCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATa   <  2:1202796/100‑1 (MQ=255)                                                                                          gtggcCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATAc  <  1:1545458/100‑1 (MQ=255)                                                                                              |                                                                                                TACCCTTCAGGGTCGCGTAAAGAATTTTACTCCCCAGATAGAGCAAATAAAACGCACCAAGATAACGTACAATGTTGAATAATATCGGGGTGGTCTTAATTAATGTCGCCACTCCAGCCCATGCCAGAAACATCAATACCGCATCGCCAATAAATACACCGCAGGCCGCAAGATAACCGCCTTTCATAC  >  NC_000913/1880456‑1880644

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