Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	2,251,239	G→A	P28P (CCC→CCT)	LY180_10740 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	2,251,239	0	G	A	96.2%	109.6 / ‑5.2	53	P28P (CCC→CCT)	LY180_10740	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (21/29);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GAGTTCGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTGGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTAC  >  CP006584/2250946‑2251473                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           gaGTTCGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGtt                                                                                                                                                                                                                                        <  1:485236/298‑1 (MQ=21)       gtgaaGGCTAATGGGATACGGTTCGGATGGCGGCTAATGAAATATCACTGCACAATCCGCCATCAGGTCATCGACCGAGAGAGTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAAGCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTtcatca                                                                                                                                                                                                                                  <  2:195689/294‑1 (MQ=11)        gTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCCCAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCCCGGTCAGCCTCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCACCGACTCCGGAGTTTGCGACCGCAGCCAGCCTGACCCCCCCCAGCGGCCCGGCACTTGGGTTccaccg                                                                                                                                                                                                                                  >  2:403927/1‑295 (MQ=11)                  tggtgtcccgTTCTGGTGGCGGCATGTGAAATTGCCCTGCCACATCCTCGCGCAGGGCATGGGCCGAGAGGTGTTTCGGGAACGTCACTGCACCGAGTTAGGCGACAAAAAACCGACCACCGTACCGGCCAGCGCCAGAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGt                                                                                                                                                                                                                       <  2:983181/289‑1 (MQ=11)                    ggagtacgtgcTGATGGCGGCCGGTGAACAATCACTGCGCAGTCCGCCATCAGGTGCTCGAACGAGAGATTTGTCGGGAACGTCCCTGCGCCGAGGTCGGCGACAGAAAACCGACCGCCGTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGGCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGGCTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTc                                                                                                                                                                                                                     <  2:959416/288‑1 (MQ=11)                            ttCTGATGGCGGCCATGGAAACTTCACCGCCACCGCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCt                                                                                                                                                                                                             <  2:539724/298‑1 (MQ=255)                                      ggcCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTCCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGCCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAGCCTGCTcaccaag                                                                                                                                                                                                   >  2:950721/1‑296 (MQ=21)                                          aGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCCTCGACTTCGGAGTTTGCGACCGCAGCCAGACTGCCCCCCCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGGGTGACCTCTTCCATACCCCTCTCACCCCGCTTg                                                                                                                                                                                               >  2:364477/1‑298 (MQ=11)                                            gtaaaTATCCCTGCCACATCCGCCAGCAGGTCATCGCCCGAGGGATTTTTCGGGAACGTCGATGCACCGAGTTCGGCGACAAAAAACCGACCCCCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGAc                                                                                                                                                                                             <  2:253532/296‑1 (MQ=255)                                                      cTGCCCCATCCGCCTTCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAGACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTcc                                                                                                                                                                                   <  1:840272/298‑1 (MQ=21)                                                        gCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAg                                                                                                                                                                                 <  1:504647/298‑1 (MQ=21)                                                              cgccacggcggggcATGGACCGAGAGATTGTTCGGGAACGTCACGGCACCGGGTTCGGCGACAACAGACCGCCCACCTTACCGGCCAGCGCAACAAGGTCAGGCAGTTCCAGGGCGGCGACCGCCGGCTCGGTCAGCATCGGTGGCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACGGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAAt                                                                                                                                                                           <  2:330905/285‑1 (MQ=11)                                                                   cATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACAt                                                                                                                                                                      <  1:950721/298‑1 (MQ=21)                                                                   cATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACAt                                                                                                                                                                      >  1:559482/1‑298 (MQ=21)                                                                    aTCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACAtt                                                                                                                                                                     <  1:130800/298‑1 (MQ=21)                                                                     tCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCCGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCCCCACGCTTTACCGGATTGTCCCGGGTAATcaacttt                                                                                                                                                                    >  2:169622/1‑293 (MQ=255)                                                                       gggttagcgcacGAGAGAGTTTGCGGGAACGACCGTCCGCCGAGTTCGGCGAAAAAAACCCGACCCCCTTGCCGGCCGGCGGCACAAGGTCAGGCAGTTGCAGTGCGGCGGCCTCCTGCTCGGTCAGCTTCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTc                                                                                                                                                                  <  2:604033/287‑1 (MQ=11)                                                                                  cGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATgg                                                                                                                                                       <  1:972691/298‑1 (MQ=21)                                                                                  cGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATgg                                                                                                                                                       >  1:100731/1‑298 (MQ=21)                                                                                   gagagaTTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGt                                                                                                                                                      <  1:224284/298‑1 (MQ=21)                                                                                    agagaTTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCCGGGTAATCACATTTTCTTTGTTCATGGtt                                                                                                                                                     >  1:625221/1‑298 (MQ=21)                                                                                              cGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTc                                                                                                                                           <  1:211138/298‑1 (MQ=18)                                                                                                ggAACGTCCCTGCACCGAGTTCGGGGACAAAAAACCCGCCACCTTACCGGCCGGCGACAAGAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTg                                                                                                                                         <  2:384490/298‑1 (MQ=255)                                                                                                 gAACGTCACTGCACCGAGTTCGGCGACAACAAACCGACCACCTTACCGGCCAGCGCCACAAGGCCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTCCATCAGCGTGACCTGTTCGATAACCTGACCACCACCCCTGACCGGATTGCCCAGGGTACTCACATTTTCTTTGTCCTTGGTTTTTTCACTTCTGa                                                                                                                                        >  2:459313/1‑298 (MQ=11)                                                                                                  aacggtcgtGCACCGAGCTGTGCGACGAAAAACCGCACGGCGTCGCGGCCAGCGCGAGAAGGGCCGGCAGTTGCGGTGCGGCGAACTGCTGTGCGGGCAGCGTCGGTGCCGTGAGGCCCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGGCTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCAGCACGCGTGACCGGATTGTCCAGTGTAATCACATTTTCTTTGTTCATGGTTGTCTCACTTCTGaa                                                                                                                                       <  2:896534/290‑1 (MQ=11)                                                                                                        aCTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCgggg                                                                                                                                 >  1:362806/1‑298 (MQ=18)                                                                                                             aCCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGGCCACGGTAATCACATTTTCTTTGTTCTTGGTTCTCTCACTTCTGCATACGGGTTAAc                                                                                                                            >  2:363410/1‑298 (MQ=255)                                                                                                              ccGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAAcc                                                                                                                           >  1:216338/1‑298 (MQ=18)                                                                                                                               aaaCCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGAcc                                                                                                          >  1:387602/1‑298 (MQ=18)                                                                                                                                             aCCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCTCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTCATCTCATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTACCCGGTCAGCCAGGCTGACCGGACGACAATccag                                                                                            >  2:903901/1‑295 (MQ=255)                                                                                                                                                         cacaAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTgc                                                                                <  1:459313/298‑1 (MQ=18)                                                                                                                                                              ggTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGcggcgg                                                                           <  1:143767/298‑1 (MQ=18)                                                                                                                                                               gTCAGGCCGTTCCAGTGCGGCGACCCCCTGCTCGGTCAGCATCGGGGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGgt                                                                          <  2:387602/298‑1 (MQ=18)                                                                                                                                                                  aGGCCGTTGCAGTGCGGCGCCCTCCTGCTCGGGCAGCATCGGGGCCGTCCTGCGCGGCAGCACCGTAATCAGGGCATCGACTTCGGAGTGTGCGACCGCAGCCCGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATGGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGtt                                                                       <  2:362806/298‑1 (MQ=11)                                                                                                                                                                  aGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGtt                                                                       <  1:12292/298‑1 (MQ=18)                                                                                                                                                                  aGGCAGTTCCAGTGCGGCCACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGGTATTGCGGCGGTGtt                                                                       <  1:120894/298‑1 (MQ=18)                                                                                                                                                                   ggCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTg                                                                      >  1:635563/1‑298 (MQ=18)                                                                                                                                                                       gTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTcc                                                                  <  1:879747/298‑1 (MQ=18)                                                                                                                                                                           cAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGcc                                                              <  1:134907/298‑1 (MQ=18)                                                                                                                                                                              tGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGt                                                           <  1:403927/298‑1 (MQ=17)                                                                                                                                                                                   cgACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACg                                                      <  1:169622/298‑1 (MQ=17)                                                                                                                                                                                           tGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTCGGTTTCATCCGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGCCCGGACTGTCCCGGGCAATCCCATTTTCTTTGTTCATGGTTTTCTCCCTTCTGAACCCGGGTTAACCGGTCCGCCAGCCTGACCGGCTGAAACTCCCCGCCCGATATCGCCGCGGTGTTCCTCCCCCCGCTCGCCGCCCTTccc                                              >  2:179451/1‑296 (MQ=11)                                                                                                                                                                                              tCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCtt                                           >  1:271549/1‑298 (MQ=17)                                                                                                                                                                                                      cATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCa                                   <  1:856875/298‑1 (MQ=17)                                                                                                                                                                                                           gTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCCCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTg                              >  2:945336/1‑298 (MQ=14)                                                                                                                                                                                                            tGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGa                             >  1:471815/1‑298 (MQ=14)                                                                                                                                                                                                            tGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGGGTTGCCCCCGCCGGCCGACCCCGTCCACCTTCTCAATCATCTTGa                             >  2:484144/1‑298 (MQ=9)                                                                                                                                                                                                              ccGTCATGCGCGGCAGCACCTTAAGCAGTGCATCGACGTCGGAGTTTGCGACCGCTGCCAGACTGGCACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGGGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAGCCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATg                           <  2:504962/298‑1 (MQ=9)                                                                                                                                                                                                               cGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCCTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATgg                          >  1:335108/1‑298 (MQ=14)                                                                                                                                                                                                                    tGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCGTGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCg                     >  1:787825/1‑298 (MQ=14)                                                                                                                                                                                                                                ccttgagccgctcATCGGCTTGTGAGTTTGCGCCCGCAGCCAGAAGGCCACCGCGCCGCGCCCCGGCA‑GTGGGTTTGCGCGGGGTGGCCTGTTAGTTAACCTGCTCACCACGCTTGACCTGATTGTCCAGGGTAATCACATGTTCTTTGTTCATGGTTTTCTCACGTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGc         <  2:906209/286‑1 (MQ=1)                                                                                                                                                                                                                                    aaTCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCGCATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCt     >  1:543102/1‑298 (MQ=14)                                                                                                                                                                                                                                       cAGTGCATCGACCTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTAGGTTTCATCAGCGTGACCTGCTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTCATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGCTCACCGGTCAGCCAGCCTGACCGGATGAAACTCACCGGCCGATATCGCCGCGGTGTTGCTCCAGCCGGTCGACTCCCTTCACTCTCTCACTCATGTTTATCGTTTCCATTTCGCCCCGCCCCTTAc  >  2:64288/1‑298 (MQ=11)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           GAGTTCGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCA‑TTGGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTAC  >  CP006584/2250946‑2251473

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