BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,502,848	T→A	10.3%	V604E (GTA→GAA)	ydcP →	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,502,848	0	T	A	10.3%	1285.5 / 105.5	436	V604E (GTA→GAA)	ydcP	putative peptidase
Reads supporting (aligned to +/- strand): ref base T (200/190); new base A (24/21); total (225/211)
Fisher's exact test for biased strand distribution p-value = 8.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 446 total aligned reads displayed.

GTATGAAGCACATGAAGAGAAGGGCGAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGT  >  CP009273/1502699‑1502997

     aaGCACATGAAGAGAAGGGCGAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTaa                                                                                                                                                  >  2:33416/1‑150 (MQ=255)

        cacaTGTAGAGTAGGGCGAAGTCCCTGTGATGATCACCAAGCATTGTCTTCGCTTTGCCATTAATCTGTGCCCTTAACAGGCGAAAGGCAATATCATAAGCTGGAAGGCGACGACAATGCAACTGGTTAACGGCGATGAAGAATTAACGc                                                                                                                                               <  2:3467013/150‑1 (MQ=255)

          caTGAAGAGAAGGGCGAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTa                                                                                                                                             >  1:4867522/1‑150 (MQ=255)

          caTGAAGAGAAGGGCGAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGAATTAACGCTa                                                                                                                                             >  2:2695486/1‑150 (MQ=255)

          caTGAAGAGAAGGGCGAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGATGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTa                                                                                                                                             >  2:4430565/1‑150 (MQ=255)

           aTGAAGAGAAGGGCGAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTaa                                                                                                                                            <  2:1786783/150‑1 (MQ=255)

             gAAGAGAAGGGCGAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAg                                                                                                                                          <  1:4654670/150‑1 (MQ=255)

                   aaGGGCGAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGAt                                                                                                                                    <  1:2272143/150‑1 (MQ=255)

                         gAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTgccgc                                                                                                                              >  1:2703155/1‑150 (MQ=255)

                         gAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGAATTAACGCTAAAGTTTGATTgccgc                                                                                                                              >  1:3553434/1‑150 (MQ=255)

                           aGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGccc                                                                                                                            >  1:830111/1‑150 (MQ=255)

                           aGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGAATTAACGCTAAAGTTTGATTgccgggg                                                                                                                            <  2:3076874/150‑4 (MQ=255)

                              cccGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATg                                                                                                                         <  2:2153376/150‑1 (MQ=255)

                                 ggTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCga                                                                                                                      <  2:4698827/150‑1 (MQ=255)

                                   tgatgaTCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCgaga                                                                                                                    >  2:483344/1‑150 (MQ=255)

                                   tgatgaTCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCgaga                                                                                                                    >  1:4404883/1‑150 (MQ=255)

                                    gatgatCACCAAGCATTGTCTGCGCGTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGAt                                                                                                                   >  2:154472/1‑150 (MQ=255)

                                      tgatCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGAATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGc                                                                                                                 >  1:2535549/1‑150 (MQ=255)

                                        atCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCAc                                                                                                               <  2:4019092/150‑1 (MQ=255)

                                        atCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCAc                                                                                                               <  1:1910412/150‑1 (MQ=255)

                                        atCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGAATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCAc                                                                                                               >  1:4040964/1‑150 (MQ=255)

                                         tCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACg                                                                                                              >  2:3704907/1‑150 (MQ=255)

                                         tCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACg                                                                                                              <  2:2923479/150‑1 (MQ=255)

                                         tCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACg                                                                                                              >  1:216933/1‑150 (MQ=255)

                                           aCCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGAATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTc                                                                                                            <  1:1080826/150‑1 (MQ=255)

                                            ccAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGAATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCa                                                                                                           <  2:5291139/150‑1 (MQ=255)

                                              aaGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCAtt                                                                                                         <  2:2433663/150‑1 (MQ=255)

                                              aaGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCAtt                                                                                                         <  2:958940/150‑1 (MQ=255)

                                              aaGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCAtt                                                                                                         >  1:2253919/1‑150 (MQ=255)

                                               aGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTg                                                                                                        <  2:568901/150‑1 (MQ=255)

                                                 cATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGc                                                                                                      >  1:1223166/1‑150 (MQ=255)

                                                  aTTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCa                                                                                                     >  2:3962467/1‑150 (MQ=255)

                                                    tGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCaaa                                                                                                   <  1:1657039/150‑1 (MQ=255)

                                                     gTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCaaaa                                                                                                  >  1:2210323/1‑150 (MQ=255)

                                                      tCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCaaaat                                                                                                 <  1:714188/150‑1 (MQ=255)

                                                           gcTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCaaaatcaaaa                                                                                            <  1:2663496/150‑1 (MQ=255)

                                                           gcTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCaaaatcaaaa                                                                                            >  2:1184609/1‑150 (MQ=255)

                                                            cTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCaaaatcaaaaa                                                                                           <  2:4545380/150‑1 (MQ=255)

                                                            cTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCaaaatcaaaaa                                                                                           >  1:2584212/1‑150 (MQ=255)

                                                             tttGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCCTGCGAGATGCACGTCATTGGCaaaatcaaaaat                                                                                          >  2:339397/1‑150 (MQ=255)

                                                                  cattAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACAt                                                                                     <  1:4242706/148‑1 (MQ=255)

                                                                      aaTCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTg                                                                                 >  2:3538009/1‑150 (MQ=255)

                                                                       aTCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGa                                                                                >  1:4519390/1‑150 (MQ=255)

                                                                       aTCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGa                                                                                >  1:4981206/1‑150 (MQ=255)

                                                                         cTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGaaa                                                                              >  2:315198/1‑150 (MQ=255)

                                                                            tgCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATg                                                                           <  2:619475/150‑1 (MQ=255)

                                                                              cccGAAACAGGCGAAAGGCAATATCAAATGCTGGAAGGTGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGcc                                                                         <  2:5125184/150‑1 (MQ=255)

                                                                              cccGAAACAGGAGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGAATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGcc                                                                         <  2:2575162/150‑1 (MQ=255)

                                                                               ccGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCg                                                                        <  2:1034377/150‑1 (MQ=255)

                                                                               ccGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCg                                                                        <  2:1281926/150‑1 (MQ=255)

                                                                               ccGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCg                                                                        <  2:3993956/150‑1 (MQ=255)

                                                                               ccGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCg                                                                        <  1:1747019/150‑1 (MQ=255)

                                                                                 gAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGAATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGtt                                                                      >  1:5279225/1‑150 (MQ=255)

                                                                                  aaaCAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTa                                                                     <  1:1348903/150‑1 (MQ=255)

                                                                                  aaaCAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTa                                                                     <  2:1651891/150‑1 (MQ=255)

                                                                                     cAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCg                                                                  >  2:1356588/1‑150 (MQ=255)

                                                                                          gAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAg                                                             <  2:4918712/150‑1 (MQ=255)

                                                                                            aaGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCg                                                           >  2:920512/1‑150 (MQ=255)

                                                                                            aaGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCg                                                           >  1:3328051/1‑150 (MQ=255)

                                                                                            aaGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCg                                                           <  1:3970110/150‑1 (MQ=255)

                                                                                               gCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAg                                                        <  1:197324/150‑1 (MQ=255)

                                                                                                cAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGt                                                       <  1:4312491/150‑1 (MQ=255)

                                                                                                  atatCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTgg                                                     >  1:2057145/1‑150 (MQ=255)

                                                                                                   tatCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGc                                                    <  2:4244742/150‑1 (MQ=255)

                                                                                                      cAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATc                                                 >  2:766614/1‑150 (MQ=255)

                                                                                                       aaaaGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCACCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATcc                                                >  2:1423769/1‑150 (MQ=255)

                                                                                                            cTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAg                                           >  1:2986394/1‑150 (MQ=255)

                                                                                                             tGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGt                                          >  1:3482271/1‑150 (MQ=255)

                                                                                                             tGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGt                                          >  1:4212237/1‑150 (MQ=255)

                                                                                                               gAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTcc                                        <  1:4623382/150‑1 (MQ=255)

                                                                                                                aaGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCg                                       >  2:4106502/1‑150 (MQ=255)

                                                                                                                aaGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGAATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCg                                       >  2:1073030/1‑150 (MQ=255)

                                                                                                                   gCGACGCCAATGCAACTGGTTAATGGCGATGAAGTATTAACGCTAAAGTTTGAGTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGAt                                    <  1:897148/150‑1 (MQ=255)

                                                                                                                   gCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGAt                                    <  2:229850/150‑1 (MQ=255)

                                                                                                                   gCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGAt                                    <  2:2265047/150‑1 (MQ=255)

                                                                                                                    cGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATg                                   <  1:1132234/150‑1 (MQ=255)

                                                                                                                     gACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGa                                  >  2:3046147/1‑150 (MQ=255)

                                                                                                                      aCGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAg                                 <  1:385466/150‑1 (MQ=255)

                                                                                                                        gCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAgct                               >  2:2958437/1‑150 (MQ=255)

                                                                                                                        gCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAgct                               <  1:2412516/150‑1 (MQ=255)

                                                                                                                          cAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAgctgc                             >  2:385134/1‑150 (MQ=255)

                                                                                                                          cAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAgctgc                             >  2:1937296/1‑150 (MQ=255)

                                                                                                                           aaTGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGACGCGTAGTGGCATCCGTAAGTCCGGATGAgctgct                            >  2:2051991/1‑150 (MQ=255)

                                                                                                                           aaTGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAgctgct                            >  1:1522969/1‑150 (MQ=255)

                                                                                                                           aaTGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAgctgct                            >  1:4123473/1‑150 (MQ=255)

                                                                                                                                aaCTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGaaaa                       >  1:323454/1‑150 (MQ=255)

                                                                                                                                 aCTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAAc                      >  1:2322623/1‑150 (MQ=255)

                                                                                                                                  cTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACa                     <  1:1016500/150‑1 (MQ=255)

                                                                                                                                   tGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACAt                    <  2:519085/150‑1 (MQ=255)

                                                                                                                                   tGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACAt                    >  1:3708287/1‑150 (MQ=255)

                                                                                                                                    ggTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACAtt                   <  2:4867522/150‑1 (MQ=255)

                                                                                                                                     gTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTg                  <  1:2941805/150‑1 (MQ=255)

                                                                                                                                        aaCGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCg               >  1:5124633/1‑150 (MQ=255)

                                                                                                                                        aaCGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGAAAGTCCGGATGAGCTGCTGAAAACATTGCCg               >  1:1179297/1‑150 (MQ=255)

                                                                                                                                             cGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTTCCGGGTAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCg          <  2:2236451/150‑1 (MQ=255)

                                                                                                                                              gATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGa         >  1:572800/1‑150 (MQ=255)

                                                                                                                                                tGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGaaa       <  2:4725488/150‑1 (MQ=255)

                                                                                                                                                  aaGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTACGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAg     >  2:3254411/1‑150 (MQ=255)

                                                                                                                                                  aaGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAg     <  1:319011/150‑1 (MQ=255)

                                                                                                                                                   aGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAgg    <  2:2812414/150‑1 (MQ=255)

GTATGAAGCACATGAAGAGAAGGGCGAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGT  >  CP009273/1502699‑1502997

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: `‑`