BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	68,880	T→G	13.9%	H390P (CAC→CCC)	araB ←	L‑ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	68,880	0	T	G	13.9%	206.0 / 2.4	86	H390P (CAC→CCC)	araB	L‑ribulokinase
Reads supporting (aligned to +/- strand): ref base T (25/43); new base G (14/4); total (39/47)
Fisher's exact test for biased strand distribution p-value = 2.86e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGAT‑CGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCC  >  NC_000913/68647‑69117

tGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGAT‑CGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGAAGGGAGGTTCAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTACAGTTCGGAGGTGGGCGGTCGTTAAACCAGTCGAGCA‑CCACCGGCAGGGGATcc                                                                                                                                                                                                                                             >  1:494706/1‑239 (MQ=255)

     tGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGAT‑CGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCCGCGGAGCGGCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGGGGCGGCCGTTAAACCACTCGAGCA‑CCACCGCCAGGGGGGCCAgaga                                                                                                                                                                                                                                        >  2:1000080/1‑239 (MQ=255)

     tGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGAT‑CGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCAACCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTGGAGCA‑CCACCGGCACGGGATCCcgagc                                                                                                                                                                                                                                        >  1:696645/1‑238 (MQ=255)

                   gtggcgccccAGTGCCCTCACGTTTGTAACGCGGATCCCCGTAA‑CGGTAAAGCCACCCATGTTTGGGCGTGGCCCAATGGGGGGGCCAGAAACCAACCCGCCGAGCAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTGCGGTGTGCGGCGGCCGTTAACCCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGccc                                                                                                                                                                                                                          <  2:789064/232‑1 (MQ=255)

                                    tCACGTTATTAACGGCGATCCCCTGAT‑CGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCCGCAATCAAACCGCCGAACAGCAGCGGAGCGCCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCC‑GGGGTTGGTTAGCGATCGGTGGGCGGGGGCCGTTAAACAAGTCCAGCA‑CCACCGGCCGGGGAGCCCAGGACGGGTTTTTGGCCCATGCTCCCGTGAGCGcc                                                                                                                                                                                                         >  2:629110/1‑239 (MQ=255)

                                    tCACGTTATTAACGGCGATCCCCTGAT‑CGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGGTGGCGTTCGGGGTGCGGCGGCCGTTAAACCCG‑CGTGGACCCACCGGCAGGGGGTCCCGAGAGCGAGTTTTGGGCCAAGGTTTGGTccgggcg                                                                                                                                                                                                         >  2:466830/1‑233 (MQ=255)

                                            ttaACGGCGATCCCCGTAG‑CGGTAAAGCACTCCATGATGGCGCGTGCGCCAAAGGCGGTGGCAACAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑g                                                                                                                                                                                                <  1:201684/239‑1 (MQ=255)

                                            tgaccGGCGATCCCCTGAC‑CGGTAAAGCCCTCCATGATTGCCCGTGGCCCAAAGGCGGTGGAAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑g                                                                                                                                                                                                <  1:466830/235‑1 (MQ=255)

                                                 cggctaccccTGGTCCGGTAAAGCCCTCCATGATTGGGGGGCCGCCAAGGGCGGTGGCAGCAATCAACCCGCCGAACAGTAGGGGGGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTT‑GGTCAGCGCCGGAAGCA‑gt‑tgtt                                                                                                                                                                                          <  2:538710/233‑1 (MQ=255)

                                                      tCCCCTGAT‑CGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGGTTTTTGGCCCATGCTTCGGTCAGGGCCGGAAGCC‑GT‑TGTTTctg                                                                                                                                                                                      >  1:353560/1‑238 (MQ=255)

                                                        cccTGAT‑CGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGccc                                                                                                                                                                                                                          <  1:859151/202‑1 (MQ=255)

                                                                 ggTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGGGTTCCAGGGACGGGTTTTTGGCCCATGCTTCGGGCCGCGCCGGAAGCG‑GT‑TGTTTTTGGCTGGCGTTGa                                                                                                                                                                           >  1:28063/1‑239 (MQ=255)

                                                                   tAAAGCCCCCCAGGATTGGGCGGGCGCCAAAGGCGGTGGCAGAAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGttt                                                                                                                                                                                         <  2:62641/223‑1 (MQ=255)

                                                                   tAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGGAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGGGATCCAACGACGGAGTTTTTGCCCATGGTTCGGGCCGCGCGGGAAACA‑GG‑TGGTTctggcggggggggatt                                                                                                                                                                         >  2:615554/1‑231 (MQ=255)

                                                                   tAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAAAAGCCGCGGGGCGTAGGGAGCGAGGTTAAGATCGGTACTCACCCCTTTCA‑GGCGGTGGTTAGCGTTCGGGGGGCGGCGGCCGTTAAACCAGTCGGGCA‑CCCCCGGCCGGGGGCCCGAGGGCGGGTTTTTGGCCCACGCCTCGGCTAGCGCCGGAAGGC‑AT‑GGGTTTTGGCTGGCgggtatt                                                                                                                                                                         >  2:986869/1‑233 (MQ=255)

                                                                      aGCCCTCCGTGATTGCGCGGGCGCCAAAGGCGGTGGACGCAATCAACCCCCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGc                                                                                                                                                                      <  1:615554/239‑1 (MQ=255)

                                                                           tCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTTGTTAGCGTTCGGGGGGCGGCGGCCGTTAAAACAGTGGAGCA‑CCACCGGCCGGGGTTCCAGAGACGGATTTTTGGCCCACGCTTTGGGCGGCGCCGGAAGCC‑GTGTTTTTCT‑GCTGGGGTTGGTTTTCGtttt                                                                                                                                                                 >  2:915841/1‑239 (MQ=255)

                                                                             cATGATTGCGCGTGCGCCAAAGGGGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCa                                                                                                                                                               <  1:629110/239‑1 (MQ=255)

                                                                               tGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGGTTTTTGGCCCATGCTTCGGTCACGGCCGGGAGCC‑GT‑TGTTTCTGGCTGGCGTTGATTTGTGGTTTCCGt                                                                                                                                                             >  1:683525/1‑239 (MQ=255)

                                                                                     cgcgTGCGCCAACGGCGGTGGCAGCAATCCAACCGCCGAACAGCAGCGGAGCGTCGGGAGCGAGGGCCAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGGGCGGCGGCCGTTAAACCCGCCGCGCA‑CCCCCCCAGGGGGATCCCGGGGCGGATGTTTTGCCCATGCCTCGGTCACCGCTGG‑CGCACGT‑TTTTTCTGGCTGGCGTGGGTTTGCGTTTTCATTTCCGGa                                                                                                                                                       >  2:338407/1‑239 (MQ=255)

                                                                                       cgTGCGCCAAAGGCGGTGGCCGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTACGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGGGGGCGGCGGCCCTTAAACCACGCGAGAA‑CCACCGGCAGGGGGTCCAGAGACGGATTTTTGGCCCCACCTTTGGTCATCGCCGGCAACA‑GG‑TGTCGCTGGCTGGGGGTTTTTTTGGTTTTCCGTTCCGGAt                                                                                                                                                      >  2:593281/1‑238 (MQ=255)

                                                                                       cgTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGCAATCACCCCTTTCA‑GGCGGTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAACCCAGTCGAGCA‑CCACCGGCAGGGGATCCAGGGACGGATTTTTGTCCCATGCTTCCGTCAGCGCCCGAACCA‑GT‑TGTTTCTGGCTGGCCGTGATTTTGGCTTTTAGTTCGGGagg                                                                                                                                                     >  2:940549/1‑237 (MQ=255)

                                                                                          tcgcCAAAGGCGGTGGCACCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCCCCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCC‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCt                                                                                                                                                   <  1:819405/237‑1 (MQ=255)

                                                                                          gcgcCAAAGGCGGGGGCAGCAAACAACCCCCCGACCAGCAGGGGGGGGTGGGTAGCGAGGTTAAGATCGGTAAGCACCCCTTTCA‑GGCGTTGGTTAGGGTTCGGTGGGCGGGGGCCGTTAAACCAGTCGAGCC‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCt                                                                                                                                                   <  2:1023049/238‑1 (MQ=255)

                                                                                                      gTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGGTTTTTGGCCCATGGTTCGGTCAGCGCCGGAATCA‑GG‑TGTTTCTGGCTGGGGGTGGTTTGCGGTTTCAGTTCCGGCTGCTgggggggcagggg                                                                                                                                      >  2:604337/1‑232 (MQ=255)

                                                                                                                tccaaCCGCCGAACAGCAGCGGGGCGTCCGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGc                                                                                                                            <  1:585894/236‑1 (MQ=255)

                                                                                                                 cAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGcc                                                                                                                           <  1:807238/239‑1 (MQ=255)

                                                                                                                   aacctgcgAACAGCCGCTGAGCGTCGGTAGCGAGGTTAAGAACGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGGTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAg                                                                                                                         <  1:230396/233‑1 (MQ=255)

                                                                                                                     ccgccgGACAGCAGCGGAGCCTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGcc                                                                                                                       <  1:612242/239‑1 (MQ=255)

                                                                                                                     ccgccgAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTCCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGcc                                                                                                                       <  1:403128/239‑1 (MQ=255)

                                                                                                                         cgAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCACCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGGCTTGATTTGCGTTTTCAGTTCCGGCTGCGGGGCGGCAAGCTGTTCCGGCGCCCACACGAga                                                                                                                   >  2:138390/1‑238 (MQ=255)

                                                                                                                          gAAAAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCCCCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGGGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTa                                                                                                                  <  1:795929/239‑1 (MQ=255)

                                                                                                                                cagcGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGAGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGGGATCCAGAGACGGGTTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAATTCCGGGGGCTGGCGGGCCAGCGGTTCCAGGGGCCGGCCGGGAAcgcgga                                                                                                            >  2:795929/1‑237 (MQ=255)

                                                                                                                                 aggggggCGTCGGTAGCGAGGTTAAGGTCGGTAATCACCCCTTTCA‑GGGGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGAc                                                                                                            <  1:623241/232‑1 (MQ=255)

                                                                                                                                  gcGGGGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGAcc                                                                                                           <  1:915841/238‑1 (MQ=255)

                                                                                                                                    ggggCGTCGTTAGCGGGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGGGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGAcc                                                                                                           <  2:324238/233‑1 (MQ=255)

                                                                                                                                      aGCGTCGGTAGCGAGGTTAAGATCGGTAATCCCCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAAc                                                                                                       <  1:565520/238‑1 (MQ=255)

                                                                                                                                       gcggcgGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGGTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAAcc                                                                                                      <  1:354776/238‑1 (MQ=255)

                                                                                                                                               agggAGGTTAAGGTCGGTAAGCACCCCTTGCA‑GGGGTTGGTTAGCGGTCGGGGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTGCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑Ga                                                                                            <  1:314997/236‑1 (MQ=255)

                                                                                                                                                  gAGGTTAAGATCGGTAATCCCCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑Gata                                                                                          <  1:296053/238‑1 (MQ=255)

                                                                                                                                                     gTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAAc                                                                                                       <  1:193175/223‑1 (MQ=255)

                                                                                                                                                        aaGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTC‑GATATCCCCa                                                                                    >  1:119160/1‑238 (MQ=255)

                                                                                                                                                               gTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGCCGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCCACCGAGTACGCGACCAAACCAGGGGAA‑AATACCACAAAACGcgcg                                                                            >  2:486456/1‑236 (MQ=255)

                                                                                                                                                               gTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGTGGCCAGCCGAGTACGGGGCCAACACAGGCGTG‑GATATCGCCAAACGCCGa                                                                            >  1:30292/1‑239 (MQ=255)

                                                                                                                                                               gTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGGGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACCCCGa                                                                            >  1:183686/1‑239 (MQ=255)

                                                                                                                                                               gTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGa                                                                            >  1:539645/1‑239 (MQ=255)

                                                                                                                                                               gTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGa                                                                            >  1:935078/1‑239 (MQ=255)

                                                                                                                                                               gTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGAA‑GATATCACCAAACGCCGa                                                                            >  1:467572/1‑239 (MQ=255)

                                                                                                                                                                tAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAACACAGGCGCA‑GATATCACCAAACGCCgct                                                                           >  1:251653/1‑237 (MQ=255)

                                                                                                                                                                tAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGGCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGa                                                                            <  1:986869/238‑1 (MQ=255)

                                                                                                                                                                 aacccccccTTTCA‑GGCGTTGGTTAGGGTTCGGTGGGGGGGGGCCGTTAAACCAGTCGGGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCAGGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGAt                                                                           <  2:183686/233‑1 (MQ=255)

                                                                                                                                                                     cccccTTTCA‑GGCGTTGGTTGGGGTTCGGTGTGCGGCGGCCTTTAAACCAGTCGAGCC‑CCACCGGCAGGTGATCCAGAGGCGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGcc                                                                      <  1:348795/238‑1 (MQ=255)

                                                                                                                                                                     cccccTTTCA‑GGCGTTGGTTAGCGTTCGGTGGGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGcc                                                                      <  1:940549/238‑1 (MQ=255)

                                                                                                                                                                     cccccTTTCA‑GGCGGTGGTTAGCGGTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGcc                                                                      <  1:794871/238‑1 (MQ=255)

                                                                                                                                                                     aCCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGcc                                                                      <  1:374965/239‑1 (MQ=255)

                                                                                                                                                                     aCCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTACAGCGGCCAGCCGAGTACTCGACCAAACCAGGCGAA‑GAAATCACCAAACCCCGCTTgg                                                                        >  1:1060182/1‑237 (MQ=255)

                                                                                                                                                                      ccccTTTCAGGGGGTTTGTTGGTGTGCGTTGGGGGGCTCGCGTTAAACCCGTCGAGAC‑ACACGGAGTGGGGATCCGGAACTGGTTTTTTGACCCATGCTGTGGTCAGCGCGTGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTGGCGTTTGCAGTTCCTGATGCTGGGCGGCTAGCTGTGCCAGCGGCCAGCCGAGTACGCGAACAAGACAGGCGTA‑GATATCACCAAACGGCGATTGGcc                                                                      <  2:1030166/239‑1 (MQ=255)

                                                                                                                                                                       cccTTTCA‑GGCTTTGGTTGGCGTTCGGTGGTCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGGCGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGcc                                                                      <  1:467124/237‑1 (MQ=255)

                                                                                                                                                                          tttCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCt                                                                  >  1:758445/1‑238 (MQ=255)

                                                                                                                                                                          tttCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGGGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGCCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGGGGCCAGCTGTTCAAGCGCCCCGCCGGGTAACCGACCAAAACCGCCGAC‑ATTATCCCCAAACCCCGCTTGGCCTgtt                                                                  >  2:932658/1‑236 (MQ=255)

                                                                                                                                                                           ttCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTC‑GATATCACCAAACGCCGATTGGCCTGCtt                                                                 >  1:787038/1‑238 (MQ=255)

                                                                                                                                                                             cA‑GGCGTTGGTTAGCGTTCGGTGTGCGGGGGCCGTTAAACCAGTCGAGCA‑CCTCCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCa                                                              <  1:486456/239‑1 (MQ=255)

                                                                                                                                                                               cggggTTGGTTAGCGTTGGGTGGGGGGGGGCCTTTAAACCAGTCGAGCA‑CCGCCGGCAGGGGAACCAGAGACGGATTTTTGGCCCATGTTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAg                                                             <  2:467572/235‑1 (MQ=255)

                                                                                                                                                                                ggCTTTGGTTAGCGTTCGGTGTGGGGCGGCCGTTAAACCAGTCGGGCA‑CCCCCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGGTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGa                                                            <  2:683525/239‑1 (MQ=255)

                                                                                                                                                                                 gggttgttTAGGGTTGGGGTTGCGGGGGCCGTTAAACCGGGCGAACC‑CCTCCGGCAGGGGATCCGGGGTCGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTGTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGa                                                            <  2:696645/236‑1 (MQ=255)

                                                                                                                                                                                  cGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGAcc                                                          >  1:472850/1‑239 (MQ=255)

                                                                                                                                                                                  cGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCAGGCCTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCACTTCCGGATGCTGGGCGGCCACCTGTTCCACCGGCCAGCAGAGTCCGCGGACCAACCCTGCGAA‑GAAAACACCAAAAGCCGAATGGCCGACTTCAAGAcc                                                          >  2:467124/1‑239 (MQ=255)

                                                                                                                                                                                     tGGTTAGGGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGa                                                        <  1:339551/238‑1 (MQ=255)

                                                                                                                                                                                      ggTTAGGGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCC‑CCCCCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATa                                                      <  1:604337/239‑1 (MQ=255)

                                                                                                                                                                                      ggTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATa                                                      <  1:1000080/239‑1 (MQ=255)

                                                                                                                                                                                      ggTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATAACCAAACGCCGATTGGCCTGCTTCCAGACCGATa                                                      >  1:130153/1‑239 (MQ=255)

                                                                                                                                                                                      gctgtgggTTCGGGGTGGGGGGGCCTTTAAACCCGTCGGGGA‑CCCCCGGCAGGTGAACCAGGGAGGGATTTTTGGCCCATGCGTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTGGCGTTTTCAGTGCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATa                                                      <  2:353560/232‑1 (MQ=255)

                                                                                                                                                                                                 ggTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGGGACGGATTTTTGGCCCATGCTTCGGTCAGCACCCGAAGCA‑GG‑TTTTTCTGGCTGGCGTTGAGTTCGGTTTTCAGTTCCGGGTGCGGGGCGGGAAGCTGGTCCTGCGGCCACCCGGGGACCAGGCCCAA‑CAGGAGGAGGAATTCACAACACACCGATTGGCCGGCTGTCAGACCGGAAAAACCaagaaa                                           >  2:970354/1‑234 (MQ=255)

                                                                                                                                                                                                 ggTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGGA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGca                                            >  1:976464/1‑238 (MQ=255)

                                                                                                                                                                                                     tgCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGCCGCA‑GATATCACCAAACGCCGATTGGCCGGCTTCCAGACCGTTAAATCCAGGcaccac                                        >  1:476107/1‑238 (MQ=255)

                                                                                                                                                                                                      gcggggGCCGTAAAACCAGGCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATTCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTGCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACg                                       <  1:779136/238‑1 (MQ=255)

                                                                                                                                                                                                       ggggggcCCTTAAACCAGTCGGCCA‑CCCCCGGGAGGTGACCCAGAGACGGATTTTTGGCTCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGc                                      <  2:912456/234‑1 (MQ=255)

                                                                                                                                                                                                       cggcggCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGc                                      <  1:1063363/238‑1 (MQ=255)

                                                                                                                                                                                                          cggcCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGCCCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCCCGCTGcc                                  >  1:912456/1‑239 (MQ=255)

                                                                                                                                                                                                                     cccgTCGGGCC‑CCCCCGGCAGGTTATCCAGTGTCGTTTTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGAc                       <  2:28063/236‑1 (MQ=255)

                                                                                                                                                                                                                        gTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAAc                                                                                                       <  1:340742/157‑1 (MQ=255)

                                                                                                                                                                                                                        gTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAAc                                                                                                       >  2:340742/1‑157 (MQ=255)

                                                                                                                                                                                                                         tCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCa                   >  1:241692/1‑239 (MQ=255)

                                                                                                                                                                                                                          cgggCA‑CCCCCGGCGGGGGTGCCAGGGTCGGGTTTTTGGCCCTTGCTTCGGGCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGAGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCaa                  <  2:476107/236‑1 (MQ=255)

                                                                                                                                                                                                                          cgcgCA‑CCGCCGGCAGGTTATCCAGGTGCGGTTTTTTGGCCCCTGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATGTGCGTTTACAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCaa                  <  2:548302/236‑1 (MQ=255)

                                                                                                                                                                                                                                       gCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAAc       >  1:580810/1‑238 (MQ=255)

                                                                                                                                                                                                                                            tGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGcc   <  1:726934/237‑1 (MQ=255)

                                                                                                                                                                                                                                            tGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGACCCACGCTGCCATCAACCTGCCCGCAAAAACCTTTAACTGccc  >  1:404863/1‑238 (MQ=255)

TGACCTGGTTTTTCCGCGCGATGCCGCCCAGTGCCATCACGTTATTAACGGCGATCCCCTGAT‑CGGTAAAGCACTCCATGATTGCGCGTGCGCCAAAGGCGGTGGCAGCAATCAAACCGCCGAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCA‑GGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCA‑CCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCA‑GT‑TGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTA‑GATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCC  >  NC_000913/68647‑69117

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`