BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,146,030	T→A	55.8%	S194C (AGT→TGT)	yegL ←	VMA domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,146,030	0	T	A	55.8%	135.0 / inf	299	S194C (AGT→TGT)	yegL	VMA domain protein
Reads supporting (aligned to +/- strand): ref base T (62/70); new base A (59/108); total (121/178)
Fisher's exact test for biased strand distribution p-value = 4.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 303 total aligned reads displayed.

CATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAACACTT  >  CP009273/2145881‑2146179

caTCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACt                                                                                                                                                       <  1:2685418/150‑1 (MQ=255)

     tGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATa                                                                                                                                                  <  1:172642/150‑1 (MQ=255)

          aGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAg                                                                                                                                             <  1:1746564/150‑1 (MQ=255)

            gCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCt                                                                                                                                           <  1:3957487/150‑1 (MQ=255)

              aaGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTaa                                                                                                                                         <  2:4233355/150‑1 (MQ=255)

               aGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGACACCGAACGCAGCGAACAGGATAACCCGCTaaa                                                                                                                                        >  2:2469190/1‑150 (MQ=255)

               aGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTaaa                                                                                                                                        >  1:2429113/1‑150 (MQ=255)

               aGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTaaa                                                                                                                                        >  2:3840760/1‑150 (MQ=255)

                    gcgcTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGtt                                                                                                                                   <  1:1232342/150‑1 (MQ=255)

                     cgcTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTc                                                                                                                                  >  2:932558/1‑150 (MQ=255)

                      gcTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCa                                                                                                                                 <  2:4294637/150‑1 (MQ=255)

                        tGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACg                                                                                                                               <  1:1095853/150‑1 (MQ=255)

                        tGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACg                                                                                                                               >  1:630580/1‑150 (MQ=255)

                        tGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACg                                                                                                                               <  1:2531648/150‑1 (MQ=255)

                        tGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACg                                                                                                                               >  1:2282848/1‑150 (MQ=255)

                         gATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACgg                                                                                                                              <  2:4216454/150‑1 (MQ=255)

                              gCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTg                                                                                                                         <  1:1937780/150‑1 (MQ=255)

                                  cGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAAt                                                                                                                     >  1:1585849/1‑150 (MQ=255)

                                   gTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATc                                                                                                                    >  2:3755278/1‑150 (MQ=255)

                                    ttCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATcc                                                                                                                   <  2:2011408/150‑1 (MQ=255)

                                     tCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGACCTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCt                                                                                                                  <  2:1051984/150‑1 (MQ=255)

                                     tCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCt                                                                                                                  <  2:363940/150‑1 (MQ=255)

                                      ccAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCtt                                                                                                                 >  2:1301951/1‑150 (MQ=255)

                                      ccAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCtt                                                                                                                 >  2:159834/1‑150 (MQ=255)

                                       cAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTg                                                                                                                >  1:3850135/1‑150 (MQ=255)

                                        aaCGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGc                                                                                                               >  1:2252675/1‑150 (MQ=255)

                                        aaCGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGc                                                                                                               >  1:4586888/1‑150 (MQ=255)

                                        aaCGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGc                                                                                                               >  2:406945/1‑150 (MQ=255)

                                        aaCGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGc                                                                                                               <  1:1889900/150‑1 (MQ=255)

                                         aCGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCa                                                                                                              <  2:4413199/150‑1 (MQ=255)

                                           ggTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCaga                                                                                                            >  2:2771109/1‑150 (MQ=255)

                                              tGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGc                                                                                                         <  2:1670366/150‑1 (MQ=255)

                                                 tGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCaaa                                                                                                      <  1:1027354/150‑1 (MQ=255)

                                                 tGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCaaa                                                                                                      <  1:2378829/150‑1 (MQ=255)

                                                  gCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAg                                                                                                     <  2:2042577/150‑1 (MQ=255)

                                                      aGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTg                                                                                                 <  1:3933235/150‑1 (MQ=255)

                                                        aCCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAGCCGAACGCAGCGAACTGGATAACCAGCTACACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGAc                                                                                               <  2:256425/150‑1 (MQ=255)

                                                        aCCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGAc                                                                                               <  2:3787929/150‑1 (MQ=255)

                                                        aCCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGAc                                                                                               <  1:389858/150‑1 (MQ=255)

                                                        aCCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGAc                                                                                               <  1:4296080/150‑1 (MQ=255)

                                                         ccAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACg                                                                                              >  2:637824/1‑150 (MQ=255)

                                                         ccAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACg                                                                                              >  2:2021402/1‑150 (MQ=255)

                                                           aGACGCCAGCTCACACTGACGTCCAACCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGaa                                                                                            >  2:770568/1‑150 (MQ=255)

                                                            gACTCCAGCTCACACTGTCGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCTTTCCCGGCGTGGAGCGGGAAACCGAACTCAGCGAACAGGATAACCATCTAATCAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAAc                                                                                           <  2:2251360/150‑1 (MQ=255)

                                                            gACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAAc                                                                                           <  2:3165482/150‑1 (MQ=255)

                                                             aCGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACg                                                                                          <  2:1195491/150‑1 (MQ=255)

                                                             aCGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACg                                                                                          <  1:4452866/150‑1 (MQ=255)

                                                                  aGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAAt                                                                                     <  2:2487480/150‑1 (MQ=255)

                                                                  aGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAAt                                                                                     <  2:3967380/150‑1 (MQ=255)

                                                                    cTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAAttt                                                                                   <  2:1955218/150‑1 (MQ=255)

                                                                    cTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAAttt                                                                                   <  1:1298525/150‑1 (MQ=255)

                                                                    cTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAAttt                                                                                   <  1:115649/150‑1 (MQ=255)

                                                                     tCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTtg                                                                                  >  1:3686666/1‑150 (MQ=255)

                                                                     tCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATAtg                                                                                  >  1:3455335/1‑150 (MQ=255)

                                                                      cacaCTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTtgt                                                                                 >  1:4100723/1‑150 (MQ=255)

                                                                      cacaCTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTtgt                                                                                 >  2:4457139/1‑150 (MQ=255)

                                                                      cacaCTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCGAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTtgt                                                                                 >  2:3576233/1‑150 (MQ=255)

                                                                      cacaCTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCATCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTtgt                                                                                 <  2:4331158/150‑1 (MQ=255)

                                                                           tGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAg                                                                            >  1:3802831/1‑150 (MQ=255)

                                                                               gTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCTAACGCAGCGAACTGGATAACCAGCTAAACAGTTCAAGGAATTGTAATACTTGCAGAGGCAAAGGCTGACGAACGATAATTTGTGACAGCGTc                                                                        >  1:671622/1‑150 (MQ=255)

                                                                                    gCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCat                                                                   <  2:1744862/150‑1 (MQ=255)

                                                                                      cTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCatat                                                                 >  2:2930829/1‑150 (MQ=255)

                                                                                      cTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCatat                                                                 >  1:876541/1‑150 (MQ=255)

                                                                                           ggCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGcc                                                            <  1:3627958/150‑1 (MQ=255)

                                                                                           ggCGCCTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGcc                                                            <  1:2697862/150‑1 (MQ=255)

                                                                                            gcgcTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGccc                                                           <  1:3158281/150‑1 (MQ=255)

                                                                                              gcTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGccccc                                                         >  1:4641989/1‑150 (MQ=255)

                                                                                                ttCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTg                                                       <  2:2802399/150‑1 (MQ=255)

                                                                                                    aGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAaca                                                   <  2:1342700/150‑1 (MQ=255)

                                                                                                       acaACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCg                                                >  1:3460339/1‑150 (MQ=255)

                                                                                                         aaCTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGAt                                              <  1:2983040/150‑1 (MQ=255)

                                                                                                         aaCTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGAt                                              >  1:1947963/1‑150 (MQ=255)

                                                                                                            ttCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGa                                           <  2:3217136/150‑1 (MQ=255)

                                                                                                            ttCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGa                                           <  1:4561775/150‑1 (MQ=255)

                                                                                                             tCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGaa                                          <  2:4529599/150‑1 (MQ=255)

                                                                                                                gTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaag                                       <  2:1157083/150‑1 (MQ=255)

                                                                                                                  tCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaagaa                                     <  1:3555460/150‑1 (MQ=255)

                                                                                                                  tCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaagaa                                     <  2:3684041/150‑1 (MQ=255)

                                                                                                                    ccGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAaaagaaag                                   <  1:2217371/150‑1 (MQ=255)

                                                                                                                        cGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCaaa                               <  2:4175909/150‑1 (MQ=255)

                                                                                                                          tGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAAcc                             >  1:3271071/1‑150 (MQ=255)

                                                                                                                          tGGAGCGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAAcc                             >  1:3979319/1‑150 (MQ=255)

                                                                                                                              gCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCtttt                         <  1:2038762/150‑1 (MQ=255)

                                                                                                                               cGGGAAACCGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTa                        >  1:388740/1‑150 (MQ=255)

                                                                                                                                   aaaCCGAACGCAGCGACATGGATAACCAGCTAAACTGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGTACACCGATGGAAAAGATAGCAAACCTTTTAtctt                    <  1:1499073/150‑1 (MQ=255)

                                                                                                                                    aaCCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtcttc                   >  2:197852/1‑150 (MQ=255)

                                                                                                                                     aCCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTAtcttct                  <  2:2429113/150‑1 (MQ=255)

                                                                                                                                      ccGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATcttctt                 <  1:3370852/150‑1 (MQ=255)

                                                                                                                                      ccGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATcttctt                 <  1:4180784/150‑1 (MQ=255)

                                                                                                                                       cGAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCttcttc                <  1:2010166/150‑1 (MQ=255)

                                                                                                                                        gAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTcc               >  2:4317298/1‑150 (MQ=255)

                                                                                                                                        gAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTcc               >  2:4465758/1‑150 (MQ=255)

                                                                                                                                        gAACGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTcc               <  1:1605827/150‑1 (MQ=255)

                                                                                                                                         aaCGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTccc              <  2:1624169/150‑1 (MQ=255)

                                                                                                                                           cGCAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTccccc            <  2:1490747/150‑1 (MQ=255)

                                                                                                                                             cAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCg          <  1:637824/150‑1 (MQ=255)

                                                                                                                                             cAGCGAACAGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCg          <  1:2021449/150‑1 (MQ=255)

                                                                                                                                                  aaCTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAaca     >  2:4165627/1‑150 (MQ=255)

                                                                                                                                                   aCTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAacac    <  1:2168412/150‑1 (MQ=255)

                                                                                                                                                    cTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAACACt   >  2:3295839/1‑150 (MQ=255)

CATCCTGACAAGGCAAGTCCGCGCTGATGTGCGACGTTCCAACGGTTGATGCATAGACCAGACGCCAGCTCACACTGACGTCCAGCCTTTTGGCGCTTCCAGCACAACTTCCGTTCCCGGCGTGGAGCGGGAAACCGAACGCAGCGAACTGGATAACCAGCTAAACAGTTCACGGAATTGTAATCCTTGCAGAGGCAAAGGCTGACGAACGCTAATTTGTGCCAGCGTCTTCATATCCGCCCCCTGAACACCGATGGAAAAGAAAGCAAACCTTTTATCTTCTTCCCCCCGAAACACTT  >  CP009273/2145881‑2146179

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