BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	134,166	A→G	8.9%	Q851R (CAG→CGG)	acnB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	134,166	0	A	G	8.9%	324.0 / 15.8	112	Q851R (CAG→CGG)	acnB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (51/51); new base G (8/2); total (59/53)
Fisher's exact test for biased strand distribution p-value = 9.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

Only 100 of 113 total aligned reads displayed.

AATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCC  >  NC_000913/134018‑134315

aaTGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAg                                                                                                                                                       <  1:298034/150‑1 (MQ=255)

aaTGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAg                                                                                                                                                       <  2:623222/150‑1 (MQ=255)

 aTGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAAGAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGt                                                                                                                                                      <  1:1172916/150‑1 (MQ=255)

  tGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGACGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTa                                                                                                                                                     <  1:1402900/150‑1 (MQ=255)

      ttCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACAcc                                                                                                                                                 >  2:1170190/1‑150 (MQ=255)

      ttCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACAcc                                                                                                                                                 <  1:6481/150‑1 (MQ=255)

      ttCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACAcc                                                                                                                                                 <  1:839957/150‑1 (MQ=255)

       tCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCg                                                                                                                                                <  2:1123500/150‑1 (MQ=255)

        ccTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCga                                                                                                                                               >  2:1078196/1‑150 (MQ=255)

          tGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCgaga                                                                                                                                             >  1:1040606/1‑150 (MQ=255)

               tCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCGGTGTACTGAGAAAGct                                                                                                                                        <  2:84533/150‑2 (MQ=255)

                cTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCg                                                                                                                                       >  2:1051955/1‑150 (MQ=255)

                  gCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCGGTGTACTGAGAAAGctggt                                                                                                                                     <  1:623911/150‑5 (MQ=255)

                  gCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGAt                                                                                                                                     >  2:877070/1‑150 (MQ=255)

                  gCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGAt                                                                                                                                     >  1:1534901/1‑150 (MQ=255)

                   cGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATg                                                                                                                                    <  1:1278766/150‑1 (MQ=255)

                      aaCTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATgggg                                                                                                                                 <  1:1286574/150‑1 (MQ=255)

                       aCTGGCGGCGGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGt                                                                                                                                <  2:1615474/150‑1 (MQ=255)

                          ggcggcTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGAt                                                                                                                             <  1:827240/150‑1 (MQ=255)

                                tGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCa                                                                                                                       <  1:87372/150‑1 (MQ=255)

                                tGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCa                                                                                                                       <  1:216230/150‑1 (MQ=255)

                                 gTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCATCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAg                                                                                                                      <  2:1530974/150‑1 (MQ=255)

                                 gTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAg                                                                                                                      <  1:697794/150‑1 (MQ=255)

                                 gTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAg                                                                                                                      <  1:1193364/150‑1 (MQ=255)

                                      ggcgCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGc                                                                                                                 <  2:303598/150‑1 (MQ=255)

                                        cgcTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCgg                                                                                                               <  1:909481/150‑1 (MQ=255)

                                              ttGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAg                                                                                                                      >  2:292132/1‑137 (MQ=255)

                                              ttGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAg                                                                                                                      <  1:292132/137‑1 (MQ=255)

                                               tGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTaaaa                                                                                                        >  1:761066/1‑150 (MQ=255)

                                                    aaCTGGCGACGTCGTAAGAGTACCAGACCTACGGGGCGCAGGTTGATATATGAGCCGTTGATACTTTCCGTTATCTGATCTTCATGCAGGTTTCTCAGTACACGGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCa                                                                                                   <  2:521662/150‑1 (MQ=255)

                                                      cTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAgc                                                                                                 <  2:1027010/150‑1 (MQ=255)

                                                      cTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAgc                                                                                                 >  1:142073/1‑150 (MQ=255)

                                                          cGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCAc                                                                                             <  1:1074681/150‑1 (MQ=255)

                                                              gCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACgcgc                                                                                         <  2:755620/150‑1 (MQ=255)

                                                                cGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTTTCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTg                                                                                       >  2:1446324/1‑150 (MQ=255)

                                                                cGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTg                                                                                       >  1:1054770/1‑150 (MQ=255)

                                                                     gagTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCAt                                                                                  <  2:1568198/150‑1 (MQ=255)

                                                                      agTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATa                                                                                 >  2:1546652/1‑150 (MQ=255)

                                                                        tACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGTTACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATaaa                                                                               <  2:1548190/150‑1 (MQ=255)

                                                                          ccAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACg                                                                             >  2:1568822/1‑150 (MQ=255)

                                                                           cAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCAGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGa                                                                            >  1:1153458/1‑150 (MQ=255)

                                                                            aGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGac                                                                           >  2:430894/1‑150 (MQ=255)

                                                                              aCCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGacac                                                                         >  1:1165767/1‑150 (MQ=255)

                                                                              aCCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGacac                                                                         >  1:282672/1‑150 (MQ=255)

                                                                                 tACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAAt                                                                      >  1:463774/1‑150 (MQ=255)

                                                                                      ggCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCg                                                                 <  1:1051955/150‑1 (MQ=255)

                                                                                       gcgcAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCgg                                                                <  2:202310/150‑1 (MQ=255)

                                                                                        cgcAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACAGCTATAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCTCGCGCTGCGCATAAACGACACAATGCCCGGt                                                               <  1:699774/150‑1 (MQ=255)

                                                                                        cgcAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGt                                                               >  2:177261/1‑150 (MQ=255)

                                                                                            ggTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAAt                                                           >  2:1393776/1‑150 (MQ=255)

                                                                                               aGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATgag                                                        >  1:953610/1‑150 (MQ=255)

                                                                                               aGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATgag                                                        >  1:1080298/1‑150 (MQ=255)

                                                                                               aGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATgag                                                        >  2:1243178/1‑150 (MQ=255)

                                                                                               aGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATgag                                                        <  2:1165767/150‑1 (MQ=255)

                                                                                               aGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTACAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATgag                                                        >  2:1369274/1‑150 (MQ=255)

                                                                                               aGATAAAACAGACGTTGATACTTACCGTTATCTGAAATTCAACCAGCTTTCTCAGTACACAGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATTCCCGGTGAATgag                                                        >  2:1590474/1‑150 (MQ=255)

                                                                                                  tAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGAtt                                                     >  1:827452/1‑150 (MQ=255)

                                                                                                  tAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGAtt                                                     >  2:612542/1‑150 (MQ=255)

                                                                                                     aaCAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGTGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTc‑cc                                                 <  2:344843/150‑1 (MQ=255)

                                                                                                                 tACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGAGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCAttttttt                                     >  2:974976/1‑150 (MQ=255)

                                                                                                                   cTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCAGGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTAt                                   <  1:343391/150‑1 (MQ=255)

                                                                                                                   cTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTAt                                   <  1:995652/150‑1 (MQ=255)

                                                                                                                     tACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTAttt                                 >  1:854971/1‑150 (MQ=255)

                                                                                                                      aCCGTTATCTGAACTACAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTc                                <  1:1369274/150‑1 (MQ=255)

                                                                                                                       ccGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCt                               <  1:430894/150‑1 (MQ=255)

                                                                                                                       ccGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGA‑TTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTa                              >  1:525859/1‑150 (MQ=255)

                                                                                                                           tATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAAc                           >  2:1424051/1‑150 (MQ=255)

                                                                                                                             tCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCa                         <  1:1607470/150‑1 (MQ=255)

                                                                                                                             tCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCa                         >  2:1416156/1‑150 (MQ=255)

                                                                                                                              cTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCAt                        >  2:1135667/1‑150 (MQ=255)

                                                                                                                                 aaCTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGc                     <  1:855407/150‑1 (MQ=255)

                                                                                                                                   cTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCg                   >  2:50539/1‑150 (MQ=255)

                                                                                                                                   cTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCg                   >  1:34778/1‑150 (MQ=255)

                                                                                                                                     cccatcggCTTTCTCGGTGTACTGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGtt                 >  2:843384/8‑150 (MQ=255)

                                                                                                                                      cAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTc                >  2:1410225/1‑150 (MQ=255)

                                                                                                                                       aaCCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTc‑c                                                  >  1:679391/1‑115 (MQ=255)

                                                                                                                                       aaCCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTc‑c                                                  <  2:679391/115‑1 (MQ=255)

                                                                                                                                       aaCCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTcc               >  2:1238008/1‑150 (MQ=255)

                                                                                                                                        atcggCTTTCTCGGTGTACTGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCg              >  1:1217453/5‑150 (MQ=255)

                                                                                                                                        atcggCTTTCTCGGTGTACTGAGAAAGCCGATGCGGTGAAATTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCg              >  2:832554/5‑150 (MQ=255)

                                                                                                                                         ccAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCgg                                                                <  1:936327/100‑1 (MQ=255)

                                                                                                                                         ccAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCgg                                                                >  2:936327/1‑100 (MQ=255)

                                                                                                                                          cAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCt            <  2:761066/150‑1 (MQ=255)

                                                                                                                                           ggCTTTCTCGGTGTACTGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCTCGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTg           >  1:86832/2‑150 (MQ=255)

                                                                                                                                           ggCTTTCTCGGTGTACTGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTg           >  1:69721/2‑150 (MQ=255)

                                                                                                                                           ggCTTTCTCGGTGTACTGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTg           >  1:502296/2‑150 (MQ=255)

                                                                                                                                           ggCTTTCTCGGTGTACTGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTg           >  2:1591432/2‑150 (MQ=255)

                                                                                                                                           ggCTTTCTCGGTGTACTGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTg           >  1:86830/2‑150 (MQ=255)

                                                                                                                                           aGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTg           <  1:1568822/150‑1 (MQ=255)

                                                                                                                                            tcTTTCTCAGTACACCGATAAAGCCGATGGGGTGATTTTCCAGACTGCGTTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCTGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGt          <  2:1153458/149‑1 (MQ=255)

                                                                                                                                              tttCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGttt        <  1:1078196/150‑1 (MQ=255)

                                                                                                                                              tttCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTATAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGCGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGttt        >  2:1444935/1‑150 (MQ=255)

                                                                                                                                                tctcAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGttttt      >  1:1546110/1‑150 (MQ=255)

                                                                                                                                                 ctcAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATTCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTtc     <  2:282672/150‑1 (MQ=255)

                                                                                                                                                 ctcAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTtc     >  1:635611/1‑150 (MQ=255)

                                                                                                                                                 ctcAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTtc     <  2:953610/150‑1 (MQ=255)

                                                                                                                                                  tcAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTtct    >  2:1123030/1‑150 (MQ=255)

                                                                                                                                                  tcAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGG‑GAATCTCATTCACCGGGCAttgtgt                                      <  1:1598010/116‑5 (MQ=37)

                                                                                                                                                   cAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCGGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTtctc   >  2:637988/1‑150 (MQ=255)

                                                                                                                                                   cAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTtctc   >  2:1372971/1‑150 (MQ=255)

AATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTGCGCATAAACGACACAATGCCCGGTGAATGAGATTC‑CCGGGCATTTTTTTATTTCTAAACCATCGCCGTTCCGCTGTTTTTCTCC  >  NC_000913/134018‑134315

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