BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,021,376	A→G	D580D (GAT→GAC)	b0960 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,021,376	0	A	G	99.5%	672.0 / ‑5.7	217	D580D (GAT→GAC)	b0960	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (96/120); total (96/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 222 total aligned reads displayed.

AAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTGGCCTC  >  NC_000913/1021286‑1021465

 aaCGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTc‑g                                                                                        <  1:126134/93‑1 (MQ=255)

     cgGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑gcggc                                                                                    <  1:2991204/93‑1 (MQ=255)

      ggCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑gcggcg                                                                                   <  1:6572172/93‑1 (MQ=255)

      ggCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑gcggcg                                                                                   <  1:2447183/93‑1 (MQ=255)

      ggCTTCGCGAATTTGCGGGGGAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGCGCGTC‑gcggcg                                                                                   >  2:2571076/1‑93 (MQ=255)

       gCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑Gcggcgg                                                                                  >  1:5900667/1‑93 (MQ=255)

         ttCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATGCCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGc                                                                                <  1:1940679/93‑1 (MQ=255)

         ttCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGc                                                                                >  1:5399994/1‑93 (MQ=255)

         ttCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGc                                                                                >  2:2523621/1‑93 (MQ=255)

           cgcgAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCaa                                                                              <  1:4340165/93‑1 (MQ=255)

           cgcgAATTTGCGGGGTAACGTTCAGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCaa                                                                              <  2:4273806/93‑1 (MQ=255)

            gcgAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAat                                                                             >  2:7517570/1‑93 (MQ=255)

             cgAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAata                                                                            <  2:7282112/93‑1 (MQ=255)

              gAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAatac                                                                           <  2:5159472/93‑1 (MQ=255)

                 tttGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAatacgat                                                                        >  2:6640653/1‑93 (MQ=255)

                 tttGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAatacgat                                                                        >  2:3314162/1‑93 (MQ=255)

                 tttGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAatacgat                                                                        >  2:3116150/1‑93 (MQ=255)

                   tGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAatacgatac                                                                      >  1:263692/1‑93 (MQ=255)

                   tGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAatacgatac                                                                      <  1:2930164/93‑1 (MQ=255)

                   tGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAatacgatac                                                                      >  1:233800/1‑93 (MQ=255)

                    gCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAatacgatacg                                                                     <  1:3642666/93‑1 (MQ=255)

                      ggggTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGcc                                                                   >  2:2283155/1‑93 (MQ=255)

                      ggggTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGcc                                                                   >  2:7651411/1‑93 (MQ=255)

                       gggTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCa                                                                  <  1:6692024/93‑1 (MQ=255)

                        ggTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAg                                                                 <  1:7486786/93‑1 (MQ=255)

                         gTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGa                                                                <  2:5854949/93‑1 (MQ=255)

                         gTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGa                                                                >  1:5757006/1‑93 (MQ=255)

                         gTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGa                                                                >  2:5891809/1‑93 (MQ=255)

                          tAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGAc                                                               <  2:3503047/93‑1 (MQ=255)

                            aCGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACgg                                                             <  2:1488577/93‑1 (MQ=255)

                            aCGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGATGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACgg                                                             >  2:1832497/1‑93 (MQ=255)

                             cGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGt                                                            <  1:1689644/93‑1 (MQ=255)

                             cGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGt                                                            >  2:6665460/1‑93 (MQ=255)

                              gTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGtt                                                           <  2:4481125/93‑1 (MQ=255)

                              gTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGtt                                                           >  2:5323402/1‑93 (MQ=255)

                                tCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTAt                                                         >  2:747901/1‑93 (MQ=255)

                                 cGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATc                                                        >  1:5048266/1‑93 (MQ=255)

                                  ggCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCa                                                       <  1:2885567/93‑1 (MQ=255)

                                  ggCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCa                                                       <  1:2684888/93‑1 (MQ=255)

                                  ggCTCGCTGCACATATTGGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCa                                                       <  1:642472/93‑1 (MQ=255)

                                       gCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGccc                                                  >  2:2430845/1‑93 (MQ=255)

                                        cTGGACATATTTGATCCCCCCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGcccc                                                 <  2:1585997/93‑1 (MQ=255)

                                        cTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGcccc                                                 >  1:7313222/1‑93 (MQ=255)

                                         tGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCt                                                >  1:7742020/1‑93 (MQ=255)

                                          ggACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTg                                               >  2:4839874/1‑93 (MQ=255)

                                           gACATATTTGATAGCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGa                                              <  2:4580338/93‑1 (MQ=255)

                                           gACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGa                                              <  1:6974657/93‑1 (MQ=255)

                                           gACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGa                                              <  1:2040196/93‑1 (MQ=255)

                                           gACATATGTGTTACCACCGACGCCAGCTCAGCATCATTGTTTTTTCAAGCAGAAGACGGC‑ATACGAGATGCCAGACGGTTATCACGCCCCTGa                                              <  1:2113835/93‑1 (MQ=255)

                                            aCATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGAt                                             <  2:55115/93‑1 (MQ=255)

                                            aCATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGAt                                             <  1:5661797/93‑1 (MQ=255)

                                             cATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATg                                            <  2:4201076/93‑1 (MQ=255)

                                             cATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATg                                            >  1:1992929/1‑93 (MQ=255)

                                                atTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTa                                         <  2:2126009/93‑1 (MQ=255)

                                                atTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTa                                         <  2:1113087/93‑1 (MQ=255)

                                                 tttGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTAt                                        >  1:7812590/1‑93 (MQ=255)

                                                 tttGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTAt                                        >  2:5688748/1‑93 (MQ=255)

                                                 tttGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTAt                                        <  1:1233200/93‑1 (MQ=255)

                                                      tACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTc                                   <  2:6737447/93‑1 (MQ=255)

                                                         caccGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAg                                >  2:2445201/1‑93 (MQ=255)

                                                         caccGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAg                                <  1:4649787/93‑1 (MQ=255)

                                                          accGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGt                               >  1:4418151/1‑93 (MQ=255)

                                                           ccGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTa                              >  1:116793/1‑93 (MQ=255)

                                                            cGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTGTCACGCCCCTGATGGTATTGCTCCAGTAt                             <  1:4241803/93‑1 (MQ=255)

                                                            cGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTAt                             >  2:4547172/1‑93 (MQ=255)

                                                            cGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTAt                             <  2:3631820/93‑1 (MQ=255)

                                                            cGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTAt                             <  1:6798103/93‑1 (MQ=255)

                                                            cGACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTAt                             >  2:4123689/1‑93 (MQ=255)

                                                             gACGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATg                            <  1:740994/93‑1 (MQ=255)

                                                              aCGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATgg                           >  2:7796412/1‑93 (MQ=255)

                                                               cGCCAGCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGc                          <  2:2664232/93‑1 (MQ=255)

                                                                    gCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTGTCACGCCCCTGATGGTATTGCTCCAGTATGGCATCga                     >  2:141286/1‑93 (MQ=255)

                                                                    gCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCga                     <  2:744205/93‑1 (MQ=255)

                                                                    gCTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCga                     <  2:5056928/93‑1 (MQ=255)

                                                                     cTCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCgag                    >  1:7560364/1‑93 (MQ=255)

                                                                      tCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCgaga                   >  2:3841460/1‑93 (MQ=255)

                                                                      tCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCgaga                   <  1:1537061/93‑1 (MQ=255)

                                                                      tCAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCgaga                   <  1:6244153/93‑1 (MQ=255)

                                                                       cAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGTCATCTAGAt                  >  2:3790094/1‑93 (MQ=255)

                                                                       cAGCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGAt                  <  2:7844736/93‑1 (MQ=255)

                                                                         gCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATAc                <  2:4153096/93‑1 (MQ=255)

                                                                         gCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATAc                <  2:3443725/93‑1 (MQ=255)

                                                                         gCATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATAc                <  1:3572474/93‑1 (MQ=255)

                                                                          cATCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATAcc               >  1:3804957/1‑93 (MQ=255)

                                                                            tCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGa             >  2:71502/1‑93 (MQ=255)

                                                                            tCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGa             >  2:4565697/1‑93 (MQ=255)

                                                                            tCACGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGATATCACGCCCCTGATGGTATTGATCCAATATGGCATCGACATACCGa             >  1:2830996/1‑93 (MQ=255)

                                                                              aCGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACa           >  2:6128904/1‑93 (MQ=255)

                                                                              aCGGTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACa           >  1:107106/1‑93 (MQ=255)

                                                                                ggTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGt         <  1:3670699/93‑1 (MQ=255)

                                                                                 gTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGtt        >  2:2874733/1‑93 (MQ=255)

                                                                                 gTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGtt        <  1:4448722/93‑1 (MQ=255)

                                                                                 gTTGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGtt        >  1:986408/1‑93 (MQ=255)

                                                                                  ttGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTg       >  1:5189035/1‑93 (MQ=255)

                                                                                  ttGTGTGCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTg       >  1:4181639/1‑93 (MQ=255)

                                                                                   tgtgtgCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTgg      <  2:7212160/93‑1 (MQ=255)

                                                                                     tgtgCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTGGcc    >  1:3285468/1‑93 (MQ=255)

                                                                                      gtgCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTGGCct   >  2:1751640/1‑93 (MQ=255)

                                                                                       tgCGTC‑GCGGCGGGCAATATGGTACGCCAGACGGGTATCACGCCCCTGATGGGACTGCTCCAGTATGGCATCGAGATACCGACAGTTGGCctc  >  2:6919884/1‑93 (MQ=255)

                                                                                       tgCGTC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTGGCctc  <  2:1038646/93‑1 (MQ=255)

AAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATC‑GCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTGGCCTC  >  NC_000913/1021286‑1021465

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