BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,757,874	T→C	H225R (CAT→CGT)	b1680 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,757,874	0	T	C	99.0%	649.6 / ‑2.9	202	H225R (CAT→CGT)	b1680	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (120/80); total (120/82)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

Only 100 of 209 total aligned reads displayed.

TCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTT  >  NC_000913/1757782‑1757962

tCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTAt                                                                                          >  1:1082067/1‑93 (MQ=255)

  aTAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGc                                                                                        <  1:8227627/93‑1 (MQ=255)

     gAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCTAATTCCGGTGGGGCCATACAGTTTACGCCCg                                                                                     >  1:2733735/1‑93 (MQ=255)

     gAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCg                                                                                     >  2:4460907/1‑93 (MQ=255)

     gAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAAGTCCGGTGGGGCCATACAGTTTATGCCCg                                                                                     <  1:6809/93‑1 (MQ=255)

     gAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGACAATTCCGGTGGGGCCATACAGTTTACGCCCg                                                                                     >  2:3625921/1‑93 (MQ=255)

       aCCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGga                                                                                   <  2:5886165/93‑1 (MQ=255)

       aCCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGga                                                                                   <  1:2127261/93‑1 (MQ=255)

        ccGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGgag                                                                                  <  1:5141452/93‑1 (MQ=255)

         cGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGgaga                                                                                 <  2:5190885/93‑1 (MQ=255)

          gCCCCCTTCCCACGGCGGCATCTCCTGCCACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGgaggc                                                                                >  1:5815618/1‑91 (MQ=255)

            ccccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGcaca                                                                              >  2:6404049/1‑93 (MQ=255)

            ccccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAaca                                                                              >  1:3291128/1‑93 (MQ=255)

             cccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGGCATACAGTTTACGCCCGGAGAacac                                                                             >  2:70534/1‑93 (MQ=255)

             cccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAacac                                                                             >  1:9546545/1‑93 (MQ=255)

             cccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAacac                                                                             >  2:9447999/1‑93 (MQ=255)

             cccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAacac                                                                             <  1:7527337/93‑1 (MQ=255)

               cTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGt                                                                           <  2:1132909/93‑1 (MQ=255)

               cTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGt                                                                           >  2:5480045/1‑93 (MQ=255)

                 tCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTaa                                                                         >  2:4147969/1‑93 (MQ=255)

                     ccggcggCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGt                                                                     <  1:2898521/92‑1 (MQ=255)

                     aCGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGt                                                                     >  1:9122328/1‑93 (MQ=255)

                      cggcggCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTc                                                                    <  2:1098956/93‑1 (MQ=255)

                        gcggcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGc                                                                  >  1:1212159/1‑93 (MQ=255)

                         cggcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCa                                                                 <  1:148627/93‑1 (MQ=255)

                          ggcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCaa                                                                >  1:7510190/1‑93 (MQ=255)

                          cgcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCaa                                                                <  2:3384653/92‑1 (MQ=255)

                            cATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATc                                                              >  2:7748875/1‑93 (MQ=255)

                              tctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCa                                                            >  1:1418472/1‑93 (MQ=255)

                              tctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCa                                                            >  2:1078591/1‑93 (MQ=255)

                              tctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCa                                                            >  2:2875075/1‑93 (MQ=255)

                              tctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCa                                                            <  2:9147712/93‑1 (MQ=255)

                              tctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCa                                                            >  2:191312/1‑93 (MQ=255)

                                   tGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTCAAAGTCGCAATCCAGCGcc                                                       >  2:1627434/1‑93 (MQ=255)

                                   tGCAACAAGGCTTCTGTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGcc                                                       <  1:9766322/93‑1 (MQ=255)

                                    gCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCt                                                      >  2:8776495/1‑93 (MQ=255)

                                     caacaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTg                                                     <  2:969329/93‑1 (MQ=255)

                                      aacaaGTCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGa                                                    <  2:4230747/93‑1 (MQ=255)

                                      aacaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGa                                                    >  2:303735/1‑93 (MQ=255)

                                       acaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGaa                                                   >  1:2416900/1‑93 (MQ=255)

                                       acaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGaa                                                   >  1:8036973/1‑93 (MQ=255)

                                       acaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGaa                                                   >  1:8192811/1‑93 (MQ=255)

                                         aaGGCTTCTTTCACATAAAGAATGCCCATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACa                                                 <  1:8266578/93‑1 (MQ=255)

                                          aGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACAt                                                >  1:2602444/1‑93 (MQ=255)

                                           ggCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATc                                               >  2:1195335/1‑93 (MQ=255)

                                            gCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATcc                                              >  2:9880229/1‑93 (MQ=255)

                                            gCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATcc                                              <  1:8797292/93‑1 (MQ=255)

                                            gCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATcc                                              <  1:9463159/93‑1 (MQ=255)

                                             cttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCa                                             <  1:230125/93‑1 (MQ=255)

                                              ttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAc                                            <  2:9254430/93‑1 (MQ=255)

                                              ttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAc                                            >  2:9050118/1‑93 (MQ=255)

                                              ttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAc                                            >  2:10019690/1‑93 (MQ=255)

                                              ttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAc                                            <  2:595743/93‑1 (MQ=255)

                                              ttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAc                                            >  1:5352945/1‑93 (MQ=255)

                                              ttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAc                                            >  2:2685413/1‑93 (MQ=255)

                                               tcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCAAGCGCCTGAACATCCAcg                                           >  1:2191658/1‑92 (MQ=255)

                                                cttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATTCACCg                                          >  2:7270398/1‑93 (MQ=255)

                                                 tttCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCgg                                         >  1:3428937/1‑93 (MQ=255)

                                                  ttCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGga                                        >  1:1302530/1‑93 (MQ=255)

                                                   tCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGgat                                       >  1:2644826/1‑93 (MQ=255)

                                                   tCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGgat                                       <  1:3068232/93‑1 (MQ=255)

                                                    cacaTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGgatg                                      <  2:8283816/93‑1 (MQ=255)

                                                     acaTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGgatga                                     <  1:7251400/93‑1 (MQ=255)

                                                     acaTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGgatga                                     <  1:3843308/93‑1 (MQ=255)

                                                           aGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGACGATGCATc                               <  1:6846291/93‑1 (MQ=255)

                                                            gAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCa                              >  2:2012622/1‑93 (MQ=255)

                                                            gAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCa                              >  2:8479914/1‑93 (MQ=255)

                                                             aaTGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCAc                             <  1:4912370/93‑1 (MQ=255)

                                                                   aaTTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTg                       >  2:6182569/1‑93 (MQ=255)

                                                                   aaTTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTg                       >  1:7849183/1‑93 (MQ=255)

                                                                   aaTTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTg                       >  2:6626008/1‑93 (MQ=255)

                                                                     ttCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAg                     >  2:8645682/1‑93 (MQ=255)

                                                                     ttCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAg                     >  1:643031/1‑93 (MQ=255)

                                                                      tCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAgc                    >  1:6931076/1‑93 (MQ=255)

                                                                        cGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAgcgc                  >  1:7018877/1‑93 (MQ=255)

                                                                         ggTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGcc                 >  2:7608133/1‑93 (MQ=255)

                                                                          gTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCa                <  2:7354320/93‑1 (MQ=255)

                                                                          gTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCa                >  2:5678784/1‑93 (MQ=255)

                                                                           tGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCAt               <  2:1684870/93‑1 (MQ=255)

                                                                           tGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCAt               <  2:240672/93‑1 (MQ=255)

                                                                            ggggCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATc              >  1:9149606/1‑93 (MQ=255)

                                                                             gggCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATcc             >  2:3182409/1‑93 (MQ=255)

                                                                             gggCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATcc             >  2:7683911/1‑93 (MQ=255)

                                                                               gCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCGCCGCCTGAGCGCCATccac           <  1:217458/93‑1 (MQ=255)

                                                                               gCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCAGCGCCTGAGCGCCATccac           >  1:6739292/1‑93 (MQ=255)

                                                                                ccATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccacc          <  2:6655956/93‑1 (MQ=255)

                                                                                 cATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccacca         >  2:7126997/1‑93 (MQ=255)

                                                                                 cATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccacca         <  1:9884576/93‑1 (MQ=255)

                                                                                 cATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccacca         <  1:191312/93‑1 (MQ=255)

                                                                                 cATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccacca         >  1:1655271/1‑93 (MQ=255)

                                                                                  aTACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAg        >  1:7249693/1‑93 (MQ=255)

                                                                                  aTACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAg        >  1:9485427/1‑93 (MQ=255)

                                                                                   tACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGc       >  2:7075946/1‑93 (MQ=255)

                                                                                    aCAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCa      <  2:6927954/93‑1 (MQ=255)

                                                                                      aGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACt    <  1:4269183/93‑1 (MQ=255)

                                                                                      aGTTTACGCCCGGAGAACACGTAAAAGTCGCAATACAGCGCCTGCACATCCACCGGCTAATGCATCACCGCCTGAGAGCCGCCCACCAACACt    >  1:5690062/1‑93 (MQ=255)

                                                                                       gTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACtt   >  2:4389696/1‑93 (MQ=255)

                                                                                       gTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACtt   >  2:9679660/1‑93 (MQ=255)

                                                                                        tttACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACttt  >  1:2383712/1‑93 (MQ=255)

                                                                                        tttACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACttt  >  1:4106708/1‑93 (MQ=255)

TCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTT  >  NC_000913/1757782‑1757962

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