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	95.3%	87.8 / NA	43	P28P (CCC→CCT)	LY180_10740	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/23);  minor base C (2/0);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTGGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCG  >  CP006584/2250956‑2251475                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             aGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATCCGCGGCAGCCCCTTCATCAGTGCATCGCCCTCGGAGTTTGCCACCGCACCCAGCCTGACCCCGCGCCGCATCCCCGCCTTCGGTTtcacccacgt                                                                                                                                                                                                                                >  2:335353/1‑291 (MQ=11)              aCAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCCCCGAGTTCGGCGACAAAAAACCGACCACCCTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCCGCAGCACCTTAACCAGTGCATCGCATTCGGACTTTGCGACCGCAGCCAGACTGACACCGCGCACCGTCCCGGCCTTACGTTTCATCACCGCGACCCGTTCGCTa                                                                                                                                                                                                                   >  2:312127/1‑298 (MQ=11)                 gTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAAcc                                                                                                                                                                                                                <  1:476962/298‑1 (MQ=21)                   tCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTg                                                                                                                                                                                                              <  1:642289/298‑1 (MQ=21)                               cAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCtt                                                                                                                                                                                                  <  1:583053/298‑1 (MQ=21)                                 gTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGtt                                                                                                                                                                                                                        >  2:234091/1‑274 (MQ=21)                                 gTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGtt                                                                                                                                                                                                                        <  1:234091/274‑1 (MQ=21)                                 gTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGa                                                                                                                                                                                                <  1:273719/298‑1 (MQ=21)                                                          aTCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACAtt                                                                                                                                                                       <  1:124779/298‑1 (MQ=21)                                                               gTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCtt                                                                                                                                                                  <  1:143997/298‑1 (MQ=21)                                                                tCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCttt                                                                                                                                                                 <  1:658678/298‑1 (MQ=21)                                                                      aCCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCAt                                                                                                                                                           <  1:91492/298‑1 (MQ=21)                                                                              aTTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTtct                                                                                                                                                   <  1:640551/298‑1 (MQ=18)                                                                                    cGGGAACGTCACTGCACCGGGTTCGGCGACAAGAAACCTACCTCCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGTCGACCTCCTGCTCGGTCAGCATCTGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTc                                                                                                                                             <  1:143434/298‑1 (MQ=255)                                                                                     gggAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCCTGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCCCCACGCTTGACCGGCTTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCCTACtttc                                                                                                                                            >  2:594435/1‑296 (MQ=16)                                                                                                     cGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGCGCCGTCATGCGCGGCAGCCCCTTAATCCGTGCATCGCCTGCGGACTTTGCGACCGCAGCCCGACTGACACCGCGCAGCGTCCCGGCATTCGGTTTCATCAGCGTGACCTCTCCGATCCCCTGCCCACCACGCTCTACCTGATTGTCCGGGGCAGTATCATTTTCTTTGTCCATCGTTTCCTCCCCTCTGACTCGGGCttcccgg                                                                                                                            >  2:607458/1‑293 (MQ=11)                                                                                                      ggctcggtcGACAAAAACCCGGCCACCTTACCGCGCAGCGCCACAAGGGCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTAAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGGGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTTCCGGCAGTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCgg                                                                                                                           <  2:219100/290‑1 (MQ=11)                                                                                                      gAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGAcc                                                                                                                                                                                                                            <  1:334096/201‑1 (MQ=21)                                                                                                      gAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGAcc                                                                                                                                                                                                                            >  2:334096/1‑201 (MQ=21)                                                                                                      gAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCgg                                                                                                                           >  1:260840/1‑298 (MQ=18)                                                                                                         ttCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCa                                                                                                                        <  2:106946/298‑1 (MQ=18)                                                                                                         ttCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCa                                                                                                                        <  1:54420/298‑1 (MQ=18)                                                                                                              cGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGAGGTTAACCGGTCAGCCAg                                                                                                                   >  1:202472/1‑298 (MQ=18)                                                                                                                   aaaaaCCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGa                                                                                                              >  1:260860/1‑298 (MQ=18)                                                                                                                          gACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATg                                                                                                       <  1:91564/298‑1 (MQ=18)                                                                                                                                 ttACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATca                                                                                                >  1:395113/1‑298 (MQ=18)                                                                                                                                       gCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTg                                                                                                                                                                                                              <  1:211418/182‑1 (MQ=25)                                                                                                                                       gCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTg                                                                                                                                                                                                              >  2:211418/1‑182 (MQ=25)                                                                                                                                               cacaAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCCAATTTTCTTTGTTCATGGTTTTCTCCCTTCTGAATCGGGGTTACCCGGCAAGCCAGGCTGACCGGGTGAAAATCACAGGCCGATATTga                                                                                  >  2:272149/1‑297 (MQ=255)                                                                                                                                                                    tGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATACTGCGACGGTGTTGCTCCAGCCGGt                                                             >  1:1024/1‑298 (MQ=17)                                                                                                                                                                              tcctGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGtt                                                   <  1:335353/298‑1 (MQ=17)                                                                                                                                                                              tcctGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGtt                                                   >  1:616462/1‑298 (MQ=17)                                                                                                                                                                                          aGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTaa                                                                                                                               <  1:298634/210‑1 (MQ=17)                                                                                                                                                                                          aGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTaa                                                                                                                               >  2:298634/1‑210 (MQ=17)                                                                                                                                                                                           gCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATc                                      <  1:389951/298‑1 (MQ=17)                                                                                                                                                                                              ttggTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATg                                   >  1:334471/3‑298 (MQ=11)                                                                                                                                                                                                 gTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTg                                >  1:586057/1‑298 (MQ=14)                                                                                                                                                                                                         aTGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACCGGCCGATATTGCGGCGGCGTTGCCCCAGCCGGTCGACGCCGTTCACCTTCTCAATCCTGCTGATGGTGTc                        >  2:513499/1‑298 (MQ=9)                                                                                                                                                                                                                  agcaCCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGACTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACGGGCCCATATTGCGGCGCAGTTGATCGCGCCGGCCGCAGCCGTTCACCTCCCCAGACATGTTATTGgttgcgtgttcgcc               >  2:395527/1‑286 (MQ=11)                                                                                                                                                                                                                    aaCCTTAATCCGTGCAACGCGTTCGGAGTTTGCGAACGCAGCCAGACGGCAACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGGGTGACCTGGTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCa             <  2:260840/297‑1 (MQ=11)                                                                                                                                                                                                                      ccTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGc           <  1:31266/298‑1 (MQ=14)                                                                                                                                                                                                                      ccTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGc           <  1:272149/298‑1 (MQ=14)                                                                                                                                                                                                                               gTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCCCAACCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCg  >  1:37483/1‑298 (MQ=11)                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                             AGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTGGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCG  >  CP006584/2250956‑2251475

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