BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,561,521	C→T	S754N (AGC→AAC)	b1489 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,561,521	0	C	T	99.4%	554.1 / ‑5.2	161	S754N (AGC→AAC)	b1489	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (73/87); total (73/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-01

Only 100 of 171 total aligned reads displayed.

TATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATT  >  NC_000913/1561431‑1561610

tatCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGttt                                                                                         >  2:3069518/1‑93 (MQ=255)

tatCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGttt                                                                                         <  2:244602/93‑1 (MQ=255)

tatCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGttt                                                                                         >  2:426527/1‑93 (MQ=255)

tatCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGttt                                                                                         <  2:2069751/93‑1 (MQ=255)

 atCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGtttt                                                                                        <  1:1342884/93‑1 (MQ=255)

 atCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGtttt                                                                                        <  2:650539/93‑1 (MQ=255)

  tCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTg                                                                                       <  2:1745440/93‑1 (MQ=255)

      tGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCa                                                                                   >  2:4277180/1‑93 (MQ=255)

      tGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCa                                                                                   >  1:2015469/1‑93 (MQ=255)

        aaTAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAAt                                                                                 >  2:1950365/1‑93 (MQ=255)

         aTAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATg                                                                                >  2:967699/1‑93 (MQ=255)

         aTAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATg                                                                                <  2:2701421/93‑1 (MQ=255)

         aTAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATg                                                                                >  2:1725159/1‑93 (MQ=255)

           aaCGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCAAATGCt                                                                              >  1:1381012/1‑93 (MQ=255)

            aCGCGACATTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTg                                                                             <  1:750085/93‑1 (MQ=255)

             cgcgACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTgg                                                                            <  1:3713515/93‑1 (MQ=255)

             cgcgACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTgg                                                                            <  2:807443/93‑1 (MQ=255)

             cgcgACAGTGGATCTTGCGTAGCATCTAAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTgg                                                                            <  2:2603203/93‑1 (MQ=255)

               cgACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTa                                                                          <  2:165656/93‑1 (MQ=255)

                gACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTaa                                                                         >  1:1438348/1‑93 (MQ=255)

                gACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTaa                                                                         <  2:3237994/93‑1 (MQ=255)

                 gcAGTGGATCTTGCGTAGCATCGCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAAt                                                                        <  1:3677186/92‑1 (MQ=255)

                 aCAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAAt                                                                        >  1:2065960/1‑93 (MQ=255)

                 aCAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAAt                                                                        <  1:4453934/93‑1 (MQ=255)

                 aCAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAAt                                                                        <  1:4602948/93‑1 (MQ=255)

                    gTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGc                                                                     >  2:2722518/1‑93 (MQ=255)

                    gTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGGTAAATTGAGGTTTTGCCCAATGCTGGTAATGGc                                                                     >  1:4566066/1‑93 (MQ=255)

                        aTCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCa                                                                 >  1:1600121/1‑93 (MQ=255)

                         tCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCaa                                                                <  2:178164/93‑1 (MQ=255)

                         tCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCaa                                                                >  2:669911/1‑93 (MQ=255)

                         tCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCaa                                                                >  2:3606874/1‑93 (MQ=255)

                          cTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAg                                                               <  2:1989052/93‑1 (MQ=255)

                            tGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTa                                                             <  1:3808374/93‑1 (MQ=255)

                             gCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTaa                                                            >  1:4219163/1‑93 (MQ=255)

                             gCGTAGCATATCAAATTGCTCTTTGGTTTCGACGCCTTCAGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTaa                                                            <  1:3296017/93‑1 (MQ=255)

                              cGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAg                                                           <  2:1889802/93‑1 (MQ=255)

                              cGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAg                                                           <  2:4595642/93‑1 (MQ=255)

                                tAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGc                                                         <  1:635208/93‑1 (MQ=255)

                                 aGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCa                                                        >  2:4537917/1‑93 (MQ=255)

                                 aGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCa                                                        >  2:3697789/1‑93 (MQ=255)

                                  gCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCaa                                                       >  1:3237776/1‑93 (MQ=255)

                                  gCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCaa                                                       <  1:2595280/93‑1 (MQ=255)

                                  gCATCTCAAATTGCTCTTTGGTTTCAACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCaa                                                       <  1:33585/93‑1 (MQ=255)

                                    aTCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAgg                                                     >  2:584114/1‑93 (MQ=255)

                                     tctcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGa                                                    >  2:2080963/1‑93 (MQ=255)

                                       tcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATg                                                  >  1:2235244/1‑93 (MQ=255)

                                       tcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATg                                                  >  1:378049/1‑93 (MQ=255)

                                        cAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGc                                                 <  1:3394847/93‑1 (MQ=255)

                                          aaTTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGt                                               <  1:2846309/93‑1 (MQ=255)

                                             tGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGtttt                                            <  2:2509625/93‑1 (MQ=255)

                                              gCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGttttt                                           <  2:2840930/93‑1 (MQ=255)

                                                tctTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCg                                         >  2:355469/1‑93 (MQ=255)

                                                tctTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCg                                         >  2:544988/1‑93 (MQ=255)

                                                    tGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCa                                     <  2:2208766/93‑1 (MQ=255)

                                                    tGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCa                                     <  2:1327595/93‑1 (MQ=255)

                                                     ggTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAg                                    >  2:1576943/1‑93 (MQ=255)

                                                     ggTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAg                                    <  1:4033302/93‑1 (MQ=255)

                                                      gTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGa                                   >  1:2501488/1‑93 (MQ=255)

                                                       tttCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGAc                                  <  2:309823/93‑1 (MQ=255)

                                                       tttCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGAc                                  <  2:4410801/93‑1 (MQ=255)

                                                         tCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACaa                                <  2:4521728/93‑1 (MQ=255)

                                                          cGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAAc                               <  1:2968601/93‑1 (MQ=255)

                                                          cGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAAc                               <  1:1242558/93‑1 (MQ=255)

                                                            aCGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGa                             >  2:4303522/1‑93 (MQ=255)

                                                            aCGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGa                             >  1:246207/1‑93 (MQ=255)

                                                            aCGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGa                             >  2:4193030/1‑93 (MQ=255)

                                                             cGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGAt                            <  1:455475/93‑1 (MQ=255)

                                                              gCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACCATc                           <  2:1746792/93‑1 (MQ=255)

                                                              gCCTTCCGCCACGACGGGTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATc                           >  1:2803240/1‑93 (MQ=255)

                                                               ccTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCg                          <  1:1576943/93‑1 (MQ=255)

                                                                 ttCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGAc                        >  1:485860/1‑93 (MQ=255)

                                                                 ttCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGAc                        >  1:1845701/1‑93 (MQ=255)

                                                                     gCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACaaaa                    >  2:4689883/1‑93 (MQ=255)

                                                                      ccACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAAc                   <  1:1119730/93‑1 (MQ=255)

                                                                       cACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACt                  >  2:1092699/1‑93 (MQ=255)

                                                                        ccgacgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtt                 <  2:1465146/92‑1 (MQ=255)

                                                                          gacgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtttt               <  2:4585599/93‑1 (MQ=255)

                                                                          gacgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtttt               <  1:2571402/93‑1 (MQ=255)

                                                                          gacgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtttt               >  1:2584301/1‑93 (MQ=255)

                                                                           acgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTg              <  1:2061874/93‑1 (MQ=255)

                                                                            cgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGt             >  1:569173/1‑93 (MQ=255)

                                                                             ggTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGGTCGACAAAACTTTTGTc            <  2:389489/93‑1 (MQ=255)

                                                                              gTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCa           >  2:914590/1‑93 (MQ=255)

                                                                              attAAATAGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCa           >  2:2883064/2‑93 (MQ=255)

                                                                               ttAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCaa          >  1:2563271/1‑93 (MQ=255)

                                                                               ttAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCaa          >  1:3382715/1‑93 (MQ=255)

                                                                                 aaaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAtt        <  2:857977/93‑1 (MQ=255)

                                                                                 aaaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAtt        <  2:527758/93‑1 (MQ=255)

                                                                                 aaaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAtt        >  2:2858145/1‑93 (MQ=255)

                                                                                  aaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAttt       <  1:3112/93‑1 (MQ=255)

                                                                                  aaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAttt       <  2:1413561/93‑1 (MQ=255)

                                                                                  aaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAttt       >  2:4493549/1‑93 (MQ=255)

                                                                                   aTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAtttt      <  2:3347837/93‑1 (MQ=255)

                                                                                   aTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAtttt      >  1:1126208/1‑93 (MQ=255)

                                                                                   aTTGAGGTGTTGCCCAATGGTGGTAATGGGTTCAAGTAAGGCAAGGATGCCTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAtttt      <  1:3168339/93‑1 (MQ=255)

                                                                                     tGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGa    <  1:4425588/93‑1 (MQ=255)

                                                                                       aGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGAtt  <  1:3965298/93‑1 (MQ=255)

                                                                                       aGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGAtt  <  2:2830047/93‑1 (MQ=255)

                                                                                       aGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGAtt  <  1:4032268/93‑1 (MQ=255)

                                                                                       aGGTTTTGACCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGAtt  >  2:4147658/1‑93 (MQ=255)

TATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATT  >  NC_000913/1561431‑1561610

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