BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,407,139	A→T	27.8%	C118* (TGT→TGA)	intR ←	Rac prophage; integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,407,139	0	A	T	27.8%	687.6 / inf	405	C118* (TGT→TGA)	intR	Rac prophage; integrase
Reads supporting (aligned to +/- strand): ref base A (147/145); new base T (54/59); total (201/204)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 410 total aligned reads displayed.

AGCAGGGCGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCT  >  CP009273/1406990‑1407288

 gCAGGGCGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTcac                                                                                                                                                      <  2:415454/150‑1 (MQ=255)

    gggCGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACaaa                                                                                                                                                   >  1:3468895/1‑150 (MQ=255)

     ggCGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAat                                                                                                                                                  <  2:2793617/150‑1 (MQ=255)

      gCGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAatt                                                                                                                                                 <  1:4913056/150‑1 (MQ=255)

      gCGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAatt                                                                                                                                                 <  2:3260206/150‑1 (MQ=255)

            taGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAattattatt                                                                                                                                           >  1:9024/1‑150 (MQ=255)

             aGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTc                                                                                                                                          >  1:5253599/1‑150 (MQ=255)

             aGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTc                                                                                                                                          <  2:2867727/150‑1 (MQ=255)

              gTTATCAACTGTGCGCACGGTTCGTCCTTGTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATtcttct                                                                                                                                         <  2:2355915/150‑1 (MQ=255)

                 aTCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAg                                                                                                                                      <  1:5134694/150‑1 (MQ=255)

                      cTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGttt                                                                                                                                 <  1:2927309/150‑1 (MQ=255)

                       tgtgCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGtttt                                                                                                                                <  2:4866647/150‑1 (MQ=255)

                       tgtgCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGtttt                                                                                                                                >  1:3005132/1‑150 (MQ=255)

                        gtgCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGttttt                                                                                                                               <  1:5050072/150‑1 (MQ=255)

                        gtgCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGttttt                                                                                                                               <  1:1981185/150‑1 (MQ=255)

                        gtgCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGttttt                                                                                                                               <  2:594926/150‑1 (MQ=255)

                        gtgCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGttttt                                                                                                                               <  1:271342/150‑1 (MQ=255)

                           cgcACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTATTAGTGTTTTTAAt                                                                                                                            >  1:318122/1‑150 (MQ=255)

                             cACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATtg                                                                                                                          >  1:569553/1‑150 (MQ=255)

                               cGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATtgtg                                                                                                                        <  2:1378285/150‑1 (MQ=255)

                                ggTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGTGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGa                                                                                                                       >  2:4978495/1‑150 (MQ=255)

                                ggTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGa                                                                                                                       <  2:867678/150‑1 (MQ=255)

                                 gTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGAt                                                                                                                      <  2:1181020/150‑1 (MQ=255)

                                 gTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGAGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGAt                                                                                                                      >  2:711860/1‑150 (MQ=255)

                                  ttcttcCTTTTTTGTTGGATCTTGGATTATCCATGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGAtt                                                                                                                     <  1:4789411/150‑1 (MQ=255)

                                  ttcttcCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGAtt                                                                                                                     <  2:2211479/150‑1 (MQ=255)

                                  ttcttcCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGAtt                                                                                                                     <  2:1081238/150‑1 (MQ=255)

                                    cttcCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTc                                                                                                                   >  1:994795/1‑150 (MQ=255)

                                        cTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGttt                                                                                                               >  2:2697658/1‑150 (MQ=255)

                                        cTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGttt                                                                                                               <  1:1197189/150‑1 (MQ=255)

                                          tttttGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTTTGATTTCGTTttc                                                                                                             >  1:1463383/1‑150 (MQ=255)

                                          tttgtGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTttc                                                                                                             <  2:2492041/150‑1 (MQ=255)

                                            tttGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTttctt                                                                                                           <  2:3568770/150‑1 (MQ=255)

                                            tttGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTttctt                                                                                                           <  2:1766556/150‑1 (MQ=255)

                                            tttGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTttctt                                                                                                           <  1:5319355/150‑1 (MQ=255)

                                             ttgttgGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTttcttc                                                                                                          <  1:4335317/150‑1 (MQ=255)

                                               gttgGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCAt                                                                                                        <  1:4995376/150‑1 (MQ=255)

                                                ttgGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCAtt                                                                                                       >  1:203098/1‑150 (MQ=255)

                                                     tCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGtgtt                                                                                                  >  2:1497357/1‑150 (MQ=255)

                                                      cTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGtgttt                                                                                                 >  2:261868/1‑150 (MQ=255)

                                                        tGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGtgtttgt                                                                                               >  1:2653939/1‑150 (MQ=255)

                                                         ggATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGtgtttgtt                                                                                              >  1:4125356/1‑150 (MQ=255)

                                                          gATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTttgttg                                                                                             >  1:2491553/1‑150 (MQ=255)

                                                           aTTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTtgttgt                                                                                            >  2:426404/1‑150 (MQ=255)

                                                            ttATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTgttgtt                                                                                           <  2:3005132/150‑1 (MQ=255)

                                                            ttATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTgttgtt                                                                                           <  1:63163/150‑1 (MQ=255)

                                                               tCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAAc                                                                                        >  2:4558295/1‑150 (MQ=255)

                                                                  tggTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAActct                                                                                     >  1:229554/2‑150 (MQ=255)

                                                                  aGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTGTGTTGTTAActct                                                                                     >  1:3244246/1‑150 (MQ=255)

                                                                    gTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGATTGTTTTTAATTTTGATTTCGTTTTCTTCATTGTTTTTGTTGTTAACTCTGt                                                                                   >  1:1820618/1‑150 (MQ=255)

                                                                     tAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGtc                                                                                  >  1:2251627/1‑150 (MQ=255)

                                                                     tAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTATTCATTGTGTTTGTTGTTAACTCTGtc                                                                                  >  2:4405815/1‑150 (MQ=255)

                                                                      acaGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGtct                                                                                 >  1:1265487/3‑150 (MQ=255)

                                                                      aaaGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGtct                                                                                 <  2:1706741/150‑1 (MQ=255)

                                                                          cGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAAGTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCttt                                                                             >  1:3732583/1‑150 (MQ=255)

                                                                              tctcCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAAtt                                                                         <  1:3511166/150‑1 (MQ=255)

                                                                                tcCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAAtttt                                                                       <  2:4290286/150‑1 (MQ=255)

                                                                                  cATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTg                                                                     >  2:398011/1‑150 (MQ=255)

                                                                                    tGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTc                                                                   <  2:4722953/150‑1 (MQ=255)

                                                                                     gcagcaGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTcc                                                                  <  1:4575877/150‑1 (MQ=255)

                                                                                      cagcagTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCa                                                                 <  1:2810011/150‑1 (MQ=255)

                                                                                       agcagTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAg                                                                >  2:185842/1‑150 (MQ=255)

                                                                                          agTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGata                                                             >  1:466134/1‑150 (MQ=255)

                                                                                          agTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGata                                                             >  2:5258351/1‑150 (MQ=255)

                                                                                                  ccGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAg                                                     >  1:3769759/1‑150 (MQ=255)

                                                                                                    ggTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGct                                                   <  2:4058788/150‑1 (MQ=255)

                                                                                                    ggTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGct                                                   >  2:2146125/1‑150 (MQ=255)

                                                                                                      tAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGCTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGctct                                                 <  2:1338750/150‑1 (MQ=255)

                                                                                                           tAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTAGTTTTCTTCAGTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGaa                                            >  2:3343229/1‑150 (MQ=255)

                                                                                                            aaGATATCTCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAAc                                           >  1:2074450/1‑150 (MQ=255)

                                                                                                            aaGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAAc                                           <  1:806817/150‑1 (MQ=255)

                                                                                                            aaGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAAc                                           <  1:4535301/150‑1 (MQ=255)

                                                                                                             aGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACg                                          >  1:3914245/1‑150 (MQ=255)

                                                                                                               atatCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGtt                                        >  2:4344670/1‑150 (MQ=255)

                                                                                                                tatCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGttt                                       >  2:3047937/1‑150 (MQ=255)

                                                                                                                 atCGCTATAACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGtttt                                      >  1:1547001/1‑150 (MQ=255)

                                                                                                                     cTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATg                                  <  2:4822762/150‑1 (MQ=255)

                                                                                                                       ataACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATgtc                                >  1:717541/1‑148 (MQ=255)

                                                                                                                        taACGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGAct                               <  1:4720436/150‑1 (MQ=255)

                                                                                                                        taACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGAct                               >  1:4850485/1‑150 (MQ=255)

                                                                                                                          aCGAATATGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGActct                             <  1:3661656/150‑1 (MQ=255)

                                                                                                                           cGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACtctc                            <  1:1014254/150‑1 (MQ=255)

                                                                                                                           cGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACtctc                            <  1:3261314/150‑1 (MQ=255)

                                                                                                                           cGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACtctc                            <  1:1720728/150‑1 (MQ=255)

                                                                                                                              atatGTGTTATCGGGTTACTTTCACATATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAATGCTCTCCGAACGTTTTTATGACTCTCGtt                         <  1:1293235/150‑1 (MQ=255)

                                                                                                                              atatGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGtt                         <  2:3119825/150‑1 (MQ=255)

                                                                                                                                atGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGt                       >  1:2380852/1‑150 (MQ=255)

                                                                                                                                atGTGATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGGCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGt                       >  1:256749/1‑150 (MQ=255)

                                                                                                                                    gATATCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCAcc                   >  1:190262/1‑150 (MQ=255)

                                                                                                                                      tatCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCAt                 >  1:3958792/1‑150 (MQ=255)

                                                                                                                                       atCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCAtt                >  1:3067884/1‑150 (MQ=255)

                                                                                                                                       atCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCAtt                >  2:4226475/1‑150 (MQ=255)

                                                                                                                                        tCGGGGTACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGATTCTCGTTGTCACCAttt               >  2:1953476/1‑150 (MQ=255)

                                                                                                                                         cGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCAtttt              <  1:2149961/150‑1 (MQ=255)

                                                                                                                                         cGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCAtttt              >  2:4377861/1‑150 (MQ=255)

                                                                                                                                          ggggTACTTTCACAAATTATTATTCTGGGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCAttttt             <  1:107451/150‑1 (MQ=255)

                                                                                                                                           gggTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTg            <  2:4530951/150‑1 (MQ=255)

                                                                                                                                           gggTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTg            >  2:3885487/1‑150 (MQ=255)

                                                                                                                                             gTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGcc          >  2:780158/1‑150 (MQ=255)

                                                                                                                                              tACTTTCTCAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGccc         <  2:5107111/150‑1 (MQ=255)

AGCAGGGCGATATAGTTATCAACTGTGCGCACGGTTCTTCCTTTTTTGTTGGATCTTGGATTATCCAGGTAAAGCGTTTCTCCATGCAGCAGTTCATTCCGGTAGTTTAAGATATCGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCT  >  CP009273/1406990‑1407288

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