BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,701,551	T→C	Q178R (CAG→CGG)	rnc ←	RNase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,701,551	0	T	C	99.5%	669.5 / ‑4.8	192	Q178R (CAG→CGG)	rnc	RNase III
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (76/115); total (77/115)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

Only 100 of 199 total aligned reads displayed.

CGGCTTGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCTGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTATCGCCTGG  >  NC_000913/2701461‑2701639

cgGCTTGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCgg                                                                                        >  2:4272710/1‑93 (MQ=255)

  gCTTGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGAc                                                                                      <  1:4492310/93‑1 (MQ=255)

    ttGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTa                                                                                    >  1:2438813/1‑93 (MQ=255)

     tGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTAc                                                                                   >  2:4407457/1‑93 (MQ=255)

       aaCCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCa                                                                                 >  1:6016549/1‑93 (MQ=255)

       aaCCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCa                                                                                 <  1:1763516/93‑1 (MQ=255)

        aCCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAg                                                                                >  1:1011253/1‑93 (MQ=255)

          cTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGAt                                                                              >  1:2227021/1‑93 (MQ=255)

             tgCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAg                                                                           <  2:5750805/93‑1 (MQ=255)

              gCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGt                                                                          <  1:1034079/93‑1 (MQ=255)

               ccAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTc                                                                         <  1:710030/93‑1 (MQ=255)

               ccAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTc                                                                         <  1:3991280/93‑1 (MQ=255)

               ccAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTc                                                                         <  1:3087838/93‑1 (MQ=255)

               ccAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTc                                                                         <  1:2065575/93‑1 (MQ=255)

               ccAACCACCGGTTCACTAAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTc                                                                         >  2:2971984/1‑93 (MQ=255)

                cAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTcg                                                                        >  1:5542196/1‑93 (MQ=255)

                  accaccGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTcggc                                                                      >  2:6464653/1‑93 (MQ=255)

                   ccaccGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTcggca                                                                     <  2:314299/93‑1 (MQ=255)

                   ccaccGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTcggca                                                                     <  2:2172516/93‑1 (MQ=255)

                   ccaccGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTcggca                                                                     >  1:4342079/1‑93 (MQ=255)

                        ggTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTcggcagcggc                                                                <  1:2948246/93‑1 (MQ=255)

                         gTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTcggcagcggca                                                               >  2:247738/1‑93 (MQ=255)

                          ttCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTcggcagcggcag                                                              >  2:481197/1‑93 (MQ=255)

                           tCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGa                                                             <  2:5631890/93‑1 (MQ=255)

                            cACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGAt                                                            >  1:6077984/1‑93 (MQ=255)

                            cACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGAt                                                            <  1:917908/93‑1 (MQ=255)

                            cACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGAt                                                            <  2:2059351/93‑1 (MQ=255)

                            cACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGAt                                                            <  1:5034051/93‑1 (MQ=255)

                             aCTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATg                                                           <  2:6036568/93‑1 (MQ=255)

                                cAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCg                                                        <  1:4640927/93‑1 (MQ=255)

                                  ggCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGGGAc                                                      >  2:1724425/1‑93 (MQ=255)

                                    ccGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGAccc                                                    >  1:3759278/1‑93 (MQ=255)

                                    ccGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGAccc                                                    >  2:2735369/1‑93 (MQ=255)

                                    ccGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGAccc                                                    <  1:6226783/93‑1 (MQ=255)

                                     cGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCt                                                   <  2:4736996/93‑1 (MQ=255)

                                      gCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTg                                                  >  1:5233522/1‑93 (MQ=255)

                                      gCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTg                                                  >  2:1397007/1‑93 (MQ=255)

                                       cTGACGTGCCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGc                                                 <  2:5012754/93‑1 (MQ=255)

                                       cTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGc                                                 <  2:143855/93‑1 (MQ=255)

                                        tGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCa                                                >  2:1571090/1‑93 (MQ=255)

                                          aCCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCaaa                                              >  1:2635979/1‑93 (MQ=255)

                                           ccTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAat                                             >  2:828528/1‑93 (MQ=255)

                                            cTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAata                                            <  1:228666/93‑1 (MQ=255)

                                            cTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAata                                            <  1:5335789/93‑1 (MQ=255)

                                             tGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAatat                                           >  2:613868/1‑93 (MQ=255)

                                             tGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAatat                                           >  2:3694967/1‑93 (MQ=255)

                                             tGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAatat                                           <  2:1231657/93‑1 (MQ=255)

                                              ggCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATAtt                                          <  1:1095819/93‑1 (MQ=255)

                                              ggCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATAtt                                          >  1:2363594/1‑93 (MQ=255)

                                                 aGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCtt                                       <  2:3752521/93‑1 (MQ=255)

                                                  gTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTg                                      >  1:3859336/1‑93 (MQ=255)

                                                   gggATAGTAAATTCCGGAGCGTGCGCTTCGCCACGGGCCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGc                                     <  1:1128194/92‑1 (MQ=255)

                                                    ggATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCa                                    <  2:5806183/93‑1 (MQ=255)

                                                     gATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCaa                                   <  2:3810346/93‑1 (MQ=255)

                                                      aTAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAg                                  <  1:6465565/93‑1 (MQ=255)

                                                       tAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAgc                                 >  2:1303492/1‑93 (MQ=255)

                                                           aaaTTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGt                             <  2:1527402/93‑1 (MQ=255)

                                                           aaaTTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGt                             <  2:2566434/93‑1 (MQ=255)

                                                            aaTTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGtt                            <  2:4540395/93‑1 (MQ=255)

                                                            aaTTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGtt                            <  1:2490428/93‑1 (MQ=255)

                                                            aaTTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGtt                            <  1:6463935/93‑1 (MQ=255)

                                                             aTTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGttt                           >  2:594087/1‑93 (MQ=255)

                                                             aTTCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGttt                           <  1:1156956/93‑1 (MQ=255)

                                                              ttCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGtttt                          <  2:735249/93‑1 (MQ=255)

                                                               tCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGCCCCTGCAAATATTCTTGCAAGCGCGTTTTc                         <  2:3342876/93‑1 (MQ=255)

                                                               tCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTc                         <  2:3689579/93‑1 (MQ=255)

                                                               tCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTc                         >  1:3690842/1‑93 (MQ=255)

                                                               tCCTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTc                         <  1:5282694/93‑1 (MQ=255)

                                                                ccTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCg                        <  1:3164292/93‑1 (MQ=255)

                                                                ccTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCg                        <  1:2977100/93‑1 (MQ=255)

                                                                 cTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCgg                       <  2:1014494/93‑1 (MQ=255)

                                                                 cTGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCgg                       <  1:1315240/93‑1 (MQ=255)

                                                                  tGATCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGa                      <  2:6059816/93‑1 (MQ=255)

                                                                     tCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCt                   <  1:1741068/93‑1 (MQ=255)

                                                                     tCGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGAAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCt                   >  2:5333215/1‑93 (MQ=255)

                                                                      cGTGCGCTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCtt                  <  1:5149801/93‑1 (MQ=255)

                                                                         gcgcTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCttttt               <  2:5586553/93‑1 (MQ=255)

                                                                          cgcTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTg              >  2:4854671/1‑93 (MQ=255)

                                                                           gcTTCGCCACGTACCTGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGt             >  2:5097315/1‑93 (MQ=255)

                                                                           gcTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGt             >  2:199626/1‑93 (MQ=255)

                                                                            cTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGtt            <  1:3570999/93‑1 (MQ=255)

                                                                            cTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGtt            >  2:782449/1‑93 (MQ=255)

                                                                            cTTCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGtt            >  2:362074/1‑93 (MQ=255)

                                                                             ttCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGttt           <  1:5889843/93‑1 (MQ=255)

                                                                             ttCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGttt           <  2:5450719/93‑1 (MQ=255)

                                                                             ttCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGAACCTGCAAATATTCTTGCAAGCGCGGTTTCGGATCTTTTTGttt           >  1:4151079/1‑93 (MQ=255)

                                                                              tCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTa          <  1:639749/93‑1 (MQ=255)

                                                                              tCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTa          <  2:1035342/93‑1 (MQ=255)

                                                                              tCGCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTa          >  1:5855481/1‑93 (MQ=255)

                                                                               cGCCACGTGCCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTAt         <  1:77277/93‑1 (MQ=255)

                                                                                gCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTATc        <  2:4696955/93‑1 (MQ=255)

                                                                                gCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTATc        >  2:4265068/1‑93 (MQ=255)

                                                                                gCCACGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTATc        >  2:3985065/1‑93 (MQ=255)

                                                                                   aCGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTATCGcc     <  1:1516497/93‑1 (MQ=255)

                                                                                   aCGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTATCGcc     <  2:4200590/93‑1 (MQ=255)

                                                                                   aCGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTATCGcc     >  2:2365698/1‑93 (MQ=255)

                                                                                    cGTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTATCGCCt    <  1:5757476/93‑1 (MQ=255)

                                                                                     gTACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTATCGCCTg   >  2:1088574/1‑93 (MQ=255)

                                                                                      tACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTATCGCCTgg  >  2:5190514/1‑93 (MQ=255)

                                                                                      tACCCGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTATCGCCTgg  >  1:744687/1‑93 (MQ=255)

CGGCTTGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGGATAGTAAATTCCTGATCGTGCGCTTCGCCACGTACCTGGACTACCAGATAAGTCGGCAGCGGCAGATGGCGACCCTGCAAATATTCTTGCAAGCGCGTTTTCGGATCTTTTTGTTTATCGCCTGG  >  NC_000913/2701461‑2701639

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