BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,542,050	A→G	L144P (CTG→CCG)	bioH ←	BioH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,542,050	0	A	G	86.4%	654.6 / 86.1	249	L144P (CTG→CCG)	bioH	BioH
Reads supporting (aligned to +/- strand): ref base A (19/15); new base G (104/111); total (123/126)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 259 total aligned reads displayed.

GTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCACG  >  NC_000913/3541959‑3542141

 tCACCCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCggg                                                                                           <  1:3347496/93‑1 (MQ=255)

  cAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGa                                                                                          >  1:5825606/1‑93 (MQ=255)

   aaCCTCCGGCATCGGTAACGCAAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACCCCAGGaa                                                                                         <  1:6537800/93‑1 (MQ=255)

    aCCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAAc                                                                                        <  2:3562267/93‑1 (MQ=255)

    aCCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAAc                                                                                        <  2:3686525/93‑1 (MQ=255)

       tCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGc                                                                                     >  1:126766/1‑93 (MQ=255)

       tCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGc                                                                                     >  2:2486150/1‑93 (MQ=255)

       tCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGc                                                                                     >  2:5046753/1‑93 (MQ=255)

       tCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGc                                                                                     >  2:4930691/1‑93 (MQ=255)

        ccGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCt                                                                                    <  1:1357202/93‑1 (MQ=255)

        ccGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCt                                                                                    >  1:1640994/1‑93 (MQ=255)

         cGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTc                                                                                   <  2:5940666/93‑1 (MQ=255)

            cATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCAc                                                                                <  2:2130311/93‑1 (MQ=255)

            cATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCAc                                                                                <  2:5796484/93‑1 (MQ=255)

            cATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCATGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCAc                                                                                >  2:3549978/1‑93 (MQ=255)

             aTCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACt                                                                               >  1:4347590/1‑93 (MQ=255)

             aTCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACt                                                                               >  2:4111661/1‑93 (MQ=255)

             aTCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCTCCACt                                                                               <  2:6615189/93‑1 (MQ=255)

               cGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGt                                                                             <  1:4752092/93‑1 (MQ=255)

               cGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGt                                                                             <  2:846768/93‑1 (MQ=255)

               cGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGt                                                                             <  2:698475/93‑1 (MQ=255)

                   aaCGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGc                                                                         <  1:6034851/93‑1 (MQ=255)

                    aCGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCt                                                                        <  2:1393453/93‑1 (MQ=255)

                    aCGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCt                                                                        <  1:2300297/93‑1 (MQ=255)

                    aCGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCt                                                                        <  1:4527037/93‑1 (MQ=255)

                    aCGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACACt                                                                        >  2:3365568/1‑93 (MQ=255)

                      gCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGa                                                                      <  1:6321129/93‑1 (MQ=255)

                      cccAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGa                                                                      <  1:4892593/92‑1 (MQ=255)

                        cAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGaaa                                                                    >  1:6167827/1‑93 (MQ=255)

                         aGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGACGTTTCAGTCCACATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGaaaa                                                                   >  1:5622127/1‑93 (MQ=255)

                          gAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAat                                                                  >  2:3915872/1‑93 (MQ=255)

                           aaCGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAatc                                                                 <  1:1477279/93‑1 (MQ=255)

                             cGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAatcat                                                               <  2:6193946/93‑1 (MQ=255)

                             cGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAatcat                                                               <  2:5927126/93‑1 (MQ=255)

                              ggTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAatcatc                                                              <  2:3991560/93‑1 (MQ=255)

                              ggTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAatcatc                                                              <  1:54312/93‑1 (MQ=255)

                               gTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAAtcatca                                                             <  1:4648498/93‑1 (MQ=255)

                                ttttCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCAc                                                            >  2:190824/1‑93 (MQ=255)

                                ttttCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCAc                                                            <  2:2364687/93‑1 (MQ=255)

                                ttttCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCAc                                                            <  2:2094881/93‑1 (MQ=255)

                                ttttCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCAc                                                            <  1:499183/93‑1 (MQ=255)

                                   tcttcAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGa                                                         <  2:4789809/93‑1 (MQ=255)

                                     ttcAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGt                                                       <  1:3533286/93‑1 (MQ=255)

                                      tcAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGtt                                                      >  2:3041112/1‑93 (MQ=255)

                                      tcAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGtt                                                      <  2:4365460/93‑1 (MQ=255)

                                       cAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTtg                                                     >  2:5864110/1‑93 (MQ=255)

                                        aaCGCCCGCGCATCCTGGCGCGCTGTTTCAGTCCCCATGGTTTGTAGCGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTtgc                                                    >  1:5340303/1‑93 (MQ=255)

                                         aCGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTtgct                                                   <  2:5502805/93‑1 (MQ=255)

                                         aCGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTtgct                                                   >  1:4514802/1‑93 (MQ=255)

                                          cGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTtgctg                                                  >  2:1684756/1‑93 (MQ=255)

                                           gCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTtgctgc                                                 <  2:1086692/93‑1 (MQ=255)

                                            cccGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTgctgct                                                <  2:2028231/93‑1 (MQ=255)

                                             ccGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGctgctg                                               >  2:293288/1‑93 (MQ=255)

                                                 gcATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAAt                                           >  2:4889420/1‑93 (MQ=255)

                                                 gcATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAAt                                           >  1:3785032/1‑93 (MQ=255)

                                                 gcATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAAt                                           >  1:5514117/1‑93 (MQ=255)

                                                  cATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGACAATCATCACTGAGTTGCTGCTGAAATc                                          >  2:1225228/1‑93 (MQ=255)

                                                  cATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGGACGCTGAAAATCATCACTGAGTTGCTGCTGAAATc                                          >  1:935810/1‑93 (MQ=255)

                                                   aTCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATcc                                         <  2:5003792/93‑1 (MQ=255)

                                                    tCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATccc                                        >  2:252065/1‑93 (MQ=255)

                                                    tCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGATGAAATccc                                        >  2:5683659/1‑93 (MQ=255)

                                                      cTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTCCGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGc                                      <  1:2354000/93‑1 (MQ=255)

                                                       tGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGcc                                     <  1:1062500/93‑1 (MQ=255)

                                                       tGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGcc                                     <  2:6567151/93‑1 (MQ=255)

                                                         gcgcgcCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAg                                   <  1:2023425/93‑1 (MQ=255)

                                                         gcgcgcCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAg                                   <  1:4114534/93‑1 (MQ=255)

                                                          cgcgcCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGc                                  <  2:2067309/93‑1 (MQ=255)

                                                           gcgcCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCGCTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCa                                 <  2:3396902/93‑1 (MQ=255)

                                                           gcgcCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCATTGAGTTGCTGCTGAAATCCCGccacca                                 >  1:203137/1‑93 (MQ=255)

                                                             gcCGTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACg                               >  2:6489865/1‑93 (MQ=255)

                                                             gcCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACg                               >  1:1877793/1‑93 (MQ=255)

                                                              ccGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGt                              >  1:783393/1‑93 (MQ=255)

                                                              ccGTTTCAGTCCCCATGGCTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGCGTTGCTGCTGACATCCCGCCAGCACGt                              >  1:4203882/1‑93 (MQ=255)

                                                                gTTTCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTcc                            <  1:6510514/93‑1 (MQ=255)

                                                                 tttCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCg                           >  1:6206251/1‑93 (MQ=255)

                                                                 tttCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCg                           >  1:2216813/1‑93 (MQ=255)

                                                                  ttCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCgg                          <  1:4919129/93‑1 (MQ=255)

                                                                  ttCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAACCCCGCCAGCACGTCCgg                          <  1:889518/93‑1 (MQ=255)

                                                                   tCAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGt                         <  1:1576128/93‑1 (MQ=255)

                                                                    cAGTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGtt                        <  1:3268682/93‑1 (MQ=255)

                                                                      gTCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGtttt                      >  2:459815/1‑93 (MQ=255)

                                                                       tCCCCATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTa                     <  1:5698370/93‑1 (MQ=255)

                                                                        ccccATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTAt                    >  2:1838444/1‑93 (MQ=255)

                                                                        ccccATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTAt                    >  2:4323356/1‑93 (MQ=255)

                                                                        ccccATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTAt                    <  1:2587143/93‑1 (MQ=255)

                                                                        ccccATGGTTTGTAACGCCGGAAGCCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTAt                    <  1:6104551/93‑1 (MQ=255)

                                                                           cATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATctc                 >  1:3851914/1‑91 (MQ=255)

                                                                           cATGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATccc                 >  2:5378820/1‑93 (MQ=255)

                                                                            aTGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATcccc                <  1:2724594/93‑1 (MQ=255)

                                                                            aTGGTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATcccc                <  1:1922155/93‑1 (MQ=255)

                                                                              ggTTTGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCgg              >  1:2510386/1‑93 (MQ=255)

                                                                              ggTTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCgg              >  1:1176419/1‑93 (MQ=255)

                                                                                tttGTAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGcc            >  2:3221569/1‑93 (MQ=255)

                                                                                    tAACGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTc        >  2:1048992/1‑93 (MQ=255)

                                                                                      aCGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGt      <  1:3503531/93‑1 (MQ=255)

                                                                                      aCGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGt      <  2:447388/93‑1 (MQ=255)

                                                                                      aCGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGt      <  2:4176280/93‑1 (MQ=255)

                                                                                       cGCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTc     >  1:5195393/1‑93 (MQ=255)

                                                                                        gCCGGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCa    <  2:5143003/93‑1 (MQ=255)

                                                                                          cAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCACg  <  1:4822620/93‑1 (MQ=255)

GTCAACCTCCGGCATCGGTAACGCCAGAACGGTTTTCTTCAACGCCCGCGCATCCTGGCGCGCCGTTTCAGTCCCCATGGTTTGTAACGCCAGGAACCGCTCCACTGTACGCTGAAAATCATCACTGAGTTGCTGCTGAAATCCCGCCAGCACGTCCGGTTTTATCCCCGGCCACTCGTCACG  >  NC_000913/3541959‑3542141

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