BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,994,247	A→T	92.9%	V57E (GTG→GAG)	yigE ←	DUF2233 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,994,247	0	A	T	92.9%	836.4 / 41.7	254	V57E (GTG→GAG)	yigE	DUF2233 family protein
Reads supporting (aligned to +/- strand): ref base A (3/15); new base T (107/129); total (110/144)
Fisher's exact test for biased strand distribution p-value = 2.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 254 total aligned reads displayed.

CGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTA  >  CP009273/3994101‑3994393

cGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTccc                                                                                                                                                 <  1:929151/150‑1 (MQ=255)

   tGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGct                                                                                                                                              <  1:1963475/150‑1 (MQ=255)

   tGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATGAATAGCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGct                                                                                                                                              <  2:4431258/150‑1 (MQ=255)

    gTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGctc                                                                                                                                             >  1:4165637/1‑150 (MQ=255)

    gTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGctc                                                                                                                                             <  1:3371953/150‑1 (MQ=255)

        aaaCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCCTTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGtt                                                                                                                                         >  2:1224637/1‑150 (MQ=255)

        aaaCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGtt                                                                                                                                         >  2:1019157/1‑150 (MQ=255)

         aaCCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGttt                                                                                                                                        >  2:2232833/1‑150 (MQ=255)

             gAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAgg                                                                                                                                    >  2:3446723/1‑150 (MQ=255)

               gcggcgCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGAt                                                                                                                                  <  2:115289/150‑1 (MQ=255)

                   cgcATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAAc                                                                                                                              >  2:3312103/1‑150 (MQ=255)

                     cATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACgg                                                                                                                            >  1:4272841/1‑150 (MQ=255)

                        aGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGtat                                                                                                                         <  1:2005092/150‑1 (MQ=255)

                                 catAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTa                                                                                                                <  2:2694369/150‑1 (MQ=255)

                                 catAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTa                                                                                                                <  2:4485159/150‑1 (MQ=255)

                                  atAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTAc                                                                                                               >  2:2284863/1‑150 (MQ=255)

                                     gATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGt                                                                                                            <  2:667673/150‑1 (MQ=255)

                                     gATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGt                                                                                                            >  2:1467296/1‑150 (MQ=255)

                                          cgccgTTCATCGCCATCTGCACCTGACCCTGACTATTAAGATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCg                                                                                                       <  2:4217967/150‑1 (MQ=255)

                                            ccgTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTc                                                                                                     <  2:3053861/150‑1 (MQ=255)

                                                tCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGAt                                                                                                 <  1:3482407/150‑1 (MQ=255)

                                                tCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGAt                                                                                                 <  1:3482625/150‑1 (MQ=255)

                                                  aTCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCt                                                                                               >  2:3161883/1‑150 (MQ=255)

                                                     gCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTgag                                                                                            <  2:583552/150‑1 (MQ=255)

                                                         tCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTg                                                                                        <  1:374131/150‑1 (MQ=255)

                                                           tGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTgcg                                                                                      <  2:1129281/150‑1 (MQ=255)

                                                                  gtccctgacTATTTATATCCGCCAGCAGAGCTTGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACTTTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTTCGGTCAGCGTCGGTTCTGAGAGTGCGCAatctt                                                                               <  2:1274938/148‑3 (MQ=255)

                                                                   accctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAatcatc                                                                              <  1:2872902/150‑1 (MQ=255)

                                                                   accctgacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAatcatc                                                                              <  1:1701993/150‑1 (MQ=255)

                                                                      ctgacTATTAATATCCGCCAGCAGAGCATGTAACGGTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGc                                                                           <  2:4275340/150‑1 (MQ=255)

                                                                        gacTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCgg                                                                         <  2:3714250/150‑1 (MQ=255)

                                                                         acTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGc                                                                        >  1:2879990/1‑150 (MQ=255)

                                                                          cTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCa                                                                       >  1:2583394/1‑150 (MQ=255)

                                                                          cTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCa                                                                       <  2:1390203/150‑1 (MQ=255)

                                                                           tATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCaa                                                                      <  1:1982322/150‑1 (MQ=255)

                                                                           tATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCaa                                                                      <  2:1564557/150‑1 (MQ=255)

                                                                             ttAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACa                                                                    <  1:4184732/150‑1 (MQ=255)

                                                                              tAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAg                                                                   <  2:244570/150‑1 (MQ=255)

                                                                               aaTATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGc                                                                  >  2:1288999/1‑150 (MQ=255)

                                                                               aaTATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTGTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGc                                                                  >  2:3399540/1‑150 (MQ=255)

                                                                                atatCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCa                                                                 >  1:1341585/1‑150 (MQ=255)

                                                                                     cGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACaa                                                            >  2:1973185/1‑150 (MQ=255)

                                                                                        cagcagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAggg                                                         >  2:2985139/1‑150 (MQ=255)

                                                                                        cagcagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAggg                                                         >  1:2720430/1‑150 (MQ=255)

                                                                                         agcagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGa                                                        >  1:3219771/1‑150 (MQ=255)

                                                                                         agcagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGa                                                        <  2:2627562/150‑1 (MQ=255)

                                                                                         agcagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGa                                                        <  2:2738742/150‑1 (MQ=255)

                                                                                         agcagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGa                                                        >  2:3386594/1‑150 (MQ=255)

                                                                                          gcagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGaa                                                       >  2:1226877/1‑150 (MQ=255)

                                                                                           cagAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGACAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGAGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAg                                                      >  2:906139/1‑150 (MQ=255)

                                                                                            agagCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGt                                                     >  1:687493/1‑150 (MQ=255)

                                                                                            agagCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGt                                                     <  1:3786458/150‑1 (MQ=255)

                                                                                            agagCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGt                                                     >  2:323580/1‑150 (MQ=255)

                                                                                                 aTGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGt                                                <  2:4069153/150‑1 (MQ=255)

                                                                                                 aTGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGt                                                >  2:424649/1‑150 (MQ=255)

                                                                                                   gTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGa                                              >  1:4418450/1‑150 (MQ=255)

                                                                                                   gTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGa                                              >  1:1640063/1‑150 (MQ=255)

                                                                                                   gTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGa                                              >  1:3897366/1‑150 (MQ=255)

                                                                                                   gTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGa                                              <  2:1596663/150‑1 (MQ=255)

                                                                                                    tAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAg                                             <  1:3486785/150‑1 (MQ=255)

                                                                                                    tAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAg                                             <  2:1167319/150‑1 (MQ=255)

                                                                                                       cGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGc                                          <  2:1401508/150‑1 (MQ=255)

                                                                                                       cGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGc                                          <  2:1526622/150‑1 (MQ=255)

                                                                                                        gTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCg                                         >  1:2941878/1‑150 (MQ=255)

                                                                                                         ttCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCga                                        <  2:1889575/150‑1 (MQ=255)

                                                                                                          tCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCgag                                       >  1:1719161/1‑150 (MQ=255)

                                                                                                           ccccACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCgaga                                      <  1:1495900/150‑1 (MQ=255)

                                                                                                                 gCTTCGCCTTTGGCTTTGTGCCTGTTCATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATTCGCCTGTACGGTCAGCTTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGTAAAATc                                <  2:565076/150‑1 (MQ=255)

                                                                                                                   ttCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTTAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCg                              >  1:4251120/1‑150 (MQ=255)

                                                                                                                   ttCGCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCg                              <  1:2232833/150‑1 (MQ=255)

                                                                                                                      gCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGttt                           <  2:4165712/150‑1 (MQ=255)

                                                                                                                      gCCATTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGttt                           >  1:1075/1‑150 (MQ=255)

                                                                                                                         aTTGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTgag                        <  1:2030486/150‑1 (MQ=255)

                                                                                                                          ttGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTgaga                       >  1:2889697/1‑150 (MQ=255)

                                                                                                                           tGGCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGAt                      <  1:3665192/150‑1 (MQ=255)

                                                                                                                             gCTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTc                    <  2:2889697/150‑1 (MQ=255)

                                                                                                                              cTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCa                   >  1:3201917/1‑150 (MQ=255)

                                                                                                                              cTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCa                   >  1:4021443/1‑150 (MQ=255)

                                                                                                                              cTTTTTGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCa                   <  1:1825601/150‑1 (MQ=255)

                                                                                                                              atttttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGTTTCa                   <  1:4077601/149‑1 (MQ=255)

                                                                                                                               tttttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCaa                  <  1:4040798/150‑1 (MQ=255)

                                                                                                                               tttttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCaa                  <  1:1789513/150‑1 (MQ=255)

                                                                                                                                ttttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAg                 >  2:4003793/1‑150 (MQ=255)

                                                                                                                                ttttGCCAGTACATGTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAg                 <  2:2707605/150‑1 (MQ=255)

                                                                                                                                 tttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAgg                <  1:1119548/150‑1 (MQ=255)

                                                                                                                                  ttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGt               >  2:2787548/1‑150 (MQ=255)

                                                                                                                                  ttGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGt               >  2:2786900/1‑150 (MQ=255)

                                                                                                                                   tGCCAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTg              >  2:2121504/1‑150 (MQ=255)

                                                                                                                                      cAGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATc           >  1:492410/1‑150 (MQ=255)

                                                                                                                                       aGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCa          >  1:1012971/1‑150 (MQ=255)

                                                                                                                                       aGTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCa          >  1:721261/1‑150 (MQ=255)

                                                                                                                                        gTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCAt         <  2:3679543/150‑1 (MQ=255)

                                                                                                                                        gTACATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCAt         >  1:4436874/1‑150 (MQ=255)

                                                                                                                                          aCATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTc       >  1:1576785/1‑150 (MQ=255)

                                                                                                                                           cATTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTcc      <  1:2451914/150‑1 (MQ=255)

                                                                                                                                            aTTTTCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCt     >  1:2509543/1‑150 (MQ=255)

                                                                                                                                             ttttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCtt    <  1:2986397/150‑1 (MQ=255)

                                                                                                                                             ttttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCtt    >  1:1048821/1‑150 (MQ=255)

                                                                                                                                              tttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCttt   >  1:723739/1‑150 (MQ=255)

                                                                                                                                               ttCTCCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTa  <  2:855398/150‑1 (MQ=255)

CGATGTACAAACCGAGCGGCGCATAGCTTTCATCATAGATGCCGCCGTTCATCGCCATCTGCACCTGACCCTGACTATTAATATCCGCCAGCAGAGCATGTAACGTTCCCCACGCTTCGCCATTGGCTTTTTGCCAGTACATTTTCACCCGCTCTGTTTGAGGATTAACGGTATACGCCTGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTA  >  CP009273/3994101‑3994393

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