BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,618,536	C→A	100%	R267R (CGG→AGG)	yhhI →	putative transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,618,536	0	C	A	98.8%	692.0 / 0.5	249	R267R (CGG→AGG)	yhhI	putative transposase
Reads supporting (aligned to +/- strand): ref base C (2/1); new base A (127/118); total (130/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 249 total aligned reads displayed.

GAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACG  >  CP009273/3618392‑3618680

gagCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGGAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGt                                                                                                                                               >  1:1954635/1‑148 (MQ=18)

gagCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGt                                                                                                                                               >  2:1229476/1‑148 (MQ=18)

gagCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGt                                                                                                                                               >  2:1654581/1‑148 (MQ=18)

gagCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGt                                                                                                                                               <  2:2449118/148‑1 (MQ=18)

gagCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGt                                                                                                                                               >  1:1278376/1‑148 (MQ=18)

   cATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCa                                                                                                                                            >  1:6746/1‑148 (MQ=17)

     tGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCata                                                                                                                                          >  2:3569542/1‑148 (MQ=17)

        cAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCataata                                                                                                                                       >  2:2122410/1‑148 (MQ=17)

         aGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAg                                                                                                                                      <  2:1382598/148‑1 (MQ=17)

         aGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAg                                                                                                                                      >  2:1191913/1‑148 (MQ=17)

          gTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGc                                                                                                                                     <  2:76563/148‑1 (MQ=17)

          gTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGc                                                                                                                                     >  1:2685998/1‑148 (MQ=17)

             aCGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGa                                                                                                                                  >  2:2677519/1‑148 (MQ=17)

                 aaTCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaa                                                                                                                              <  2:2378310/148‑1 (MQ=17)

                  aTCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaaa                                                                                                                             <  2:1800970/148‑1 (MQ=17)

                  aTCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTACCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaaa                                                                                                                             <  1:3433240/148‑1 (MQ=17)

                   tCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACaaaa                                                                                                                            >  2:2513817/1‑148 (MQ=14)

                                 cACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGaaa                                                                                                              <  1:1612139/148‑1 (MQ=21)

                                 cACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGaaa                                                                                                              <  1:2093136/148‑1 (MQ=21)

                                 cACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGaaa                                                                                                              <  2:638788/148‑1 (MQ=21)

                                   cGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATg                                                                                                            <  2:1215259/148‑1 (MQ=21)

                                   cGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATg                                                                                                            <  1:927549/148‑1 (MQ=21)

                                   cGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATg                                                                                                            >  1:319791/1‑148 (MQ=21)

                                   cGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATg                                                                                                            <  2:1257923/148‑1 (MQ=21)

                                   cGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATg                                                                                                            <  1:2477238/148‑1 (MQ=21)

                                   cGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATg                                                                                                            <  2:3427800/148‑1 (MQ=21)

                                     gCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGAc                                                                                                          <  1:1654581/148‑1 (MQ=21)

                                       agagAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACgg                                                                                                        <  2:3249489/148‑1 (MQ=21)

                                        gagaAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGt                                                                                                       >  1:2868551/1‑148 (MQ=21)

                                        gagaAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGt                                                                                                       >  1:1874435/1‑148 (MQ=21)

                                        gagaAGAAATCAGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACTTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTACATAATAGCAGAACAAAAGAAAGAGCCTGAAATGTCGGt                                                                                                       <  2:2205991/148‑1 (MQ=11)

                                             gaaATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGat                                                                                                  >  2:2564319/1‑148 (MQ=21)

                                               aaTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGatat                                                                                                <  1:244487/148‑1 (MQ=21)

                                                aTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                                               <  1:2764379/148‑1 (MQ=21)

                                                aTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                                               <  1:728483/148‑1 (MQ=21)

                                                aTCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtatt                                                                                               <  2:3691722/148‑1 (MQ=21)

                                                  ccGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGAtattat                                                                                             <  2:450630/148‑1 (MQ=21)

                                                      cTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCa                                                                                         <  1:1447041/148‑1 (MQ=21)

                                                      cTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCa                                                                                         >  1:1111156/1‑148 (MQ=21)

                                                        tCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGt                                                                                       >  2:568590/1‑148 (MQ=21)

                                                           tatTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTct                                                                                    >  2:985039/1‑148 (MQ=21)

                                                           tatTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTct                                                                                    >  1:1194329/1‑148 (MQ=21)

                                                            atTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctg                                                                                   <  1:2362503/148‑1 (MQ=21)

                                                               gTTTGCGATGTCCCTGATGCACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTctgctg                                                                                <  2:2655172/148‑1 (MQ=21)

                                                                tttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGa                                                                               <  2:2853271/148‑1 (MQ=21)

                                                                tttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGa                                                                               <  2:1452436/148‑1 (MQ=21)

                                                                tttGAGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGa                                                                               <  1:3461450/148‑1 (MQ=21)

                                                                tttGAGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGa                                                                               <  1:3461319/148‑1 (MQ=21)

                                                                 ttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAt                                                                              <  2:2974860/148‑1 (MQ=21)

                                                                 ttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAt                                                                              <  2:842494/148‑1 (MQ=21)

                                                                 ttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAt                                                                              <  2:506563/148‑1 (MQ=21)

                                                                  tGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAtt                                                                             <  2:3072778/148‑1 (MQ=21)

                                                                    cGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTa                                                                           <  2:2077582/148‑1 (MQ=21)

                                                                     gATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTaa                                                                          <  1:506013/148‑1 (MQ=34)

                                                                     gATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTaa                                                                          <  2:942763/148‑1 (MQ=21)

                                                                     gATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTaa                                                                          <  1:558211/148‑1 (MQ=21)

                                                                      aTGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAAc                                                                         >  1:2593443/1‑148 (MQ=21)

                                                                         tCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGc                                                                      <  2:3149524/148‑1 (MQ=21)

                                                                         tCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGc                                                                      >  1:430851/1‑148 (MQ=21)

                                                                           ccTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTg                                                                    >  1:1361681/1‑148 (MQ=25)

                                                                           ccTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTg                                                                    >  1:3265307/1‑148 (MQ=25)

                                                                            cTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTga                                                                   <  2:1973422/148‑1 (MQ=37)

                                                                              gatgaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTgaga                                                                 >  2:1868459/1‑148 (MQ=37)

                                                                              gatgaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTgaga                                                                 >  1:793696/1‑148 (MQ=37)

                                                                               atgaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGaa                                                                <  2:2888165/148‑1 (MQ=37)

                                                                               atgaACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGaa                                                                >  2:1440657/1‑148 (MQ=37)

                                                                                  aaCTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAAtt                                                             >  1:860459/1‑148 (MQ=255)

                                                                                     ttATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGc                                                          <  1:1134272/148‑1 (MQ=255)

                                                                                      tATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGcc                                                         >  1:2511390/1‑148 (MQ=255)

                                                                                        ttGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCcac                                                       <  1:948022/148‑1 (MQ=255)

                                                                                         tGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCcaca                                                      >  2:27169/1‑148 (MQ=255)

                                                                                             ttCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGa                                                  >  1:1922017/1‑148 (MQ=255)

                                                                                                  gTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGa                                             >  2:3451485/1‑148 (MQ=255)

                                                                                                   tttGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaa                                            >  2:1555592/1‑148 (MQ=255)

                                                                                                    ttGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaaa                                           >  1:3115258/1‑148 (MQ=255)

                                                                                                    ttGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGaaa                                           >  2:1119639/1‑148 (MQ=255)

                                                                                                     tGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAAt                                          <  2:1049164/148‑1 (MQ=255)

                                                                                                       aaTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGCGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCa                                        <  1:3588161/148‑1 (MQ=255)

                                                                                                       aaTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCa                                        <  1:170818/148‑1 (MQ=255)

                                                                                                          ggAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTg                                     >  1:2182854/1‑148 (MQ=255)

                                                                                                              aGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCAc                                 >  2:1325033/1‑148 (MQ=255)

                                                                                                               gggCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACg                                >  1:1092187/1‑148 (MQ=255)

                                                                                                                ggCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGt                               >  1:1971757/1‑148 (MQ=255)

                                                                                                                 gCTGAAGAAATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTg                              >  1:389792/1‑148 (MQ=255)

                                                                                                                     aagaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgaga                          >  1:3356177/1‑148 (MQ=255)

                                                                                                                     aagaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgaga                          >  1:2814773/1‑148 (MQ=255)

                                                                                                                     aagaaATTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgaga                          >  1:3023179/1‑148 (MQ=255)

                                                                                                                        aaaTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATa                       >  2:3792814/1‑148 (MQ=255)

                                                                                                                          aTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAg                     >  1:2337336/1‑148 (MQ=255)

                                                                                                                          aTTATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAg                     <  1:2502691/148‑1 (MQ=255)

                                                                                                                           ttATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTTATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGc                    <  2:3208652/148‑1 (MQ=255)

                                                                                                                            tATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCt                   >  1:3103562/1‑148 (MQ=255)

                                                                                                                            tATGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCt                   >  1:1196533/1‑148 (MQ=255)

                                                                                                                              tGCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGc                 <  1:3495730/148‑1 (MQ=255)

                                                                                                                               gCGTGGCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCa                <  1:1694099/148‑1 (MQ=255)

                                                                                                                                   ggCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTgg            <  1:1792976/148‑1 (MQ=255)

                                                                                                                                    gCAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGc           >  1:732406/1‑148 (MQ=255)

                                                                                                                                     cAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCg          >  1:543294/1‑148 (MQ=255)

                                                                                                                                     cAGTCTCCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCg          >  1:374607/1‑148 (MQ=255)

                                                                                                                                          tcCTTTAGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTgg     >  2:3547526/1‑148 (MQ=255)

GAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACG  >  CP009273/3618392‑3618680

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