BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,533,355	G→A	A176A (GCC→GCT)	yddE ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,533,355	0	G	A	99.5%	798.9 / NA	212	A176A (GCC→GCT)	yddE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (102/109); minor base C (1/0); total (103/109)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

GTTTTTGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGGGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAATATC  >  NC_000913/1533265‑1533444

  ttttGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGt                                                                                       <  1:4071438/93‑1 (MQ=255)

   tttGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTc                                                                                      <  2:911276/93‑1 (MQ=255)

   tttGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTc                                                                                      <  1:1570188/93‑1 (MQ=255)

   tttGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTc                                                                                      <  1:4290603/93‑1 (MQ=255)

     tGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGa                                                                                    >  2:1289465/1‑93 (MQ=255)

      gCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGat                                                                                   <  2:330664/93‑1 (MQ=255)

       ccTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGata                                                                                  <  1:1026510/93‑1 (MQ=255)

          gggCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATcc                                                                               <  2:2768948/93‑1 (MQ=255)

          gggCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATcc                                                                               <  1:5542931/93‑1 (MQ=255)

          gggCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATcc                                                                               >  1:1424833/1‑93 (MQ=255)

            gCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCAc                                                                             <  2:2912750/93‑1 (MQ=255)

            gCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGACATCCAc                                                                             <  2:4256529/93‑1 (MQ=255)

               gATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTc                                                                          >  2:4157629/1‑93 (MQ=255)

               gATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTc                                                                          >  1:5109949/1‑93 (MQ=255)

               gATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTc                                                                          >  2:5875418/1‑93 (MQ=255)

                aTCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTcc                                                                         >  1:243565/1‑93 (MQ=255)

                 tCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCg                                                                        >  1:2630699/1‑93 (MQ=255)

                  cTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCgg                                                                       <  2:662150/93‑1 (MQ=255)

                      aaCGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTc                                                                   <  1:3857138/93‑1 (MQ=255)

                      aaCGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTc                                                                   <  1:355415/93‑1 (MQ=255)

                      aaCGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTc                                                                   >  2:5738322/1‑93 (MQ=255)

                      aaCGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTc                                                                   <  2:6259273/93‑1 (MQ=255)

                      aaCGAGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTc                                                                   <  1:6258577/93‑1 (MQ=255)

                       aCGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCa                                                                  <  1:5039927/93‑1 (MQ=255)

                       aCGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCa                                                                  >  2:3018395/1‑93 (MQ=255)

                         gggAATAAACCATGGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGt                                                                <  2:5252270/93‑1 (MQ=255)

                          ggAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGTGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTg                                                               >  1:764460/1‑93 (MQ=255)

                          ggAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTg                                                               >  1:3126414/1‑93 (MQ=255)

                          ggAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTg                                                               <  1:2260620/93‑1 (MQ=255)

                          ggAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTg                                                               <  1:1535651/93‑1 (MQ=255)

                           gaagaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTgg                                                              >  1:5355692/1‑93 (MQ=255)

                           gaagaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTgg                                                              <  2:5327198/93‑1 (MQ=255)

                           gaagaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTgg                                                              <  2:1464163/93‑1 (MQ=255)

                            aagaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTggg                                                             >  1:2642600/1‑93 (MQ=255)

                            aagaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTggg                                                             >  2:4328936/1‑93 (MQ=255)

                              gaaACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGat                                                           >  1:4459961/1‑93 (MQ=255)

                               aaaCCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGatc                                                          >  2:4045678/1‑93 (MQ=255)

                               aaaCCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGatc                                                          <  1:2535902/93‑1 (MQ=255)

                                aaCCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGatca                                                         >  2:1594291/1‑93 (MQ=255)

                                 aCCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGatcat                                                        >  2:1748176/1‑93 (MQ=255)

                                  ccATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGatcatc                                                       <  2:1371128/93‑1 (MQ=255)

                                     ttGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACt                                                    >  2:262569/1‑93 (MQ=255)

                                      tGCGGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACtt                                                   <  1:3649411/93‑1 (MQ=255)

                                      tGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACtt                                                   <  1:5726147/93‑1 (MQ=255)

                                        cAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACtttt                                                 <  1:433588/93‑1 (MQ=255)

                                         aGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTg                                                <  2:1865572/93‑1 (MQ=255)

                                         aGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTg                                                <  2:4412493/93‑1 (MQ=255)

                                         aGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTg                                                >  2:1485309/1‑93 (MQ=255)

                                            cAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAgt                                             <  2:5624122/93‑1 (MQ=255)

                                             aaTCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAgtg                                            <  1:2694476/93‑1 (MQ=255)

                                             aaTCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAgtg                                            >  1:322711/1‑93 (MQ=255)

                                              aTCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAgtgt                                           >  1:5869027/1‑93 (MQ=255)

                                              aTCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAgtgt                                           >  2:53500/1‑93 (MQ=255)

                                               tCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTc                                          <  1:1238421/93‑1 (MQ=255)

                                               tCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTc                                          >  1:4193177/1‑93 (MQ=255)

                                               tCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTc                                          >  1:5224255/1‑93 (MQ=255)

                                               tCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTc                                          <  1:1283925/93‑1 (MQ=255)

                                                cTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTcc                                         >  2:606049/1‑93 (MQ=255)

                                                 tttttGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCt                                        >  2:2099966/1‑93 (MQ=255)

                                                 tttttGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCt                                        <  2:5572202/93‑1 (MQ=255)

                                                   tatGCTGATAGCGGTAAGCGCATTAAGGTCAGGAGAAAGAGCGTCGATATCCACTACCGGTTTCAGTGGGATCATAACTTTTGAGTGTCAtgt                                      <  1:6297801/91‑1 (MQ=255)

                                                    ttGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCtgtt                                     <  1:2126836/93‑1 (MQ=255)

                                                      gCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCtgttgt                                   >  1:4745870/1‑93 (MQ=255)

                                                       cTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTgttgtt                                  <  1:4200106/93‑1 (MQ=255)

                                                       cTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTgttgtt                                  <  1:2363857/93‑1 (MQ=255)

                                                         gATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGc                                >  1:2023584/1‑93 (MQ=255)

                                                           tAGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCa                              >  1:18574/1‑93 (MQ=255)

                                                            aGCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCAc                             >  1:2741990/1‑93 (MQ=255)

                                                             gCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCAcc                            <  2:3999896/93‑1 (MQ=255)

                                                             gCGGTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCAcc                            >  2:2460201/1‑93 (MQ=255)

                                                                gTCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGa                         >  1:5707027/1‑93 (MQ=255)

                                                                 tCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGaa                        >  1:5887377/1‑93 (MQ=255)

                                                                 tCAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGaa                        <  1:3029854/93‑1 (MQ=255)

                                                                  cAGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAAt                       >  1:6059286/1‑93 (MQ=255)

                                                                   aGCGCATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAAtt                      >  1:5582678/1‑93 (MQ=255)

                                                                    gcgcATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTg                     >  2:5294237/1‑93 (MQ=255)

                                                                    gcgcATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTg                     >  1:2995349/1‑93 (MQ=255)

                                                                    gcgcATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTACTGTTGTTGCCACCTGAATTg                     >  2:4731905/1‑93 (MQ=255)

                                                                      gcATTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGc                   <  2:2217813/93‑1 (MQ=255)

                                                                        aTTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCaa                 >  1:2226593/1‑93 (MQ=255)

                                                                        aTTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCaa                 <  2:5080563/93‑1 (MQ=255)

                                                                        aTTAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCaa                 >  1:6244855/1‑93 (MQ=255)

                                                                         ttAAGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAg                >  2:2879411/1‑93 (MQ=255)

                                                                           aaGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGcc              <  1:3356319/93‑1 (MQ=255)

                                                                            aGGTCAGGCGAAATAGCTTCGATATCCACGTCCGGTTTCAGTGGGATCATCGCTTTTCAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCt             <  2:1979011/93‑1 (MQ=255)

                                                                            aGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCt             <  2:1986936/93‑1 (MQ=255)

                                                                            aGGTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCt             <  1:3455210/93‑1 (MQ=255)

                                                                             ggTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTg            <  2:3729830/93‑1 (MQ=255)

                                                                              gTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTgg           <  1:4559350/93‑1 (MQ=255)

                                                                              gTCAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTgg           >  2:3181118/1‑93 (MQ=255)

                                                                                cAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCa         <  1:3693137/93‑1 (MQ=255)

                                                                                cAGGCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCa         <  1:4430415/93‑1 (MQ=255)

                                                                                 aGGCGAAAGCGCGTCGATATCCACTTCCGGTTTCAGAGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAg        >  2:6029892/1‑93 (MQ=255)

                                                                                  ggCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGa       >  1:2867741/1‑93 (MQ=255)

                                                                                   gCGAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGaa      >  1:884821/1‑93 (MQ=255)

                                                                                     gAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAata    >  1:1118864/1‑93 (MQ=255)

                                                                                     gAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAata    >  1:1836778/1‑93 (MQ=255)

                                                                                     gAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAata    <  1:3316678/93‑1 (MQ=255)

                                                                                     gAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAata    >  2:3682952/1‑93 (MQ=255)

                                                                                     gAAAGAGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAata    >  2:2315140/1‑93 (MQ=255)

GTTTTTGCCTGGGCGGATCTGGAACGGGAAGAAACCATTGCAGCCAATCTTTTTGCTGATAGCGGTCAGCGCATTAAGGTCAGGCGAAAGGGCGTCGATATCCACTTCCGGTTTCAGTGGGATCATCACTTTTGAGTGTCCTGTTGTTGCCACCTGAATTGGCAAGCCTGGCAGAATATC  >  NC_000913/1533265‑1533444

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