BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	54,157	C→T	P182P (CCG→CCA)	surA ←	survival protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	54,157	0	C	T	99.7%	1061.0 / ‑6.0	298	P182P (CCG→CCA)	surA	survival protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (151/146); total (151/147)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 301 total aligned reads displayed.

GAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCAG  >  NC_000913/54067‑54246

gAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTgg                                                                                         <  2:751077/93‑1 (MQ=255)

 aaaTCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTggg                                                                                        >  1:407281/1‑93 (MQ=255)

  aaTCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGa                                                                                       <  1:7329656/93‑1 (MQ=255)

   aTCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGAt                                                                                      <  1:581/93‑1 (MQ=255)

    tCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATc                                                                                     >  2:1758745/1‑93 (MQ=255)

     cAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCa                                                                                    <  2:3259356/93‑1 (MQ=255)

      aGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAg                                                                                   <  1:1257177/93‑1 (MQ=255)

      aGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAg                                                                                   >  1:6000770/1‑93 (MQ=255)

      aGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGACTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAg                                                                                   <  1:4852547/93‑1 (MQ=255)

        cgcCGTTACCCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGa                                                                                 <  1:335791/93‑1 (MQ=255)

          ccGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGAtg                                                                               >  1:3595871/1‑93 (MQ=255)

            gTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGAtgtg                                                                             <  1:2297809/93‑1 (MQ=255)

             ttACGCGCCTGATCGAGAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCAGGTTTTCCGGCAGTGGGATCAGGATGTgg                                                                            >  2:5216332/1‑93 (MQ=255)

             ttACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTgg                                                                            >  2:6972444/1‑93 (MQ=255)

             ttACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTgg                                                                            <  2:566776/93‑1 (MQ=255)

              tACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGc                                                                           <  1:998130/93‑1 (MQ=255)

               aCGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCt                                                                          >  1:5194121/1‑93 (MQ=255)

                 gcgcCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCa                                                                        <  1:2395098/93‑1 (MQ=255)

                  cgcCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAg                                                                       >  2:5087243/1‑93 (MQ=255)

                    ccTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGt                                                                     <  1:3102358/93‑1 (MQ=255)

                     cTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGtt                                                                    <  1:3051165/93‑1 (MQ=255)

                     cTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGtt                                                                    <  2:530871/93‑1 (MQ=255)

                     cTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGtt                                                                    >  1:6191021/1‑93 (MQ=255)

                     cTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGtt                                                                    >  1:7883149/1‑93 (MQ=255)

                      tGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTc                                                                   >  2:1272450/1‑93 (MQ=255)

                      tGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTc                                                                   <  2:4837012/93‑1 (MQ=255)

                      tGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTc                                                                   <  1:6808220/93‑1 (MQ=255)

                       gATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCa                                                                  <  2:389633/93‑1 (MQ=255)

                        aTCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAg                                                                 >  1:95437/1‑93 (MQ=255)

                          cGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCTCCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCt                                                               >  2:6012878/1‑93 (MQ=255)

                          cGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCt                                                               <  2:3956117/93‑1 (MQ=255)

                           gACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTc                                                              <  1:6686544/93‑1 (MQ=255)

                           gACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTc                                                              >  1:7262821/1‑93 (MQ=255)

                              aaTGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAgg                                                           >  2:3330306/1‑92 (MQ=255)

                                tGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGc                                                         <  1:818559/93‑1 (MQ=255)

                                tGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGc                                                         <  1:4736548/93‑1 (MQ=255)

                                tGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGc                                                         <  1:774616/93‑1 (MQ=255)

                                tGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGc                                                         >  1:756396/1‑93 (MQ=255)

                                 ggCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCt                                                        <  2:3785476/93‑1 (MQ=255)

                                   cgcgcgcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTgg                                                      <  1:6308399/93‑1 (MQ=255)

                                    gcgcgcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGc                                                     <  1:1315094/93‑1 (MQ=255)

                                      gcgcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgt                                                   <  1:6526237/93‑1 (MQ=255)

                                      gcgcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGGTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGTcgt                                                   >  2:2797182/1‑93 (MQ=255)

                                        gcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgtcg                                                 >  2:2003637/1‑93 (MQ=255)

                                        gcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgtcg                                                 >  1:577053/1‑93 (MQ=255)

                                         ccTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgtcgt                                                >  2:6647413/1‑93 (MQ=255)

                                          cTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtt                                               >  1:1048343/1‑93 (MQ=255)

                                          cTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtt                                               >  2:6689318/1‑93 (MQ=255)

                                           tGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGttt                                              >  2:4646790/1‑93 (MQ=255)

                                           tGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGttt                                              >  2:6888540/1‑93 (MQ=255)

                                           tGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGttt                                              >  1:3879599/1‑93 (MQ=255)

                                           tGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGttt                                              <  1:4957208/93‑1 (MQ=255)

                                            ggCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtttt                                             <  2:579273/93‑1 (MQ=255)

                                            ggCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtttt                                             <  2:4887857/93‑1 (MQ=255)

                                            ggCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtttt                                             >  2:7773238/1‑93 (MQ=255)

                                             gCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTg                                            >  2:4777476/1‑93 (MQ=255)

                                              cTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTgg                                           <  2:1524701/93‑1 (MQ=255)

                                                ttCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGtt                                         <  1:4638546/93‑1 (MQ=255)

                                                 tCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTa                                        <  2:6987054/93‑1 (MQ=255)

                                                  ccGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGGTCAGCTCAGTGCTGGCGTCGTTTTGGTTAc                                       >  2:2973431/1‑93 (MQ=255)

                                                   cGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTAcc                                      >  1:7317573/1‑93 (MQ=255)

                                                    gCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTAccc                                     >  2:842501/1‑93 (MQ=255)

                                                      ttCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCAc                                   >  1:2869674/1‑93 (MQ=255)

                                                      ttCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCAc                                   <  2:7342126/93‑1 (MQ=255)

                                                       tCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCGcc                                  <  1:1878606/93‑1 (MQ=255)

                                                       tCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCAcc                                  <  2:2396228/93‑1 (MQ=255)

                                                        cGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACct                                 <  1:7117424/93‑1 (MQ=255)

                                                         gTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctg                                >  1:2970460/1‑93 (MQ=255)

                                                         gTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctg                                <  2:3252213/93‑1 (MQ=255)

                                                          gtCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctgc                               <  1:3757047/92‑1 (MQ=255)

                                                           tCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctgct                              >  2:5494026/1‑93 (MQ=255)

                                                           tCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctgct                              >  1:3911038/1‑93 (MQ=255)

                                                             aCCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCtgctgc                            <  2:2459611/93‑1 (MQ=255)

                                                               cTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTgcgc                          >  2:170930/1‑93 (MQ=255)

                                                                tGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGcc                         <  2:4975626/93‑1 (MQ=255)

                                                                tGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGcc                         >  1:4713312/1‑93 (MQ=255)

                                                                tGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGcc                         >  1:4254492/1‑93 (MQ=255)

                                                                 gATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCa                        <  1:6176440/93‑1 (MQ=255)

                                                                 gATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCa                        >  1:3225884/1‑93 (MQ=255)

                                                                  aTCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAg                       <  1:6574933/93‑1 (MQ=255)

                                                                  aTCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAg                       <  1:7684553/93‑1 (MQ=255)

                                                                   tCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAgg                      >  2:4536525/1‑93 (MQ=255)

                                                                   tCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAgg                      <  1:1656056/93‑1 (MQ=255)

                                                                     agagGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGa                    >  2:4040209/1‑93 (MQ=255)

                                                                      gagGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAt                   >  1:2782304/1‑93 (MQ=255)

                                                                      gagGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAt                   <  2:6499773/93‑1 (MQ=255)

                                                                      gagGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAt                   >  2:2082675/1‑93 (MQ=255)

                                                                       gggTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAtt                  <  1:1184893/92‑1 (MQ=255)

                                                                       agGTCGGGTTTTCCGGCAGTGGGCTCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAtt                  <  2:4301546/93‑1 (MQ=255)

                                                                       agGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAtt                  >  1:2125692/1‑93 (MQ=255)

                                                                        ggTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTc                 >  2:2564260/1‑93 (MQ=255)

                                                                         gTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCg                >  1:6408346/1‑93 (MQ=255)

                                                                         gTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCg                >  1:2748942/1‑93 (MQ=255)

                                                                          tCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGa               >  2:742906/1‑93 (MQ=255)

                                                                                tttCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCt         <  1:5537330/93‑1 (MQ=255)

                                                                                   ccGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCg      >  2:6821502/1‑93 (MQ=255)

                                                                                      gCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGcg   >  2:2172860/1‑92 (MQ=255)

                                                                                      gCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCa   >  1:2597721/1‑93 (MQ=255)

                                                                                      gCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCa   >  2:4236799/1‑93 (MQ=255)

                                                                                      gCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCa   <  2:1842225/93‑1 (MQ=255)

GAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCAG  >  NC_000913/54067‑54246

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 23 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`