Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,018,810	(T)7→6	coding (109/810 nt)	mgt ←	transglycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,018,804	0	T	.	97.1%	154.1 / ‑2.2	35	coding (115/810 nt)	mgt	transglycosylase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (14/20);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAA  >  CP000730/2018667‑2018940                                                                                                                                          |                                                                                                                                         tGAGACAATACTACCTTTATTTTCAATTGTTGTTTGTTCATCCACATTACCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGCATGGTTAATAATTTTTTTAATTGTATTCGTTTACtttt                                                                                                                                       <  2:250198/141‑1 (MQ=255)  gACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAAGATTTTTAATAGTATTCGTTTACttttt                                                                                                                                      <  1:348566/141‑1 (MQ=255)   acacAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACtttttt                                                                                                                                     >  1:248175/1‑141 (MQ=255)     acaAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACttttttct                                                                                                                                   <  2:281100/141‑3 (MQ=255)                   aTTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTg                                                                                                                    <  2:339312/141‑1 (MQ=255)                   aTTTTCAATTTTTCTTAGTTCATCCACATTATCGCAAGAAAAAAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTg                                                                                                                    <  2:39752/141‑1 (MQ=255)                    ttttCAATTTTTCTTAGTTCATCCACATTATCTCGTGAAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGc                                                                                                                   <  2:196595/141‑1 (MQ=255)                      ttCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCa                                                                                                                 <  1:149150/141‑1 (MQ=255)                        cAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAAc                                                                                                               >  2:308854/1‑141 (MQ=255)                         aaTTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACt                                                                                                              <  1:93650/141‑1 (MQ=255)                                   aGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGataat                                                                                                    <  1:253258/141‑1 (MQ=255)                                    gTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGataata                                                                                                   <  1:334702/141‑1 (MQ=255)                                          ccACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGtatt                                                                                             <  2:298019/141‑1 (MQ=255)                                              aTTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAg                                                                                         >  2:338902/1‑141 (MQ=255)                                               ttATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGt                                                                                        <  1:29002/141‑1 (MQ=255)                                                   cgcgTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAgg                                                                                    <  2:203964/141‑1 (MQ=255)                                                    gcgTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGt                                                                                   >  2:210144/1‑141 (MQ=255)                                                        gtAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCgt                                                                               <  1:214797/141‑1 (MQ=255)                                                         tAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCgtg                                                                              <  2:69398/141‑1 (MQ=255)                                                          aGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGTCGGTTTTCCAACTGGTTGATAATACGTATTATAGTGAGGTTCgtgt                                                                             >  2:275557/1‑141 (MQ=255)                                                           gATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGtt                                                                            >  2:61264/1‑141 (MQ=255)                                                                aaaaTACATGATAACAATAAACAATTCGATAATGATTAGAATTGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCata                                                                       <  2:248175/141‑1 (MQ=255)                                                                      cATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTAATATGCTCa                                                                 >  1:95008/1‑141 (MQ=255)                                                                                aTAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAt                                                       >  1:310611/1‑141 (MQ=255)                                                                                  aaaCAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAttt                                                     >  1:42168/1‑141 (MQ=255)                                                                                  aaaCAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAttt                                                     >  1:374932/1‑141 (MQ=255)                                                                                        tGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTTGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAtttgagttt                                               >  2:204875/1‑141 (MQ=255)                                                                                             tAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTa                                          >  1:257879/1‑141 (MQ=255)                                                                                              aaTGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTAAAAATGTTCAATTGAGTTTGAGTAc                                         >  1:132105/1‑141 (MQ=255)                                                                                                 gATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTa                                      <  1:210144/141‑1 (MQ=255)                                                                                                  aTTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTAt                                     <  1:228421/141‑1 (MQ=255)                                                                                                           gTTAATAATATTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCttt                            >  1:79561/1‑141 (MQ=255)                                                                                                                 aataTTTTTAATAGTATTCGTTTACTTTTTTT‑TTTTTTTTCTGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGc                      >  2:334453/1‑141 (MQ=255)                                                                                                                 aataTTTTTAATAGTATTCGTTTAC‑TTTTTTCTTTTTTTTTTTTTTTTCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGc                      >  2:41460/1‑141 (MQ=255)                                                                                                                      ttttAATTGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTACATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTg                 <  2:282003/141‑1 (MQ=255)                                                                                                                            tAGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCctt           <  2:310611/141‑1 (MQ=255)                                                                                                                             aGTATTCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCcttc          <  1:360662/141‑1 (MQ=255)                                                                                                                                 ttCGTTTAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCtttt      <  1:338902/141‑1 (MQ=255)                                                                                                                                     tttAC‑TTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTaaaa  <  1:308854/141‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         TGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAA  >  CP000730/2018667‑2018940

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