Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,563,642	G→A	R308H (CGC→CAC)	LY180_07445 →	transposase IS609

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,563,642	0	G	A	91.3%	28.2 / ‑0.8	24	R308H (CGC→CAC)	LY180_07445	transposase IS609
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (13/8);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

GACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTA  >  CP006584/1563346‑1563933                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                    gaCGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGcacc                                                                                                                                                                                                                                                                                                   >  1:77261/1‑299 (MQ=17)         cGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTCTCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCCAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGCCAAAGTCCGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGCCAAACTCAGGTTTAAACCGTTCCATCCGGGACCGGGGCGGGAACGAACTGCCCGGCCCGCta                                                                                                                                                                                                                                                                                           >  2:160872/1‑298 (MQ=11)         agaaaCTTCAGCCACAGGTCAGCCGCAAGGTCAAATTCAGCAACAACTGCCAGAAGCAGAAACGCAATATTCAGCGCCTGCATTCCGGTATCGCAAAGATCCGCAGGTGAGACCTGCACAAAGTCACAACGACCGCCAGCAAAAGCCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGGCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCtt                                                                                                                                                                                                                                                                                           <  2:601161/295‑1 (MQ=11)                   acGAAAGTTAAACCGCAAGGTAAAAATCCGCAACAACTGGCAGAAGAAGAAACCCAAAATACAGCGACGGCATTCCTGGATCGAAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCAGTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAagc                                                                                                                                                                                                                                                                                 <  2:41598/298‑1 (MQ=11)                       cAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGct                                                                                                                                                                                                                                                                             <  1:416035/299‑1 (MQ=17)                               ccGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCCGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCCGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCtggcgtg                                                                                                                                                                                                                                                                     >  1:48065/1‑299 (MQ=17)                                           aTTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGtgc                                                                                                                                                                                                                                                         >  1:265976/1‑299 (MQ=14)                                               aGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAAGTTTCAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCCCTGGCGTGGCGGTCAGGtgctgcc                                                                                                                                                                                                                                                     >  2:397522/1‑296 (MQ=9)                                                        tGGCAGAAGCAGAACCGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCAGTCAGCAAAAACCACGCAAAGATAGTCATCGAGGATGTGAAGGGCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTACACCGTTCGATACAGGATCAAGGCCGGTCTGACAGTCACCACACGCTTGAGTACAACAAGCCCTGCCCGTGCGGTCAGGGCCTTCCTGtcacccc                                                                                                                                                                                                                                            >  2:580632/1‑293 (MQ=11)                                                            agaagcagaaACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGACGGGGCGGAATGTCCGGGCAAAATCAGGTTTCAACCGTTCGATACTGGATCAGGGCTGGTCTGCAATGCACCGCCAGCTTGAGTCCAAGCAGCCCTGGCTTGGCGGTCAGGTGCTTGCGTTTCCGCCAGCg                                                                                                                                                                                                                                        >  2:181462/1‑299 (MQ=11)                                                             gaagcagaaACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTAATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCCGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACACGCCGCTTCGGCGTGGCGGTCAGGTGCTTGCTTTTCCCCCCGCGt                                                                                                                                                                                                                                       >  2:24726/1‑299 (MQ=11)                                                               agcagaaACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTat                                                                                                                                                                                                                                     >  1:145149/1‑298 (MQ=14)                                                               agcagaaACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCAAGCCGGGGCGCAATGTCCGGGCAAACTCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTAGGAAAGGCACCGGCAGCTTGAGTACAAGCAGCCCTGGCTTGGGGGTCAGGTTCTTGCTGTTCCGCCacggtgt                                                                                                                                                                                                                                     >  2:282153/1‑293 (MQ=11)                                                                gcagaaACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTata                                                                                                                                                                                                                                    <  1:350800/299‑3 (MQ=14)                                                                 cagaaACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTatac                                                                                                                                                                                                                                   <  1:481731/299‑1 (MQ=14)                                                                  agaaACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTAtaca                                                                                                                                                                                                                                  >  1:210926/1‑299 (MQ=14)                                                                     aaCGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGc                                                                                                                                                                                                                               >  1:570004/1‑299 (MQ=14)                                                                                     aCTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACt                                                                                                                                                                                                               <  1:247441/299‑1 (MQ=14)                                                                                        gCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACtgtg                                                                                                                                                                                                            <  1:282153/299‑1 (MQ=14)                                                                                                                                                                                                                                                       aaaaTCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGa                                             <  1:160872/299‑1 (MQ=9)                                                                                                                                                                                                                                                               gTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGaa                                     >  1:604760/1‑299 (MQ=9)                                                                                                                                                                                                                                                               gTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGCTCTTGCCTGTGGGGAGATGGTGCAGTCAGGCGGCCCGGTGAAGcgtgaa                                     >  2:190749/1‑294 (MQ=11)                                                                                                                                                                                                                                                                      ccGTTTGATACTGGATCAAGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCg                              <  1:407284/299‑1 (MQ=9)                                                                                                                                                                                                                                                                                            tGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTg        >  1:591223/1‑299 (MQ=9)                                                                                                                                                                                                                                                                                                  gAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGCGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCCGGCACCCACCGAAAAGAATTAGGCAAACGCCCGAACgca  >  2:262585/1‑297 (MQ=11)                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                    GACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTA  >  CP006584/1563346‑1563933

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