BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,413,507	(G)_5→4	coding (1904/2601 nt)	recE ←	exonuclease VIII

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,503	0	G	.	99.4%	715.8 / ‑3.2	158	coding (1908/2601 nt)	recE	exonuclease VIII
Reads supporting (aligned to +/- strand): ref base G (1/0); new base . (70/86); total (71/87)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 178 total aligned reads displayed.

ACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTAAT  >  NC_000913/1413411‑1413591

 cGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACgg                                                                                         <  2:2784390/93‑1 (MQ=255)

  gATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACggg                                                                                        <  2:3527263/93‑1 (MQ=255)

  gATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACggg                                                                                        >  1:2005822/1‑93 (MQ=255)

   aTAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACgggg                                                                                       <  1:1977298/93‑1 (MQ=255)

   aTAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACgggg                                                                                       <  2:1688105/93‑1 (MQ=255)

   aTAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACgggg                                                                                       >  1:4130615/1‑93 (MQ=255)

    tAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACggggc                                                                                      >  2:3437371/1‑92 (MQ=255)

     aaaGCGGTTACTGAATTCTTCCGGTTCAAGTACCCCGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACggggca                                                                                     >  2:4542250/1‑91 (MQ=255)

       aGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACggggcatt                                                                                   >  1:4347428/1‑89 (MQ=255)

        gCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCAttt                                                                                 >  2:89724/1‑93 (MQ=255)

         cGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCAtttt                                                                                <  2:1307697/93‑1 (MQ=255)

           gTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTa                                                                              <  1:5642270/93‑1 (MQ=255)

            ttACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTAc                                                                             <  1:1079138/93‑1 (MQ=255)

            ttACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTAc                                                                             <  2:1366869/93‑1 (MQ=255)

             tACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTTCTAAATCGAGCGTTTTTGTCGTGGTGGTGTCCAC‑GGGGCATTTTTACg                                                                            >  1:293616/1‑93 (MQ=255)

              aCTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGCGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGc                                                                           >  2:6188259/1‑93 (MQ=255)

              aCTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGc                                                                           >  2:1146884/1‑93 (MQ=255)

              aCTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCAGTTTTACGc                                                                           <  2:6122998/93‑1 (MQ=255)

               cTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGcc                                                                          >  1:5809338/1‑93 (MQ=255)

                   aTTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACaa                                                                      <  2:5525909/93‑1 (MQ=255)

                    ttcttcCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACaaa                                                                     <  2:2402220/93‑1 (MQ=255)

                    ttcttcCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACaaa                                                                     <  1:1204960/93‑1 (MQ=255)

                     tcttcCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAat                                                                    <  1:5004098/93‑1 (MQ=255)

                      cttcCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTACATCGAGCGTTTTTGTCTTTGTGGGGTCCAC‑GGGGCAGTTCTAGGCCACAAata                                                                   >  2:2999398/1‑93 (MQ=255)

                      cttcCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAata                                                                   <  1:840737/93‑1 (MQ=255)

                        tcCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAAtata                                                                 <  1:761500/93‑1 (MQ=255)

                         ccGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAg                                                                <  2:255002/93‑1 (MQ=255)

                          cGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGt                                                               <  1:1765624/93‑1 (MQ=255)

                          cGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGt                                                               <  2:536214/93‑1 (MQ=255)

                            gTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGc                                                             >  1:4227539/1‑93 (MQ=255)

                             ttCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGcc                                                            >  2:2683334/1‑93 (MQ=255)

                              tCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCg                                                           >  1:4367843/1‑93 (MQ=255)

                              tCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCg                                                           <  1:3097607/93‑1 (MQ=255)

                               cAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCgg                                                          <  1:105355/93‑1 (MQ=255)

                               cAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCgg                                                          <  1:2253629/93‑1 (MQ=255)

                                aaGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGa                                                         >  1:4943844/1‑93 (MQ=255)

                                   tACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTa                                                      <  2:5232755/93‑1 (MQ=255)

                                     cccGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATc                                                    <  1:5274141/93‑1 (MQ=255)

                                     cccGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATc                                                    <  2:2902371/93‑1 (MQ=255)

                                      ccGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCa                                                   >  1:4601043/1‑93 (MQ=255)

                                      ccGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCa                                                   <  2:4501474/93‑1 (MQ=255)

                                      ccGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCa                                                   >  1:3917732/1‑93 (MQ=255)

                                      ccGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCa                                                   <  2:6278538/93‑1 (MQ=255)

                                       cGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAg                                                  <  1:5695491/93‑1 (MQ=255)

                                          cAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCaa                                               <  2:1787554/93‑1 (MQ=255)

                                           aGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAAt                                              <  1:4240926/93‑1 (MQ=255)

                                            gTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATg                                             >  2:3923388/1‑93 (MQ=255)

                                            gTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATg                                             <  1:3049418/93‑1 (MQ=255)

                                             tGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGt                                            <  1:765067/93‑1 (MQ=255)

                                              ggAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTc                                           >  1:223049/1‑93 (MQ=255)

                                              ggAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTc                                           >  1:6200320/1‑93 (MQ=255)

                                               gAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCa                                          <  2:3110936/93‑1 (MQ=255)

                                               gAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCa                                          >  2:4491821/1‑93 (MQ=255)

                                               gAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCa                                          >  2:5982167/1‑93 (MQ=255)

                                               gAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCa                                          <  2:6144932/93‑1 (MQ=255)

                                                 aaGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATc                                        <  1:5412370/93‑1 (MQ=255)

                                                 aaGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATc                                        <  1:2107705/93‑1 (MQ=255)

                                                  aGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCg                                       >  2:488245/1‑93 (MQ=255)

                                                  aGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCg                                       >  1:3229785/1‑93 (MQ=255)

                                                   gCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGa                                      >  2:1225844/1‑93 (MQ=255)

                                                   gCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGa                                      <  2:5493956/93‑1 (MQ=255)

                                                   gCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCAAAAATATAGTGCCGGAGTATCAGCAATGTCATCGa                                      <  2:1549763/93‑1 (MQ=255)

                                                     aGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGc                                    >  2:3059375/1‑93 (MQ=255)

                                                      gTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCt                                   >  2:1457666/1‑93 (MQ=255)

                                                      gTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCt                                   >  2:526636/1‑93 (MQ=255)

                                                      gTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGGCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCt                                   >  1:3067593/1‑93 (MQ=255)

                                                          cTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTgaga                               <  1:5369181/93‑1 (MQ=255)

                                                           tAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGTGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGAc                              <  1:3742903/93‑1 (MQ=255)

                                                           tAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGAc                              >  2:4022147/1‑93 (MQ=255)

                                                           tAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGAc                              >  2:2122943/1‑93 (MQ=255)

                                                           tAAATCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGAc                              <  2:1067917/93‑1 (MQ=255)

                                                           tAAATCGAGCGTTTTTGTCTTTGTGGTGTCCA‑TGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGAc                              <  1:731416/93‑1 (MQ=255)

                                                              aTCGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTa                           <  2:4153028/93‑1 (MQ=255)

                                                                cGAGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACt                         <  2:4710202/93‑1 (MQ=255)

                                                                  aGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGa                       <  1:2937397/93‑1 (MQ=255)

                                                                  aGCGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGa                       >  2:2507793/1‑93 (MQ=255)

                                                                    cGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATa                     >  2:1462705/1‑93 (MQ=255)

                                                                    cGTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATa                     <  2:5604653/93‑1 (MQ=255)

                                                                     gTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATAc                    <  1:3979295/93‑1 (MQ=255)

                                                                     gTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATAc                    >  1:981095/1‑93 (MQ=255)

                                                                     gTTTTTGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATAc                    <  2:4286076/93‑1 (MQ=255)

                                                                      tttttGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATAcc                   <  2:5111442/93‑1 (MQ=255)

                                                                      tttttGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATAcc                   >  2:1645429/1‑93 (MQ=255)

                                                                      tttttGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATAcc                   >  2:1817188/1‑93 (MQ=255)

                                                                      tttttGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATAcc                   >  2:2327979/1‑93 (MQ=255)

                                                                         ttGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCggg                >  1:1281594/1‑93 (MQ=255)

                                                                         ttGTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCggg                <  1:5078799/93‑1 (MQ=255)

                                                                           gTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGAc              <  1:2320992/93‑1 (MQ=255)

                                                                           gTCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGAc              >  1:593696/1‑93 (MQ=255)

                                                                            tCTTTGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGAcc             >  2:5709935/1‑93 (MQ=255)

                                                                              tttGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGAcccc           >  1:2424638/1‑92 (MQ=255)

                                                                              tttGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCcg           >  1:6074271/1‑93 (MQ=255)

                                                                              tttGTGGTGTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCcg           <  1:327226/93‑1 (MQ=255)

                                                                                 gtggtgTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGt        <  2:2829828/93‑1 (MQ=255)

                                                                                 gtggtgTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGt        <  2:1510858/93‑1 (MQ=255)

                                                                                  tggtgTCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTg       >  1:5542889/1‑93 (MQ=255)

                                                                                    gtgtCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGt     <  1:5555158/93‑1 (MQ=255)

                                                                                      gtCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTaa   <  1:1406934/93‑1 (MQ=255)

                                                                                       tCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTAAt  >  1:393150/1‑93 (MQ=255)

                                                                                       tCCAC‑GGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTAAt  >  2:3382076/1‑93 (MQ=255)

ACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTAAT  >  NC_000913/1413411‑1413591

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