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	97.2%	78.9 / NA	36	P28P (CCC→CCT)	LY180_10740	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/25);  minor base C (1/0);  total (11/25)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTTCGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTGGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACC  >  CP006584/2250946‑2251463                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                 gaGTTCGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGttt                                                                                                                                                                                                                             <  1:240727/299‑1 (MQ=21)         gcgtttccgagctaTACCGGTCGGATGGCGGCCAGTGAAAGATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCGTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCg                                                                                                                                                                                                                     <  2:581342/286‑1 (MQ=11)                cATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTg                                                                                                                                                                                                              <  1:397328/299‑1 (MQ=21)                      aTACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCgagatcggaagagc                                                                                                                                                                                                        >  1:581342/1‑289 (MQ=19)                           gTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCt                                                                                                                                                                                                   <  1:91673/299‑1 (MQ=21)                            ttCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTg                                                                                                                                                                                                  <  1:652329/299‑1 (MQ=21)                                 aTGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTcac                                                                                                                                                                                             >  1:70402/1‑299 (MQ=21)                                                 tatCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGAtt                                                                                                                                                                             <  1:594513/299‑1 (MQ=21)                                                     aCTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTcc                                                                                                                                                                         <  1:296651/299‑1 (MQ=21)                                                          cgcaccccccgccAGGGCAGCGGCCGAGGTATTTTTCGGGAACGTCACCGCACCGAGTGCGGCCACAAAAAACCGCCCACCTTACCGGGCAGCGCCACAAGGTCAGGCAGTTCCAGGGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGt                                                                                                                                                                    <  2:212712/287‑1 (MQ=11)                                                                   cATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACAtt                                                                                                                                                           >  1:221665/1‑299 (MQ=21)                                                                         gTCATCGACCGAGAGAGTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCttt                                                                                                                                                     <  1:122041/299‑1 (MQ=21)                                                                          tCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTg                                                                                                                                                    <  1:243530/299‑1 (MQ=21)                                                                           cATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGt                                                                                                                                                   >  1:275267/1‑299 (MQ=21)                                                                           cATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGt                                                                                                                                                   <  1:169707/299‑1 (MQ=21)                                                                              cgcccGAGAGATTTTCAGGGAACGTCACTCCACAGAGTTCGGAGACAAAAACCCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTAACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCa                                                                                                                                                <  1:85902/296‑1 (MQ=255)                                                                               gACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCAt                                                                                                                                               >  1:99013/1‑299 (MQ=21)                                                                                   gagagaTTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGtt                                                                                                                                           >  1:222265/1‑299 (MQ=21)                                                                                        aTTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTtctc                                                                                                                                      <  1:519004/299‑1 (MQ=18)                                                                                                  aaCGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAAt                                                                                                                            >  1:308939/1‑299 (MQ=18)                                                                                                                 atgttgccGACAAAAACCCTACCCCCCTACCGCCCAGAGCCCCAATGGCAGGCAGTTCCAGGGCGGCGACCTCCTGCTCGGTCAGCATCGGGGCCGGCATGGGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTc                                                                                                             <  2:308939/292‑1 (MQ=11)                                                                                                                        cGACAAAAAACCGACAACCTTACCGGCCAGCGCCACAAGGTCAGGCAGCTCCAGTGCGGCGACCTCCTGCTGCGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCGTAATCAGTGCATCGACTGAGGATTTGGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGCTTCAGCAGCGTGACCTGTTAGCTAACCTGCTACCCAAGCCTGCACGGAACTTCCAGGTGAATCACATTTTTTTCGGTCGTGGTTTTCTCCGTTCTGCATCGGGTTTACCCGGTCCGTCAgg                                                                                                      >  2:617109/1‑299 (MQ=11)                                                                                                                          aCAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCt                                                                                                    <  1:467786/299‑1 (MQ=18)                                                                                                                            aaaaaaCCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGCAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCCTTCGGTTGCATCAGCGTGCCCTGTTCTAAACCCTGCTCACCCCGCTTGACCGGATTGTCCAGGGTCACTACATTTTCTTTTTTCATGGCTGTCTCACTTATTAATAGGGGCTAAACGGTCCGCCCGGCTgc                                                                                                  >  2:463491/1‑298 (MQ=11)                                                                                                                                          cgtcacgccctgcgcCAGTATGTCCGGCGGTGCCTGTGCGGCCACCGCCTGGTGGGTCGGCTGCGGTGCGGTCTGGTGCGGCAGCACCTTCAGCGGGGGAGCGGCTTCGGAGTTTGCGACCGCGGACAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATGAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATca                                                                                    <  2:338928/288‑1 (MQ=11)                                                                                                                                                cgCCAGCGCCACAAGGGCAGGCAGTTCCCGGGCGGCGCCCTCCTGCTCGGTCAGCCTAGGGGCCGTCAGGCGCGGCCGCACCTTCATCCGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGcc                                                                              <  2:99013/298‑1 (MQ=11)                                                                                                                                                      cgcCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAgg                                                                                                                                                                      >  2:4993/1‑205 (MQ=21)                                                                                                                                                      cgcCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAgg                                                                                                                                                                      <  1:4993/205‑1 (MQ=21)                                                                                                                                                       gcCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAACCAGGCTGACCGGATGAAAATCACAGGCCGATATTg                                                                       <  1:96029/299‑1 (MQ=18)                                                                                                                                                       cccAAAAGGGCAGGCGTTTCCAGTGCGGCCACCTCCTACTCGGTAAGCCTGGGGGCCGTCATGCGCGGCAGCACCTTCATCAGTGCATCGACTTCGGAGTTGTCGACCGCAGCCAGACTGGCACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTg                                                                       <  2:188137/298‑1 (MQ=11)                                                                                                                                                          acaAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTgcgg                                                                    <  1:48222/299‑1 (MQ=18)                                                                                                                                                                                          ctGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCAc                                    <  1:212195/299‑1 (MQ=17)                                                                                                                                                                                             cTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCtt                                 >  1:334253/1‑299 (MQ=17)                                                                                                                                                                                                       aTCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATg                       <  1:22374/299‑1 (MQ=17)                                                                                                                                                                                                        tCGGGGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGt                      <  1:463491/299‑1 (MQ=14)                                                                                                                                                                                                                            agcaCCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGAcc  <  1:198329/299‑1 (MQ=14)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                 GAGTTCGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTGGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACC  >  CP006584/2250946‑2251463

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