Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,880,851	T→A	13.6%	intergenic (‑92/+35)	leuE ← / ← dmlR	leucine efflux protein/transcriptional activator of dmlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,880,851	0	T	A	13.6%	127.5 / 12.2	44	intergenic (‑92/+35)	leuE/dmlR	leucine efflux protein/transcriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base T (18/20);  new base A (2/4);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

GCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCA  >  NC_000913/1880710‑1880999                                                                                                                                              |                                                                                                                                                     gCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTg                                                                                                                                               >  2:185501/1‑149 (MQ=255)  ccccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGa                                                                                                                                              <  4:287875/149‑1 (MQ=255)   cccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGat                                                                                                                                             <  2:373144/149‑1 (MQ=255)   cccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGat                                                                                                                                             <  1:108623/149‑1 (MQ=255)    ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGatt                                                                                                                                            <  1:38966/149‑2 (MQ=255)    ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGa                                                                                                                                              >  2:92132/1‑147 (MQ=255)    ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTAAAACCTGa                                                                                                                                              <  1:92132/147‑1 (MQ=255)        ccAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCaa                                                                                                                                        <  4:445008/149‑1 (MQ=255)         cAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAc                                                                                                                                       <  4:342660/149‑1 (MQ=255)         cAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAc                                                                                                                                       <  3:247954/149‑1 (MQ=255)         cAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAc                                                                                                                                       <  2:67755/149‑1 (MQ=255)          aGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAAcc                                                                                                                                      >  1:108369/1‑149 (MQ=255)                    aGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGa                                                                                                                            <  1:281141/149‑1 (MQ=255)                    aGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGa                                                                                                                            <  2:108369/149‑1 (MQ=255)                      tAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGACTATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAAt                                                                                                                          <  2:367094/149‑1 (MQ=255)                      tAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAAt                                                                                                                          <  4:377182/149‑1 (MQ=255)                            aGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTaa                                                                                                                    <  3:123279/149‑1 (MQ=255)                              aaCCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATc                                                                                                                  >  1:266296/1‑149 (MQ=255)                                ccccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAg                                                                                                                <  3:136383/149‑1 (MQ=255)                                 cccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAgg                                                                                                               >  1:116723/1‑149 (MQ=255)                                                   aTTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAAt                                                                                             >  1:131020/1‑149 (MQ=255)                                                       aCCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCgg                                                                                         <  2:221608/149‑1 (MQ=255)                                                          tctTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTa                                                                                      >  3:68108/1‑149 (MQ=255)                                                            tttAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGt                                                                                    >  1:378528/1‑149 (MQ=255)                                                                 ttATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCg                                                                               >  2:199949/1‑149 (MQ=255)                                                                      tttCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCaaa                                                                          >  4:126677/1‑149 (MQ=255)                                                                           aTCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACt                                                                     <  4:158929/149‑1 (MQ=255)                                                                           aTCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACt                                                                     >  1:462430/1‑149 (MQ=255)                                                                                    gATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTaaaa                                                            >  2:11082/1‑149 (MQ=255)                                                                                      ttAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAAt                                                          <  2:462430/149‑1 (MQ=255)                                                                                          aGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTa                                                      <  3:126677/149‑1 (MQ=255)                                                                                          aGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTa                                                      >  3:125311/1‑149 (MQ=255)                                                                                           ggTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTAc                                                     <  1:11082/149‑1 (MQ=255)                                                                                               aaTATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGtt                                                 >  4:274583/1‑149 (MQ=255)                                                                                                     cTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGc                                           >  3:394508/1‑149 (MQ=255)                                                                                                      tAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTc                                                        <  3:16559/135‑1 (MQ=255)                                                                                                      tAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTc                                                        >  4:16559/1‑135 (MQ=255)                                                                                                       aaCCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCAc                                         <  2:378528/149‑1 (MQ=255)                                                                                                        aCCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACt                                        >  2:360025/1‑149 (MQ=255)                                                                                                          cAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACttt                                      >  4:428311/1‑149 (MQ=255)                                                                                                           aaTCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTc                                     <  1:185501/149‑1 (MQ=255)                                                                                                             tCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGc                                   >  4:224850/1‑149 (MQ=255)                                                                                                                           gAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGcc                     <  2:381212/149‑1 (MQ=255)                                                                                                                                       tCAGGTTAAAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCACTGACGTCGCCAGCCTTGAAACATaaa         >  3:454267/1‑149 (MQ=255)                                                                                                                                            ttttAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGcc    <  2:116723/149‑1 (MQ=255)                                                                                                                                              ttAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCa  <  1:360025/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                     GCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGAAACATAAACGGCCCA  >  NC_000913/1880710‑1880999

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