BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,491,905	A→T	6.0%	intergenic (+55/+17)	ydcA → / ← hokB	protein YdcA/toxin HokB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,905	0	A	T	6.0%	320.8 / 7.7	100	intergenic (+55/+17)	ydcA/hokB	protein YdcA/toxin HokB
Reads supporting (aligned to +/- strand): ref base A (42/52); new base T (3/3); total (45/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 104 total aligned reads displayed.

GTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACACAC  >  NC_000913/1491756‑1492050

gTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGaaa                                                                                                                                                   >  1:341303/1‑150 (MQ=255)

 ttCGAAAGGGGGGATCTGACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCTTCGAAATAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGaaaa                                                                                                                                                  <  2:510005/150‑1 (MQ=255)

 ttCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGaaaa                                                                                                                                                  >  1:1037167/1‑150 (MQ=255)

   cGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATc                                                                                                                                                >  2:403563/1‑150 (MQ=255)

   cGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATc                                                                                                                                                >  2:1088090/1‑150 (MQ=255)

     aaaGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATccc                                                                                                                                              >  1:1314085/1‑150 (MQ=255)

      aaGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATcccc                                                                                                                                             <  2:72298/150‑1 (MQ=255)

          ggggATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGATTTTCCCCgcgc                                                                                                                                         <  1:1090595/150‑1 (MQ=255)

          ggggATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCgcgc                                                                                                                                         >  2:1553141/1‑150 (MQ=255)

                ctcACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGccc                                                                                                                                   >  2:959344/1‑150 (MQ=255)

                ctcACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGccc                                                                                                                                   <  2:1414538/150‑1 (MQ=255)

                 tcACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATcc                                                                                                                                               >  2:1007047/1‑137 (MQ=255)

                 tcACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATcc                                                                                                                                               <  1:1007047/137‑1 (MQ=255)

                 tcACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTACTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCt                                                                                                                                  <  1:794331/150‑1 (MQ=255)

                  cacaCTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCttttgact                                                                                                                                 <  1:1138921/150‑7 (MQ=255)

                  cacaCTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCtt                                                                                                                                 >  1:285752/1‑150 (MQ=255)

                  cacaCTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCtt                                                                                                                                 <  2:1574719/150‑1 (MQ=255)

                    cacTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTAc                                                                                                                               >  2:1277198/1‑150 (MQ=255)

                    cacTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCACGCGCTTGCCCTTAc                                                                                                                               >  1:732058/1‑150 (MQ=255)

                       tGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTg                                                                                                                            >  1:507943/1‑150 (MQ=255)

                       tGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTg                                                                                                                            >  1:126407/1‑150 (MQ=255)

                          aCGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGAc                                                                                                                         >  1:1220858/1‑150 (MQ=255)

                           cGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACg                                                                                                                        >  1:12196/1‑150 (MQ=255)

                             gCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTg                                                                                                                      >  2:436256/1‑150 (MQ=255)

                             gCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTg                                                                                                                      >  2:535406/1‑150 (MQ=255)

                              cAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGc                                                                                                                     <  2:109631/150‑1 (MQ=255)

                              cAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGc                                                                                                                     <  2:206063/150‑1 (MQ=255)

                                 gcggcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAgg                                                                                                                  <  2:1336368/150‑1 (MQ=255)

                                   ggcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGcc                                                                                                                >  1:1140973/1‑150 (MQ=255)

                                    gcAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCa                                                                                                               <  1:1277198/150‑1 (MQ=255)

                                     cAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCAt                                                                                                              <  1:1612314/150‑1 (MQ=255)

                                      aaaTTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATg                                                                                                             <  2:1329451/150‑1 (MQ=255)

                                          ttGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAgcg                                                                                                         <  1:403563/150‑1 (MQ=255)

                                             tCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATCCACTAACTGAAGTGTAAAAGGGGGGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgcag                                                                                                      <  2:517594/150‑1 (MQ=255)

                                              cTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGATTGCCCTTACATGGACGTGCAGGCCATGAGCgcagc                                                                                                     >  1:860887/1‑150 (MQ=255)

                                               tGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCTCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgcagca                                                                                                    <  2:1291512/150‑1 (MQ=255)

                                               tGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgcagca                                                                                                    >  2:1516516/1‑150 (MQ=255)

                                               tGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCgcagca                                                                                                    >  2:1391123/1‑150 (MQ=255)

                                                gTATTGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCaa                                                                                                   <  2:1072724/150‑1 (MQ=255)

                                                gTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGc                                                                                                                     <  2:367872/132‑1 (MQ=255)

                                                gTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGc                                                                                                                     >  1:367872/1‑132 (MQ=255)

                                                 tAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAc                                                                                                  >  2:1566486/1‑150 (MQ=255)

                                                     gatgGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAAcctcc                                                                                              <  1:242855/150‑1 (MQ=255)

                                                            cTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCAc                                                                                       <  1:1361903/150‑1 (MQ=255)

                                                            cTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCAc                                                                                       >  1:111254/1‑150 (MQ=255)

                                                                      aTCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGaa                                                                             >  2:800796/1‑150 (MQ=255)

                                                                       tCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAAc                                                                            <  1:800796/150‑1 (MQ=255)

                                                                         gAAAAAAACATGGACTAACTGAAGTGTAAAAGCGGTGCCATGAGTTAATGAATCTGCTAGTCAAATGCGCACGGAAAATCCCCGCGCTTGCCCTAACCTGGACGTGCAGGGCATGAGCGCGGCAACCTCCTTATCACCGTCCCGGAAccc                                                                          >  2:1601281/1‑149 (MQ=255)

                                                                         gAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCg                                                                          >  2:1392265/1‑150 (MQ=255)

                                                                         gAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCg                                                                          <  2:1314085/150‑1 (MQ=255)

                                                                          aaaaaaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGc                                                                         <  1:557677/150‑1 (MQ=255)

                                                                             aaaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGt                                                                      >  1:499870/1‑150 (MQ=255)

                                                                              aaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGtt                                                                     >  2:316149/1‑150 (MQ=255)

                                                                              aaaCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGtt                                                                     <  1:331397/150‑1 (MQ=255)

                                                                                aCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCg                                                                   <  2:1624855/150‑1 (MQ=255)

                                                                                aCATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCg                                                                   <  1:1474777/150‑1 (MQ=255)

                                                                                  aTGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTa                                                                 <  1:436256/150‑1 (MQ=255)

                                                                                   tGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTag                                                                <  2:1272881/150‑1 (MQ=255)

                                                                                     cACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTagag                                                              >  2:817468/1‑150 (MQ=255)

                                                                                      aCTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGc                                                             <  2:1140973/150‑1 (MQ=255)

                                                                                             gAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGCCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTAGTc                                                      >  1:503086/1‑150 (MQ=255)

                                                                                               aGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTTTCgg                                                    >  1:1481188/1‑150 (MQ=255)

                                                                                                gtgtAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCggg                                                   >  2:53946/1‑150 (MQ=255)

                                                                                                gtgtAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCggg                                                   >  2:60033/1‑150 (MQ=255)

                                                                                                gtgtAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCggg                                                   <  1:694833/150‑1 (MQ=255)

                                                                                                gtgtAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCggg                                                   <  1:1279974/150‑1 (MQ=255)

                                                                                                   tAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCa                                                <  1:412782/150‑1 (MQ=255)

                                                                                                          ggTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGGCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAgtgt                                         <  1:1088090/150‑1 (MQ=255)

                                                                                                          ggTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAgtgt                                         >  2:1554539/1‑150 (MQ=255)

                                                                                                            tGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACGTGGACGTGCAGGCGATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGa                                       <  1:1184507/150‑1 (MQ=255)

                                                                                                                  gagaAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCa                                 <  1:1553141/150‑1 (MQ=255)

                                                                                                                   agaAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTTTCGGGTCAGGAGTGTGAATGTCAg                                <  2:1481188/150‑1 (MQ=255)

                                                                                                                   agaAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAg                                <  2:499870/150‑1 (MQ=255)

                                                                                                                     aaaaTGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGaa                              <  2:90822/150‑1 (MQ=255)

                                                                                                                       aaTGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATc                            <  2:285752/150‑1 (MQ=255)

                                                                                                                        aTGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCg                           <  1:535406/150‑1 (MQ=255)

                                                                                                                          gAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTa                         <  2:111254/150‑1 (MQ=255)

                                                                                                                            aTCTGCTAGTCAAATGCGCTGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCATGCCATGAGCGCAGCAACCTCCTTATCACCTTCCCGGTACCGCAGTTCTTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTatt                       <  1:1128856/150‑3 (MQ=255)

                                                                                                                              cTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGc                     >  1:708007/1‑150 (MQ=255)

                                                                                                                               tGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCa                    >  1:1311889/1‑150 (MQ=255)

                                                                                                                                gCTAGTCAAATGCGCTGGGACAATCCCCGCGCGGGCCCTTACCTGGCCGGGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTAGGCGGGTCAGCAGGGTGAATGTCACAATCGTAATGCAg                   <  2:506877/150‑1 (MQ=255)

                                                                                                                                gCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAg                   <  1:53946/150‑1 (MQ=255)

                                                                                                                                gCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAg                   >  2:1454877/1‑150 (MQ=255)

                                                                                                                                gCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAg                   <  1:60033/150‑1 (MQ=255)

                                                                                                                                  tAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGAt                 >  1:219114/1‑150 (MQ=255)

                                                                                                                                    gTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATaa               >  1:591758/1‑150 (MQ=255)

                                                                                                                                     gggcaagcgcGGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAAt              >  1:1341916/7‑150 (MQ=255)

                                                                                                                                      cAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATg             >  2:873928/1‑150 (MQ=255)

                                                                                                                                       gcaagcgcGGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGa            <  2:1341916/146‑1 (MQ=255)

                                                                                                                                       aaaTGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATg                      >  1:1009029/1‑140 (MQ=255)

                                                                                                                                       aaaTGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATg                      <  2:1009029/140‑1 (MQ=255)

                                                                                                                                        caagcgcGGGGATTTTCCCCGCGCTTGCACTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGa                               <  2:233161/127‑1 (MQ=255)

                                                                                                                                        caagcgcGGGGATTTTCCCCGCGCTTGCACTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGa                               >  1:233161/4‑130 (MQ=255)

                                                                                                                                        aaTGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAg           <  1:316149/150‑1 (MQ=255)

                                                                                                                                         aTGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTGGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGc          <  1:1015740/150‑1 (MQ=255)

                                                                                                                                          agcgcGGGGATTTTCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCa         >  2:1151560/2‑150 (MQ=255)

                                                                                                                                           gcgcGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAg        >  1:125099/1‑150 (MQ=255)

                                                                                                                                           gcgcGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAg        <  2:398423/150‑1 (MQ=255)

                                                                                                                                            cgcgGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGTTCAGTAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGa       <  2:1304661/150‑1 (MQ=255)

                                                                                                                                                 ggAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGacacac  >  1:691118/1‑150 (MQ=255)

GTTCGAAAGGGGGGATCTCACACTGTACGGCAGGCGGCAAATTTGTCTGTAATGATGGTTCTATTAGTGCATCGAAAAAAACATGCACTAACTGAAGTGTAAAAGGGGTGCCATGAGAAAATGAATCTGCTAGTCAAATGCGCGGGGAAAATCCCCGCGCTTGCCCTTACCTGGACGTGCAGGCCATGAGCGCAGCAACCTCCTTATCACCGTCCCGGAACCGCAGTTCGTAGAGCGTTTGTCGGGTCAGGAGTGTGAATGTCAGAATCGTAATGCAGATAATGAGCAGACACAC  >  NC_000913/1491756‑1492050

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