Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,613,506	A→T	H9L (CAC→CTC)	yjjV →	putative DNase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,613,506	0	A	T	98.5%	220.7 / NA	67	H9L (CAC→CTC)	yjjV	putative DNase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base T (30/35);  minor base C (1/0);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAAC  >  NC_000913/4613410‑4613602                                                                                                 |                                                                                                 gATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTg                                                                                               >  1:342582/1‑100 (MQ=255)         ggTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTg                                                                                       <  2:712808/100‑1 (MQ=255)          gTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGa                                                                                      <  1:1720519/100‑1 (MQ=255)                gAAGTGAGCGACGCTCCGCAGGCGAACGCCCCGACATTTAACAATCGGGAGCTGGCTTGCTTCGCCGTTTTATCGACACTCCCTGCCCTTCTCACTTccc                                                                                >  2:1939791/1‑100 (MQ=255)                 aaGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCccg                                                                               <  2:342580/100‑1 (MQ=255)                   gTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCccgcc                                                                             >  1:101835/1‑100 (MQ=255)                      aGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGtt                                                                          >  1:647234/1‑100 (MQ=255)                      aGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGtt                                                                          >  1:2292890/1‑100 (MQ=255)                         gACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAg                                                                       >  2:667315/1‑100 (MQ=255)                           cGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTg                                                                     <  1:1179029/100‑1 (MQ=255)                           cGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTg                                                                     >  1:256709/1‑100 (MQ=255)                           cGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCACTTTGATTTCCCGCCGTTTAGTg                                                                     <  2:1360681/100‑1 (MQ=255)                               ccGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGa                                                                 >  2:814111/1‑100 (MQ=255)                               ccGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGa                                                                 >  2:1322943/1‑100 (MQ=255)                                cGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGAt                                                                >  1:189122/1‑100 (MQ=255)                                 gCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGa                                                                 <  2:976843/98‑1 (MQ=255)                                  cAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGa                                                              >  2:1073747/1‑100 (MQ=255)                                  cAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGa                                                              <  2:308498/100‑1 (MQ=255)                                   aGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGaa                                                             >  1:1028581/1‑100 (MQ=255)                                    ggCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAag                                                            <  2:1201082/100‑1 (MQ=255)                                      cGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAagag                                                          <  2:518733/100‑1 (MQ=255)                                        tctgACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGc                                                        <  1:375368/97‑1 (MQ=255)                                        aaCGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGc                                                        >  2:807913/1‑100 (MQ=255)                                         aCGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGcc                                                       <  1:2170294/100‑1 (MQ=255)                                           gACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGcaga                                                     >  2:375366/1‑97 (MQ=255)                                             cccGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGcc                                                   >  2:427402/1‑100 (MQ=255)                                               cGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAagagatcggaaga                                                 >  1:1826054/1‑91 (MQ=255)                                               cGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTg                                                 <  1:933824/100‑1 (MQ=255)                                                gACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGc                                                <  1:1628534/100‑1 (MQ=255)                                                    tttAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCa                                               >  2:372803/1‑97 (MQ=255)                                                    tttAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAAcg                                            >  1:217145/1‑100 (MQ=255)                                                    tttAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAAcg                                            >  1:2324592/1‑100 (MQ=255)                                                     ttAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAAcgc                                           >  2:461464/1‑100 (MQ=255)                                                      tAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAAcgcg                                          >  2:281201/1‑100 (MQ=255)                                                      tAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAAcgcg                                          >  2:114509/1‑100 (MQ=255)                                                         cAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCgg                                       >  2:2160764/1‑100 (MQ=255)                                                          aaTGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGc                                      <  1:427404/100‑1 (MQ=255)                                                          aaTGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGc                                      <  1:1466524/100‑1 (MQ=255)                                                           aTGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGca                                     <  1:1715245/100‑1 (MQ=255)                                                            tGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGcac                                    >  1:1788010/1‑100 (MQ=255)                                                                  tCTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCgg                              <  1:812535/100‑1 (MQ=255)                                                                   cTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCggg                             <  1:2094859/100‑1 (MQ=255)                                                                   cTGGCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCggg                             <  1:2095030/100‑1 (MQ=255)                                                                     ggCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCg                           >  2:1494327/1‑100 (MQ=255)                                                                      gCTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGt                          >  1:1906969/1‑100 (MQ=255)                                                                       cTTGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGt                          <  1:799209/99‑1 (MQ=255)                                                                         tGATTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAgg                       >  1:1036116/1‑100 (MQ=255)                                                                           aTTTGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCa                     <  1:1946974/100‑1 (MQ=255)                                                                            tttGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCaa                    >  2:121947/1‑100 (MQ=255)                                                                             ttGCCGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAg                   <  1:1020138/100‑1 (MQ=255)                                                                                 cGTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCa               <  2:256707/100‑1 (MQ=255)                                                                                  gTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCAt              >  1:765522/1‑100 (MQ=255)                                                                                  gTTTTATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCAt              >  1:476782/1‑100 (MQ=255)                                                                                    gatctCGACACCCTCTGCCATTTTGAGTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTg            <  2:2113683/96‑1 (MQ=255)                                                                                     ttATCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGt           <  2:795637/100‑1 (MQ=255)                                                                                       aTCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGt           <  1:1939798/98‑1 (MQ=255)                                                                                        tCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTcc        <  2:2010615/100‑1 (MQ=255)                                                                                        tCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTcc        >  1:540547/1‑100 (MQ=255)                                                                                        tCGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTcc        <  2:647233/100‑1 (MQ=255)                                                                                         cGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTgggatc       >  1:2113690/1‑95 (MQ=255)                                                                                         cGACACCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCg       <  2:1020612/100‑1 (MQ=255)                                                                                           acacCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGc     <  1:114511/100‑1 (MQ=255)                                                                                           acacCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGc     <  1:312387/100‑1 (MQ=255)                                                                                           acacCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGc     <  1:711366/100‑1 (MQ=255)                                                                                           acacCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGc     <  1:807915/100‑1 (MQ=255)                                                                                            cacCCTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGc     <  1:667316/99‑1 (MQ=255)                                                                                              cccTCTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAAc  <  1:7425/100‑1 (MQ=255)                                                                                                 |                                                                                                 GATACGCTGGTTGCCGAAGTGAGCGACGCTCCGCAGGCGAACGACCCGACATTTAACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAAC  >  NC_000913/4613410‑4613602

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