BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,398,840	+C	coding (964/1533 nt)	ydaH ←	putative pump protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,398,838	1	.	C	99.5%	634.4 / ‑5.7	198	G322G (GGC→GGG)	ydaH	putative pump protein
Reads supporting (aligned to +/- strand): ref base . (0/1); major base C (69/127); minor base T (1/0); total (70/128)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 217 total aligned reads displayed.

AAAACCATCACGATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATG‑CCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATT  >  NC_000913/1398746‑1398924

aaaaCCATCACGATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATg                                                                                         >  2:2688052/1‑93 (MQ=255)

    ccATCACGATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATG‑ccca                                                                                    >  2:1301465/1‑91 (MQ=255)

     cATCACGATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATG‑cccat                                                                                   <  1:5300637/93‑4 (MQ=255)

     cATCACGATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATG‑cccat                                                                                   <  1:3325758/93‑4 (MQ=255)

      aTCACGATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATG‑cccata                                                                                  >  2:3778140/1‑89 (MQ=255)

      aTCACGATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATG‑cccata                                                                                  >  1:3165896/1‑89 (MQ=255)

        cACGATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAg                                                                                 >  2:5672356/1‑93 (MQ=255)

         aCGATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGc                                                                                <  1:1300739/93‑1 (MQ=255)

         aCGATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGc                                                                                >  1:4006996/1‑93 (MQ=255)

          cGATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGcc                                                                               <  1:4156588/93‑1 (MQ=255)

           gATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCa                                                                              >  1:6077424/1‑93 (MQ=255)

           gATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCa                                                                              <  1:1694125/93‑1 (MQ=255)

            aTAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAg                                                                             <  2:3035166/93‑1 (MQ=255)

            aTAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAg                                                                             <  1:5393466/93‑1 (MQ=255)

              aaaTCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCg                                                                           <  1:4572915/93‑1 (MQ=255)

               aaTCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCga                                                                          <  1:2513091/93‑1 (MQ=255)

                 tCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCgaga                                                                        >  1:2158241/1‑93 (MQ=255)

                  cccGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGac                                                                       <  2:4872115/93‑1 (MQ=255)

                  cccGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGac                                                                       <  2:5327437/93‑1 (MQ=255)

                  cccGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGac                                                                       >  1:4251124/1‑93 (MQ=255)

                   ccGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGaca                                                                      <  2:472948/93‑1 (MQ=255)

                    cGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGacaa                                                                     <  1:3607194/93‑1 (MQ=255)

                    cGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGacaa                                                                     <  2:5536548/93‑1 (MQ=255)

                     gCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGacaac                                                                    <  2:2122738/93‑1 (MQ=255)

                     gCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGacaac                                                                    >  2:5105943/1‑93 (MQ=255)

                     gCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGacaac                                                                    <  2:4225688/93‑1 (MQ=255)

                      ccATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGacaaca                                                                   >  2:3961797/1‑93 (MQ=255)

                      ccATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGacaaca                                                                   >  2:662984/1‑93 (MQ=255)

                      ccATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGacaaca                                                                   <  1:2787629/93‑1 (MQ=255)

                        aTCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACaaa                                                                 <  2:1814151/93‑1 (MQ=255)

                        aTCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACaaa                                                                 >  2:1874715/1‑93 (MQ=255)

                          ctctTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGa                                                               <  2:381064/93‑1 (MQ=255)

                            ctTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGaaa                                                             >  2:3716981/1‑93 (MQ=255)

                             tttCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGaaaa                                                            >  2:1725280/1‑93 (MQ=255)

                               tCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGaaaaaa                                                          <  1:4815412/93‑1 (MQ=255)

                                cATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGaaaaaaa                                                         <  2:4937394/93‑1 (MQ=255)

                                 aTCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGaaaaaaaa                                                        >  2:1125590/1‑92 (MQ=255)

                                 aTCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAg                                                        <  1:3092369/93‑1 (MQ=255)

                                 aTCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAg                                                        <  1:234496/93‑1 (MQ=255)

                                  tCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGa                                                       <  1:6214225/93‑1 (MQ=255)

                                  tCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGa                                                       <  1:1084891/93‑1 (MQ=255)

                                  tCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGa                                                       >  1:2108664/1‑93 (MQ=255)

                                   cGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGaa                                                      <  2:3466752/93‑1 (MQ=255)

                                     gTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATg                                                    >  1:2327817/1‑93 (MQ=255)

                                     gTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATg                                                    <  2:1431699/93‑1 (MQ=255)

                                     gTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATg                                                    >  1:466791/1‑93 (MQ=255)

                                     gTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGCAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAAta                                                    >  1:5611424/1‑92 (MQ=255)

                                      ttCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGa                                                   >  1:5770511/1‑93 (MQ=255)

                                      ttCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGa                                                   <  1:5195393/93‑1 (MQ=255)

                                      ttCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGa                                                   <  2:3880621/93‑1 (MQ=255)

                                       tCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGATACAACAAAGAAAAAAAGAATGAt                                                  >  2:1907833/1‑93 (MQ=255)

                                        cAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATc                                                 <  1:4592739/93‑1 (MQ=255)

                                         aaTCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCa                                                <  1:4418839/93‑1 (MQ=255)

                                         aaTCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCa                                                >  1:4746147/1‑93 (MQ=255)

                                         aaTCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCa                                                <  1:5114313/93‑1 (MQ=255)

                                          aTCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAg                                               >  2:2393008/1‑93 (MQ=255)

                                           tCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGt                                              <  1:4697511/93‑1 (MQ=255)

                                           tCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGt                                              <  2:795964/93‑1 (MQ=255)

                                             aTTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTgg                                            >  2:4527003/1‑93 (MQ=255)

                                                  aTGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGa                                       <  2:4402537/93‑1 (MQ=255)

                                                   tGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGAt                                      >  1:3021885/1‑93 (MQ=255)

                                                   tGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGAt                                      <  1:311924/93‑1 (MQ=255)

                                                    gCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATa                                     <  1:5001242/93‑1 (MQ=255)

                                                     cGGTAAACCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATAc                                    <  2:4093058/93‑1 (MQ=255)

                                                      gttAAACCGGCTTGACGTCGAATGGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATAcc                                   <  1:5962616/91‑1 (MQ=255)

                                                      ggTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATAcc                                   <  1:218001/93‑1 (MQ=255)

                                                      ggTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATAcc                                   <  1:748251/93‑1 (MQ=255)

                                                       gTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCt                                  >  2:5141711/1‑93 (MQ=255)

                                                         aaaTCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCttt                                <  2:1794032/93‑1 (MQ=255)

                                                         aaaTCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGACACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCttt                                <  2:6111745/93‑1 (MQ=255)

                                                           aTCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTtaa                              <  2:5981515/93‑1 (MQ=255)

                                                              cGCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTtaataa                           >  2:2313604/1‑93 (MQ=255)

                                                               gCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATaaa                          <  2:5333522/93‑1 (MQ=255)

                                                               gCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATaaa                          >  2:4270670/1‑93 (MQ=255)

                                                               gCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATaaa                          <  2:817644/93‑1 (MQ=255)

                                                               gCCTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATaaa                          >  1:4294011/1‑93 (MQ=255)

                                                                ccTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAg                         >  2:1393664/1‑93 (MQ=255)

                                                                ccTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAg                         <  1:588033/93‑1 (MQ=255)

                                                                ccTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAg                         >  1:940870/1‑93 (MQ=255)

                                                                 cTGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAgg                        <  1:1021454/93‑1 (MQ=255)

                                                                  tGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAggg                       <  1:4985357/93‑1 (MQ=255)

                                                                  tGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAggg                       <  2:5877821/93‑1 (MQ=255)

                                                                  tGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAggg                       <  1:4388572/93‑1 (MQ=255)

                                                                  tGACGTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAggg                       <  1:956698/93‑1 (MQ=255)

                                                                      gTCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGAt                   <  1:4257102/93‑1 (MQ=255)

                                                                       tCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATg                  <  2:2546204/93‑1 (MQ=255)

                                                                       tCGAATTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATg                  <  2:5584765/93‑1 (MQ=255)

                                                                          aaTTGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCa               <  2:1374119/93‑1 (MQ=255)

                                                                            ttGTGCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATc             <  1:873092/93‑1 (MQ=255)

                                                                              gtgCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATACAGGGTGATGGCATCAc           >  1:3862055/1‑93 (MQ=255)

                                                                              gtgCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCAc           <  2:4797696/93‑1 (MQ=255)

                                                                               tgCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACg          <  1:4341219/93‑1 (MQ=255)

                                                                               ggCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACg          <  2:2798236/92‑1 (MQ=255)

                                                                                gCGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACgg         <  1:5816456/93‑1 (MQ=255)

                                                                                 cGGGTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGgt        <  1:1808307/93‑1 (MQ=255)

                                                                                  gggTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGgtg       >  2:1586459/1‑93 (MQ=255)

                                                                                  gggTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGgtg       <  1:1987278/93‑1 (MQ=255)

                                                                                   ggTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGgtgt      <  1:1413614/93‑1 (MQ=255)

                                                                                    gTAGCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGtgtg     <  2:3591621/93‑1 (MQ=255)

                                                                                       gCGATGCCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGAtt  >  1:3415634/1‑93 (MQ=255)

AAAACCATCACGATAAATCCCGCCATCTCTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATG‑CCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATT  >  NC_000913/1398746‑1398924

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