Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP028101	693,938	C→T	pseudogene (890/1244 nt)	C7Y58_03400 ←	IS3 family transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP028101	693,938	0	C	T	97.9%	115.8 / NA	48	pseudogene (890/1244 nt)	C7Y58_03400	IS3 family transposase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/16);  total (32/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTAAACTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGCTTTTCTCCTCTTAAATTAAATTCTGTTACATCTGTAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAATATGATTATCA  >  CP028101/693695‑694028                                                                                                                                                                                                                                                    |                                                                                           aTCTAAACTATTTCCTTTTCTTGACATCCTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTtct                                                                                      >  2:307749/1‑250 (MQ=25) aTCTAAACTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTtct                                                                                      >  1:398343/1‑250 (MQ=25) aTCTAAACTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTtct                                                                                      >  1:13540/1‑250 (MQ=25) aTCTAAACTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTtct                                                                                      >  1:356770/1‑250 (MQ=25)  tCTAAACTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTctcc                                                                                    <  2:100502/251‑1 (MQ=21)        cTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCATGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCttaa                                                                              >  2:531167/1‑251 (MQ=21)            ttCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCttaaatta                                                                          >  2:115719/1‑251 (MQ=21)               cTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTtaaattaaat                                                                       >  1:146038/1‑251 (MQ=21)                 tttCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTaaattaaatt                                                                      >  1:467484/1‑250 (MQ=25)                    cTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTaaattaaatt                                                                      >  2:409778/1‑247 (MQ=25)                    cTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTaaattaaatt                                                                      >  1:11502/1‑247 (MQ=25)                      tGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTa                                                                <  2:628032/251‑1 (MQ=21)                         cATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACAt                                                             <  2:694298/251‑1 (MQ=21)                          aTACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATc                                                            <  2:467484/251‑1 (MQ=21)                           tACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACAt                                                             >  1:366698/1‑249 (MQ=25)                           tACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATc                                                            >  1:312761/1‑250 (MQ=25)                           tACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATc                                                            <  2:312761/250‑1 (MQ=25)                           tACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTATCCTCTTAAATTAAATTCTGTTACAt                                                             >  2:549707/1‑249 (MQ=25)                            aCTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTctcc                                                                                    <  2:272525/225‑1 (MQ=25)                            aCTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTctcc                                                                                    >  1:272525/1‑225 (MQ=25)                            aCTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTc                                                                     >  1:207611/1‑240 (MQ=25)                            aCTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTc                                                                     <  2:207611/240‑1 (MQ=25)                            aCTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCta                                                          >  1:428927/1‑250 (MQ=21)                            aCTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCta                                                          >  1:367803/1‑250 (MQ=21)                            aCTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCta                                                          >  2:572890/1‑250 (MQ=21)                            aCTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCta                                                          >  1:123532/1‑250 (MQ=21)                                 gAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATaaa                                                      >  1:379523/1‑250 (MQ=25)                                  aGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAAc                                                     <  2:330989/250‑1 (MQ=25)                                    ttATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAAcc                                                    <  2:75181/249‑1 (MQ=25)                                      aTCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCAttt                                                >  2:143167/1‑251 (MQ=21)                                        cttcttctCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCAttttt                                              >  2:219413/1‑251 (MQ=21)                                                       tctttcttGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTg                               <  1:115719/251‑1 (MQ=18)                                                          ttcttGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCtt                            >  2:705546/1‑251 (MQ=18)                                                           tcttGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATCTCTTCCATAAGCATCTAATACTGGAGATAAGTATAGCTTTTCTCCTCTTAAATTAGATTc                                                                     >  1:594800/1‑209 (MQ=16)                                                           tcttGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTc                                                                     <  2:594800/209‑1 (MQ=18)                                                               gatctgactactacTGATACTGTCCTCCTTGATCACTATGAAATATCACATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGAATTATCTGCTCTATGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAATATAGTTTTTCTCCTCTTAAATTAACTTCTGTTACATCTATAAACTATTTTTGCTTTGGTGCTGTTGCTTCaaaa                       >  2:268796/14‑251 (MQ=11)                                                               gATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCaaaa                       >  2:529396/1‑251 (MQ=18)                                                                     aataatGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATAATGGAGATAAGTATAGATTTTCTCCTCTTAAAGTAAATTCTGTTACATCTATAGACCATTTTTGATTTGGAGCTGTTGCTTCAAAAtctct                  >  1:532961/1‑250 (MQ=255)                                                                         atGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTaa             >  1:125834/1‑251 (MQ=14)                                                                          tGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAGTTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGc                                    >  1:594419/1‑227 (MQ=9)                                                                          tGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGc                                    <  2:594419/227‑1 (MQ=14)                                                                          tGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAat            >  1:290364/1‑251 (MQ=14)                                                                          tGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAat            >  2:250399/1‑251 (MQ=14)                                                                            cTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAata           <  1:307749/250‑1 (MQ=14)                                                                               aTACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAATATGAt       <  1:250399/251‑1 (MQ=14)                                                                                 aCTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAATATGATTa     <  2:16853/251‑1 (MQ=14)                                                                                   tGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAATATGATTATc   <  2:228608/251‑1 (MQ=14)                                                                                    gTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGTTTTTCTCCTCTTAAATTAAATTCTGTTACATCTATAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAATATGATTATCa  >  1:374484/1‑251 (MQ=11)                                                                                                                                                                                                                                                    |                                                                                           ATCTAAACTATTTCCTTTTCTTGACATACTTTGAGTTATCTTCTTCTCTTTCAGTCTTTCTTGATATGAATAATGCTGATACTGTCATCCTTGATCACTATGAAATATCAAATTTTCATAGTTTTCATTTTCTTTAAATGCTAAATTTAACATATGATTTATCTGCTCCAAGTTAGGACTGCGCGAAATATCATACGAAACTATATATCTTCCATAAGCATCTAATATTGGAGATAAGTATAGCTTTTCTCCTCTTAAATTAAATTCTGTTACATCTGTAAACCATTTTTGATTTGGAGCTGTTGCTTCAAAATCTCTCTTAATATGATTATCA  >  CP028101/693695‑694028

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