Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	1,397,926	T→C	10.4%	W292R (TGG→CGG) ‡	tolA →	colicin S4 and filamentous phage transport system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,397,926	0	T	C	10.4%	143.2 / 7.8	48	W292R (TGG→CGG) ‡	tolA	colicin S4 and filamentous phage transport system
Reads supporting (aligned to +/- strand):  ref base T (13/30);  new base C (3/2);  total (16/32)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

GCAGCTGAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACC  >  AE015451/1397713‑1398121                                                                                                                                                                                                                      |                                                                                                                                                                                                    gCAGCTGAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCcgcgcg                                                                                                                                                                               <  1:451834/236‑1 (MQ=255)       gAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAg                                                                                                                                                                           <  1:725464/234‑1 (MQ=255)        aaGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAggg                                                                                                                                                                         <  1:673754/235‑1 (MQ=255)         agagGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGAGCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTg                                                                                                                                                                                                    >  1:789357/1‑207 (MQ=255)                          ccGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGGGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGGGGGCTggtc                                                                                                                                                                                                 >  1:240548/1‑191 (MQ=255)                                cAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCAt                                                                                                                                                                      >  2:1089195/1‑214 (MQ=255)                                    tGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCcgcg                                                                                                                                                                                 <  2:749072/199‑1 (MQ=255)                                    tGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATg                                                                                                                                             <  1:848239/235‑1 (MQ=255)                                    ggCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACggtg                                                                                                                                                               <  2:486993/216‑1 (MQ=255)                                     gCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCatgcg                                                                                                                                                                         <  1:152419/206‑5 (MQ=255)                                       ccAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGAAGGCGGACGTGCCCAGCCCTCGGCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTg                                                                                                                                          <  1:407518/235‑1 (MQ=255)                                        cAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGAAGGCGGACGTGCCCAGCCCTCGGCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGt                                                                                                                                            >  2:407517/1‑232 (MQ=255)                                          gaagaagaaGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCg                                                                                                                                       <  1:849978/235‑1 (MQ=255)                                              aagaagGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCcgcgc                                                                                                                                                                                <  2:446118/190‑1 (MQ=255)                                                 aagGCCCAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTcc                                                                                                                                                                                    <  2:106850/183‑1 (MQ=255)                                                      acaGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGTTCCTGCAGATCAACATGTTGCCGGa                                                                                                                                     >  1:378121/1‑225 (MQ=255)                                                        aGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGc                                                                                                                                                        >  2:1101435/1‑204 (MQ=255)                                                        aGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCAc                                                                                                                         >  1:11555/1‑235 (MQ=255)                                                          gAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCatgcgg                                                                                                                                                                        <  1:604029/186‑6 (MQ=255)                                                          gAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCa                                                                                                                       <  1:549520/235‑1 (MQ=255)                                                                 cccGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGAc                                                                                                                                    <  1:814703/215‑1 (MQ=255)                                                                 cccGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCa                                                                                                                 <  1:477611/234‑1 (MQ=255)                                                                    gCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTg                                                                                                            <  1:12098/236‑1 (MQ=255)                                                                     cAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTg                                                                                                            <  1:522963/235‑1 (MQ=255)                                                                            gcAGAAGCCAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGa                                                                                                                                     <  1:864911/203‑1 (MQ=255)                                                                              agaagaCAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCACCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGGGGGCTGGGCACGTCCGCCTTcc                                                                                                                                                                                    >  2:690357/1‑154 (MQ=255)                                                                                         caGCCCAGGCACTGGCCGACGTGTCGTCCGATGCCACCGAGCGGCAGCCGGCGCTGGCCGACGAGCAGGGTGACGAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGcc                                                                                                                                        <  2:125749/186‑1 (MQ=255)                                                                                           gCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGc                                                                                       <  1:1034221/234‑1 (MQ=255)                                                                                                         cGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGAAGGCGGACGTGCCCAGCCCTCGGCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTc                                                                         >  1:994573/1‑234 (MQ=255)                                                                                                          gAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTcggc                                                                      >  1:674307/1‑236 (MQ=255)                                                                                                                             ccGAGCGGCAGCGGGCGCTGGCCGACGAGCAGGGGGACCAGGTGGCCGGCGACTTCGACGACCTGATCAGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAggg                                                                                                                                                                         <  2:53287/118‑1 (MQ=255)                                                                                                                              cGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGAc                                                                                                                                    <  1:247581/154‑1 (MQ=255)                                                                                                                                           gcggTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGAc                                                                                          <  1:150836/183‑1 (MQ=255)                                                                                                                                              cTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTg                                       <  1:1067374/231‑1 (MQ=255)                                                                                                                                               tGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGAAGGCGGACGTGCCCAGCCCTCGGCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTgg                                                                 <  2:994572/204‑1 (MQ=255)                                                                                                                                                  ccGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGc                                                                                                          <  2:65228/160‑1 (MQ=255)                                                                                                                                                                gACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGcc                <  1:681944/236‑1 (MQ=255)                                                                                                                                                                      gTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAgtggtg                                                                                            <  2:1075121/154‑1 (MQ=255)                                                                                                                                                                      gTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGAGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCg           <  1:700217/235‑1 (MQ=255)                                                                                                                                                                              cGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGt                                                                                   >  2:247885/1‑155 (MQ=255)                                                                                                                                                                              cGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACa                                                                             >  2:914803/1‑161 (MQ=255)                                                                                                                                                                              cGACTTCGACGACCTGATCCGCATGCGCGCGGAAGGCGGACGTGCCCAGCCCTCGGCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAAcc  >  1:934950/1‑236 (MQ=255)                                                                                                                                                                                aCTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTg                                                                  >  2:324759/1‑170 (MQ=255)                                                                                                                                                                                  ctCGGCGACCTGCTCCGCAGGCGCGCGGCGGAGGGCTGGGCACGTGCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAg                                                                                                                                                      <  2:157734/83‑1 (MQ=255)                                                                                                                                                                                   tCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGc                                                                                       >  1:269021/1‑146 (MQ=255)                                                                                                                                                                                      gcgacCTGAGCCGCAGGCGGGCGGCCGAGGGCTGGGCACGGCGGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTg                                                                                                            <  2:23208/121‑1 (MQ=255)                                                                                                                                                                                                     tGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTgg                                                                 >  2:26512/1‑150 (MQ=255)                                                                                                                                                                                                       cgcgcgGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGa                                      <  2:892630/175‑1 (MQ=255)                                                                                                                                                                                                                      |                                                                                                                                                                                                    GCAGCTGAAGAGGCCAAGAAAAAGGCCGCTGCAGATGCCCAGAAGAAGAAGGCACAGGAAGCGGCCCGCAAGGCGGCAGAAGACAAGAAAGCCCAGGCACTGGCCGAGCTGTTGTCCGATACCACCGAGCGGCAGCAGGCGCTGGCCGACGAGCAGGGTGACCAGGTGGCCGGCGACTTCGACGACCTGATCCGCATGCGCGCGGCCGAGGGCTGGGCACGTCCGCCTTCCGCGCGCAAGGGCATGACGGTGGTCCTGCAGATCAACATGTTGCCGGACGGTACCATCACCAATGTCAGCGTGGCCCGGTCCAGTGGTGACGGCCCGTATGACAGTTCGGCGGTGGCTGCGGTGAAGAACATTGGTCGTTTGACCGAGATGCAGGGTATGAAGCCGAGCGATTTCAACC  >  AE015451/1397713‑1398121

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