BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,672,460	G→A	H83H (CAC→CAT)	yphB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,672,460	0	G	A	86.5%	822.9 / 74.8	281	H83H (CAC→CAT)	yphB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (20/18); new base A (132/111); total (152/129)
Fisher's exact test for biased strand distribution p-value = 8.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

TAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGC  >  NC_000913/2672368‑2672549

tAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCg                                                                                           >  1:6415576/1‑93 (MQ=255)

tAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCg                                                                                           <  2:5977133/93‑1 (MQ=255)

  gTGATAGACCCCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCatg                                                                                         <  1:3962733/93‑1 (MQ=255)

   tGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCatgt                                                                                        <  1:5853415/93‑1 (MQ=255)

   tGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCatgt                                                                                        >  1:7097002/1‑93 (MQ=255)

    gATAGACACCGCTGCAATGTTCATACCCCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTa                                                                                       <  1:5904606/93‑1 (MQ=255)

       aGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGAt                                                                                    >  1:7069364/1‑93 (MQ=255)

       aGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGAt                                                                                    >  1:4923374/1‑93 (MQ=255)

       aGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGAt                                                                                    <  2:4168425/93‑1 (MQ=255)

        gACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATa                                                                                   >  2:3958743/1‑93 (MQ=255)

         acacCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAg                                                                                  <  2:1269229/93‑1 (MQ=255)

          cacCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGt                                                                                 >  2:5218310/1‑93 (MQ=255)

           acCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTg                                                                                >  2:7268572/1‑93 (MQ=255)

              gCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTgcgc                                                                             >  1:1806102/1‑93 (MQ=255)

              gCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTgcgc                                                                             >  1:4087710/1‑93 (MQ=255)

               cTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCa                                                                            <  2:2002358/93‑1 (MQ=255)

                tGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCAt                                                                           >  2:6686736/1‑93 (MQ=255)

                tGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCAt                                                                           >  2:2241403/1‑93 (MQ=255)

                tGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCAt                                                                           >  1:3042257/1‑93 (MQ=255)

                 gCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATc                                                                          >  2:1022110/1‑93 (MQ=255)

                 gCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATc                                                                          <  1:4730425/93‑1 (MQ=255)

                   gATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATccc                                                                        >  1:1161593/1‑93 (MQ=255)

                     tGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCAc                                                                      >  1:3381540/1‑93 (MQ=255)

                     tGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCAc                                                                      <  1:6254555/93‑1 (MQ=255)

                     tGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCAc                                                                      <  1:6350611/93‑1 (MQ=255)

                      gTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACt                                                                     >  2:4581368/1‑93 (MQ=255)

                      gTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACt                                                                     >  1:4958146/1‑93 (MQ=255)

                       ttCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTc                                                                    <  2:5756674/93‑1 (MQ=255)

                          aTACCCCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAAc                                                                 <  1:2550123/93‑1 (MQ=255)

                          aTACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAAc                                                                 <  1:4755245/93‑1 (MQ=255)

                           tACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACa                                                                >  2:257141/1‑93 (MQ=255)

                           tACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACa                                                                <  1:1061880/93‑1 (MQ=255)

                           tACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACa                                                                >  2:3141335/1‑93 (MQ=255)

                           tACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACa                                                                >  1:1507050/1‑93 (MQ=255)

                             cacCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACAtt                                                              >  1:1278513/1‑93 (MQ=255)

                              acCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCACTCAACATTc                                                             >  2:4432227/1‑93 (MQ=255)

                              acCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTc                                                             >  1:3876533/1‑93 (MQ=255)

                               ccAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCg                                                            <  2:953751/93‑1 (MQ=255)

                               ccAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCg                                                            <  2:2787577/93‑1 (MQ=255)

                                cAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCgg                                                           >  2:3851103/1‑93 (MQ=255)

                                 aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGt                                                          >  1:677952/1‑93 (MQ=255)

                                 aGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGt                                                          >  1:1442769/1‑93 (MQ=255)

                                  gACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGtt                                                         <  1:4480485/93‑1 (MQ=255)

                                   acacAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTg                                                        <  2:3570945/93‑1 (MQ=255)

                                    cacaAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGc                                                       <  2:5271816/93‑1 (MQ=255)

                                       aaaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGt                                                    <  1:6676603/93‑1 (MQ=255)

                                       aaaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGt                                                    >  1:5964505/1‑93 (MQ=255)

                                        aaCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGtt                                                   <  2:2177804/93‑1 (MQ=255)

                                         aCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTg                                                  >  1:2696101/1‑93 (MQ=255)

                                         aCTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTg                                                  >  2:6082595/1‑93 (MQ=255)

                                          cTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGa                                                 >  2:4323957/1‑93 (MQ=255)

                                            atcatcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATa                                               <  1:4198378/93‑1 (MQ=255)

                                             tcatcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATAc                                              <  2:301571/93‑1 (MQ=255)

                                               atcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTc                                            >  2:4256091/1‑93 (MQ=255)

                                                tcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCa                                           <  1:6191154/93‑1 (MQ=255)

                                                tcGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCa                                           >  2:168933/1‑93 (MQ=255)

                                                    tgtgCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGcc                                       >  1:989509/1‑93 (MQ=255)

                                                    tgtgCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGcc                                       >  1:249041/1‑93 (MQ=255)

                                                    tgtgCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGcc                                       >  1:5534860/1‑93 (MQ=255)

                                                       gCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCt                                    <  1:5555595/93‑1 (MQ=255)

                                                         gagaCACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGc                                  >  2:252044/1‑93 (MQ=255)

                                                          agaCACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATGCTCACGCCCCTGcc                                 >  2:4942352/1‑93 (MQ=255)

                                                          agaCACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGcc                                 >  1:6077979/1‑93 (MQ=255)

                                                          agaCACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGcc                                 <  1:3610992/93‑1 (MQ=255)

                                                           gaCACATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCa                                <  2:790837/93‑1 (MQ=255)

                                                            acacATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAg                               >  1:7328097/1‑93 (MQ=255)

                                                            acacATTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAg                               >  2:5420848/1‑93 (MQ=255)

                                                              acaTTGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGAc                             <  2:4973404/93‑1 (MQ=255)

                                                                 ttGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACaaa                          >  2:1517190/1‑93 (MQ=255)

                                                                 ttGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACaaa                          <  2:1823689/93‑1 (MQ=255)

                                                                  tGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACaaaa                         <  2:6538978/93‑1 (MQ=255)

                                                                  tGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACaaaa                         <  2:1155565/93‑1 (MQ=255)

                                                                  tGCCACTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACaaaa                         <  2:2127470/93‑1 (MQ=255)

                                                                    ccACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACg                       >  1:5199878/1‑93 (MQ=255)

                                                                      aCTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGAt                     <  1:7386869/93‑1 (MQ=255)

                                                                      aCTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGAt                     <  2:6790791/93‑1 (MQ=255)

                                                                       cTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGAtt                    <  1:5462470/93‑1 (MQ=255)

                                                                       cTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGAtt                    >  2:3171377/1‑93 (MQ=255)

                                                                       cTCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGAtt                    <  2:5455935/93‑1 (MQ=255)

                                                                        tCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTg                   >  1:2181683/1‑93 (MQ=255)

                                                                        tCGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTg                   >  2:6860982/1‑93 (MQ=255)

                                                                         cGCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGc                  <  2:6287607/93‑1 (MQ=255)

                                                                          gCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGcc                 >  1:3500952/1‑93 (MQ=255)

                                                                          gCCCAGCCAGCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGcc                 <  1:5702378/93‑1 (MQ=255)

                                                                            ccagccagCCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGc               <  1:2961189/93‑1 (MQ=255)

                                                                             cagccagcCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCt              >  1:3110221/1‑93 (MQ=255)

                                                                             cagccagcCATCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCt              >  2:3687162/1‑93 (MQ=255)

                                                                               gccagccaTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCa            >  1:830480/1‑93 (MQ=255)

                                                                               gccagccaTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCa            >  1:4350404/1‑93 (MQ=255)

                                                                                ccagccaTCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCAc           <  2:2135782/93‑1 (MQ=255)

                                                                                  agccaTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCAccc         <  1:7236865/93‑1 (MQ=255)

                                                                                  agccaTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCAccc         <  1:3284634/93‑1 (MQ=255)

                                                                                   gccaTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCg        >  1:536905/1‑93 (MQ=255)

                                                                                    ccaTCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCgg       >  1:4604263/1‑93 (MQ=255)

                                                                                     caTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGt      <  1:5271088/93‑1 (MQ=255)

                                                                                     caTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGt      <  2:3417526/93‑1 (MQ=255)

                                                                                     caTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGt      <  1:4555841/93‑1 (MQ=255)

                                                                                      aTCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGtt     >  2:7065743/1‑93 (MQ=255)

                                                                                       tCGCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGttt    >  2:7440387/1‑93 (MQ=255)

                                                                                         gCCATGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGc  >  1:2891304/1‑93 (MQ=255)

TAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCACTCAACATTCGGTTGCAGTTGATACTCACGCCCCTGCCAGACAAAACGATTGCCGCTCACCCGGTTTGC  >  NC_000913/2672368‑2672549

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