Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,964,784	A→C	T151T (ACT→ACG)	LY180_09345 ←	transposase IS609

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,964,784	0	A	C	98.0%	84.6 / ‑4.3	49	T151T (ACT→ACG)	LY180_09345	transposase IS609
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (28/20);  minor base T (1/0);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

GTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAATCTTTAAGCCATTGCGTTTCAGTGGCATTTTTCCACTCAACCAGCCAGGAAG  >  CP006584/1964490‑1965070                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                gTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGa                                                                                                                                                                                                                                                                                             >  1:225186/1‑299 (MQ=17)   aGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGATTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTTACa                                                                                                                                                                                                                                                                                           >  1:272278/1‑299 (MQ=12)          gCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTca                                                                                                                                                                                                                                                                                    >  1:163322/1‑299 (MQ=17)                           gggAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTCACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGAGTGTTTCGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTTATCCGGAGTTGCTACTTCACTTTTTGTCTGAACACTGATGTACCACTTACGCCGGCACTGGCTGACCTTGACATTTTTAACCACAATCGTCAccatc                                                                                                                                                                                                                                                                   >  2:494565/1‑296 (MQ=11)                            tgAATGCAGTTGCTGTTTTTTTCGGTTCGGCTTCTGCACTTTATTGATGGATTTGACCTTGTGGCTTAAGTGGCGCTGAAGTCGCGCCAGCGGCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGTCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGAc                                                                                                                                                                                                                                                                  <  2:598427/298‑1 (MQ=11)                            ggAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGGTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGAATCCAGCCCGCTCATTGATGCTGAAGGGTGCACCGGAGTTGATACTTCACTTTCTGTCTGATTACTGATGTACCACTTACCGCCGGATTGGCTGCCCGTGACATATTTCACAACCACCGTGTCCTGAc                                                                                                                                                                                                                                                                  >  2:439868/1‑299 (MQ=11)                               aTGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGc                                                                                                                                                                                                                                                               >  1:606812/1‑299 (MQ=14)                                  cAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATCCTTCCCTTTCTGTCTGAATATTGATCTACCACTTCCCGCAGGACCGGCTGACCGTCACATTTTTCACAAAACCCGTTACCTGACgccgat                                                                                                                                                                                                                                                            >  2:444166/1‑294 (MQ=11)                                     tCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTcc                                                                                                                                                                                                                                                         <  1:284372/299‑1 (MQ=14)                                           tATTTTGCTTTTCTTCTTCGGCGAGTTATTGCTGAATTTGGCCTTGCGGCTTAACTGTCGCTGAAGTCTCGGCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGAGCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACGTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAgg                                                                                                                                                                                                                                                   <  2:174626/299‑2 (MQ=11)                                                       ctgcttctgcCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCa                                                                                                                                                                                                                                       >  1:189598/1‑299 (MQ=14)                                                          cttctgcCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑Gtt                                                                                                                                                                                                                                   >  1:81237/1‑299 (MQ=14)                                                             ctgcCTGTTATTGCTTAATTTGACCTTGCGGCTTAACTGTCTCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCg                                                                                                                                                                                                                                <  2:81237/299‑1 (MQ=9)                                                                 cAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCCGCCCGATCATTTATGCGGAAGGGTGAACCGGATTTGATACTTCACTTTCTGTCTGCATACTGATGTACCACTTACCGCAGGACTCGCGGACCGTGACATTTTTCACAACCCCCGTGGCCGTCCGGCTAGTCCGGGCCTGCACCCAGCCCA‑CTCTCgggct                                                                                                                                                                                                                            >  2:358555/1‑296 (MQ=11)                                                                     tattgCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGaaaa                                                                                                                                                                                                                        >  1:318858/3‑299 (MQ=14)                                                                              aTTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCCAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGTGTGAACCGGAGTTGATACTTCACTTACTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGCCCGTTACATTTTTACAAACACCCTGGACCTTACGGTAGTTCCGTTAGCACATCCCGCCCA‑GCTTTGCAAGAAAAATActcgtc                                                                                                                                                                                                               >  2:296801/1‑294 (MQ=11)                                                                                     cTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCCCTTACCGCAGGACTGGCTGACCGGGACATTTTTCACCACCCCCTGTACCTGCCGGCTATTCCGGTAACCACTCCCGCCCA‑GTTTCGCAACAAAAATACGGCTGTGTTCCt                                                                                                                                                                                                        >  2:604678/1‑299 (MQ=11)                                                                                                 cTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTa                                                                                                                                                                                            <  1:453976/299‑1 (MQ=14)                                                                                                                 gCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCg                                                                                                                                                                            >  1:458942/1‑299 (MQ=17)                                                                                                                    aGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGaa                                                                                                                                                                         >  1:307263/1‑299 (MQ=17)                                                                                                                                 ttttCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGt                                                                                                                                                            <  1:444166/299‑1 (MQ=17)                                                                                                                                   ttCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTcc                                                                                                                                                          <  1:329916/299‑1 (MQ=17)                                                                                                                                      tGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCTTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCg                                                                                                                                                       >  1:93106/1‑299 (MQ=17)                                                                                                                                        aaaaCTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCt                                                                                                                                                     <  1:328890/299‑1 (MQ=17)                                                                                                                                         aaaCTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCtt                                                                                                                                                    <  1:447686/299‑1 (MQ=17)                                                                                                                                                ttACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAAt                                                                                                                                             >  1:500233/1‑299 (MQ=17)                                                                                                                                                           aaGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAagc                                                                                                                                  >  1:160925/1‑299 (MQ=14)                                                                                                                                                              aCTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCGTTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAagcagc                                                                                                                               >  1:321973/1‑299 (MQ=14)                                                                                                                                                               cTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACTGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGcc                                                                                                                              >  1:138113/1‑299 (MQ=14)                                                                                                                                                                      ccgtgctaGCGTGGCGAGTTTGGCCACGCCAGCTGCCAGGCCGTTCATTGTTGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTcc                                                                                                                       <  2:75852/292‑1 (MQ=11)                                                                                                                                                                             aGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCAgatcgg                                                                                                                >  1:75852/1‑295 (MQ=9)                                                                                                                                                                                        gttAGCCACGCCAGCTTCCATCCCGATCATTGATGCTGAAGGGTGAACCGTAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAg                                                                                                     <  2:189598/298‑1 (MQ=9)                                                                                                                                                                                                cGCCAGCATCCAGACGGAGTATTGATGCTGAAGGGTAAACCGGAGTTGATACTTCACTATCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGCCCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCa                                                                                             <  1:152783/299‑1 (MQ=255)                                                                                                                                                                                                           aGCCCGATCATTGATGCTGAAGGGTGAACCTGATTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGTCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAg                                                                                  <  1:50248/299‑1 (MQ=9)                                                                                                                                                                                                                      tGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCa                                                                       >  1:86324/1‑299 (MQ=14)                                                                                                                                                                                                                        aTGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGACAAATACGGCTGTTTTCCTGATCGAGCTGTACGCCCTGCGGGTAGCGCAATGCATCATTCTGTCCCCGCTTTTTGAATCCGGGAAAACCAGCCCGCTCCCGGAAGAAGTTTTTGTCAGCCCGCTCAAGGTCCTTCAGTTCCTGTTGCaac                                                                     >  2:372218/1‑298 (MQ=11)                                                                                                                                                                                                                                         ggAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGG‑AAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAATCt                                                   <  1:391746/299‑1 (MQ=17)                                                                                                                                                                                                                                              tgATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGcc                                                                                                   <  2:338331/246‑1 (MQ=14)                                                                                                                                                                                                                                              tgATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGcc                                                                                                   >  1:338331/2‑247 (MQ=14)                                                                                                                                                                                                                                               gATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAATCTTTAAg                                              >  1:193171/1‑299 (MQ=17)                                                                                                                                                                                                                                                     tCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCACTCTGCCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCCTCCGGAAGAAGTTTTTGTAAGACCCATCAAGGTCTTTCACTGACTGTTCCCATGGCTGTGGGTGCTAATATATAAGCAATTg                                        >  2:584793/1‑299 (MQ=11)                                                                                                                                                                                                                                                       aCTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGATTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAATCTTTAAGCCATTGCg                                      <  1:439868/299‑1 (MQ=17)                                                                                                                                                                                                                                                             tgcttgagtctggATGGACCCCTTACCTGAGGACTGGCCGACCGTTACACTTTTTACGACCCCCGTTGCCTTACGGCTATGCCCGTAGTG‑ACACAGGCCACGGTTCGGCAGAAAAAGACGGGTGTTGTCCTGATCGAGCTGTCCGCCCGGGGGGGAGCGGAATGCATCATTCGGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTGTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTTTTGCTATGGCTGTGAGGGCGAATCTTTAAGCCATTGCGTTTCAg                                <  2:272278/287‑1 (MQ=2)                                                                                                                                                                                                                                                                    aTACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGcatc                                                                                                                                                                   <  2:31966/161‑1 (MQ=14)                                                                                                                                                                                                                                                                    aTACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGcatc                                                                                                                                                                   >  1:31966/1‑161 (MQ=14)                                                                                                                                                                                                                                                                         gATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAATCTTTAAGCCATTGCGTTTCAGTGGCATTTTTcc                    >  1:435099/1‑299 (MQ=17)                                                                                                                                                                                                                                                                             tACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCATGCGGGTAGCGGAATGCATCATTCTTTCCCCGCTTTTTGACTCGGGGAAAAGAAGGCCGCTGTCGGAAGAATGTTTTGTAAGCCCGCCCCAGCTCTTTCAGTGACCGTTGCAATGGCTGTGAGGGGGAATCTTTAAGCCATTCCGTTTCAGTGGCCTTTTTCCAccc                >  2:541556/1‑297 (MQ=11)                                                                                                                                                                                                                                                                              acaccTTCCCGCAGGACGGGAGCACCGGTACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGGGCTTTACGCCCTGAGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTGTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAATCTTTAAGCCATTGCGTTTCAGTGGCATTTTTCCACTCa               <  2:458942/295‑1 (MQ=11)                                                                                                                                                                                                                                                                                           ggACTGGCTGACCTGTACATTTTTCACACCACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAATCTTTAAGCCATTGCGTTTCAGTGGCATTTTTCCACTCAACCAGCCAGGAAg  <  2:193171/299‑1 (MQ=12)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                GTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA‑GTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTGACTGTTGCAATGGCTGTGAGGGCGAATCTTTAAGCCATTGCGTTTCAGTGGCATTTTTCCACTCAACCAGCCAGGAAG  >  CP006584/1964490‑1965070

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